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A.  Report  Summary 

The  primary  objective  of  this  project  was  to  investigate  the  feasibility  of  realizing  a  large- 
scale,  Type-II  quantum  computer  using  NV -Diamond  (NVD).  The  resulting  device  was 
then  to  be  used  to  demonstrate  Type-II  quantum  computing.  The  key  goal  is  to  use  such 
a  quantum  computer  for  efficient  simulation  of  turbulent  flows,  via  the  quantum  lattice 
gas  dynamics  model.  In  the  proposed  approach,  the  qubits  are  spins  of  color  centers  in 
NVD,  located  in  a  single,  diffraction  limited  spot,  distinguished  from  one  another  via 
spectral  signatures.  The  inter-qubit  coupling  is  mediated  via  a  combination  of  external 
lasers  and  single-photons  in  photonic-band-gap  (PBG)  based  cavities.  The  cavities  are 
etched  directly  on  to  the  diamond  substrate.  The  set  of  coupled  qubits  in  each  spot 
represents  an  individual  quantum  computer  (QC),  with  a  large  number  of  QC’s  located  on 
a  single  substrate. 

The  key  technical  challenge  in  this  project  was  to  realize  a  high-quality  PBG 
structure  in  a  crystalline  thin  film  of  diamond.  This  task  proved  to  be  extremely  difficult, 
due  to  severe  technical  constraints,  given  in  part  to  the  limited  time  and  resources 
allocated  to  this  project.  Our  initial  approach  was  to  perform  lithography  on  a  single 
crystal  diamond.  However,  the  orientation  of  the  crystal  made  it  essentially  inert  against 
the  range  of  etching  chemicals  we  tried.  In  principle,  it  would  be  possible  to  make  use  of 
a  single  crystal  diamond  with  a  different  orientation  that  would  be  etchable.  However, 
another  constraint  of  the  type  of  PBG  structure  needed  for  this  project  is  that  the  substrate 
must  be  very  thin,  ideally  close  to  2  microns,  backed  by  a  silicon  substrate.  Single  crystal 
diamonds  this  thin  are  not  available.  A  possible  alternative  is  to  polish  a  thicker  crystal 
down  to  this  thickness.  We  were  unable  to  find  any  vendor  that  would  be  able  to  perform 
this  polishing.  Harris  offered  to  make  a  film  with  a  thickness  down  to  18  microns, 
without  guaranteeing  good  surface  quality,  and  stated  that  they  do  not  know  how  to  make 
it  any  thinner.  We  then  moved  to  using  the  so-called  UNCD  (Ultrananocrystalline 
Diamond)  film,  invented  and  developed  by  Krauss  et  al.  at  Argonne  National  Laboratory. 
Our  E-beam  lithographic  process,  after  some  mdoification  in  the  steps,  was  successful  in 
etching  the  requisite  patterns.  However,  the  holes  produced  turned  out  to  be  of  non- 
uniform  depths.  This  was  attributable  primarily  to  the  roughness  of  the  surface  of  the 
UNCD  film. 

We  then  collaborated  with  Steven  Prawer  of  Melbourne  University  in  Australia  to 
pursue  an  alternative  approach:  using  FIB  (focused  ion  beam)  etching  of  a  single  crystal. 
This  approach  was  useful  in  producing  a  PBG  pattern.  However,  the  FIB  process 
produced  microscopic  damages  that  severly  degraded  the  quality  of  the  pattern,  making 
this  process  unsuitable  for  a  high-Q  PBG  cavity  of  the  type  needed  for  our  project. 
Annealing  after  the  FIB  etching  process  improved  the  quality  of  the  PBG  pattern 
somewhat,  but  not  enough. 

In  parallel  to  our  effort  for  producing  a  PBG  structure  in  diamond,  we  pursued 
several  other  lines  of  work,  in  order  to  develop  complimentary  as  well  as  alternative 
technologies  for  realizing  a  quantum  computer.  These  include  (a)  development  of  a  new 
technology  of  arbitrary  pattern  lithography,  using  Bose-condensate  based  atomic 
interferometry,  (b)  developmen  of  a  quantum  computing  architecture  using  atomic 
ensembles  as  quantum  bits,  (c)  demonstrataion  of  ultra-low  power  non-linear  optic  effects 
via  coupliing  between  the  evanescent  field  of  a  tapered  nanofiber  and  surrounding  atoms, 
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as  an  enabling  technlogy  for  optically  linking  ensemble-based  quantum  bits,  and  (d) 
demonstration  of  measurement  of  absolute  phase  of  a  localized  field  via  the  Bloch- 
Siegert  oscillation,  as  a  technique  for  suppressing  driver  phase  correlated  fluctuations  in 
the  rotation  of  a  strongly  driven  quantum  bit. 

It  is  our  conclusion  that  significant  technical  hurdles  in  the  context  of  lithography 
have  to  be  overcome  before  a  PBG  based  QC  can  be  realized  in  NVD.  The  most 
promising  approach  would  be  to  use  the  polished  single  crystal  and  conventional 
lithography.  This  would  require  research  into  developing  techniques  for  polishing 
diamond  while  maintaining  surface  smoothness,  and  incorporating  this  film  into  the  E- 
beam  lithography  process.  In  figure  1 1  of  section  B.l,  we  offer  an  explicit  plan  for  how 
such  a  scheme  could  be  carried  out  in  the  future.  We  also  believe  that,  the  alternative 
technology  of  atomic  ensembles  as  quantum  bits,  linked  with  tapered  nanofibers,  as 
developed  under  this  grant,  is  a  more  promising  approach  for  realizing  a  quantum 
computer. 

Full  or  partial  support  from  this  grant  led  to  two  Ph.D.  theses,  and  a  total  of  21 
publications,  including  5  journal  papers. 
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B.  TECHNICAL  DETAILS 


B.l.  Lithographic  Effort  for  Producing  a  PBG  Structure  in  a  Diamond  Eilm 


The  concept  of  quantum  computing  using  NV-diamond  and  a  PBG  structure  is  recalled 
briefly  in  figure  1 .  The  left-most  diagram  shows  the  global  view  of  the  whole  substrate, 
containing  an  array  of  PBG  cavities,  indicated  by  red  circles.  Each  red-circled-cavity, 


Figure  1:  Schematic  illustration  of  the  concept  of  quantum  computing  using  NV-diamond  and  a  PBG 
structure.  See  text  for  details. 

along  with  the  active  molecules  (i.e.,  diamond  color  centers)  it  houses,  represent  a  mini 
quantum  computer  (QC).  Thus,  the  diagram  shows  16  mini  QC’s  on  the  same  substrate. 
Of  course,  a  larger  substrate  (say  12  cm  diameter)  could  host  a  much  larger  number  of 
mini  QC’s.  The  next  diagram  shows  a  zoomed-in  view  of  a  single  mi  ni-QC,  with  the  red 
circle  (top  view  of  a  donut)  showing  the  volume  occupied  by  the  cavity  mode.  The  third 
diagram  shows  schematically  tha  active  molecules  inside  the  cavity  mode. 

The  fourth  diagram  shows  all  of  these  molecules,  but  arranged  according  to  their 
spectral  response.  Because  of  local  non-uniformities,  otherwise  identical  molecules  have 
a  range  of  resonance  frequencies  for  any  given  transition  —  the  so-called  inhomogeneous 
broadening.  Our  model  couples  a  molecule  at  one  spectral  point  to  its  spectral  neighbor; 
the  coupling  is  mediated  by  the  cavity  photon,  and  controlled  by  two  external  laser 
frequencies.  Each  mini-computer  can  have  a  large  number  of  distinct  quantum  bits; 
however,  our  goal  was  to  work  with  a  more  modest  number  (about  10  to  start  with, 
eventually  reaching  a  100,  in  each  mini-QC).  Thus,  the  whole  chip  would  consist  of  say 
about  10000  mini-QC’s,  each  of  which  will  hold  about  100  qubits.  Inside  each  mini-QC, 
the  coupling  is  quantum  mechanical.  Among  different  mini-QC’s,  the 
communication/coupling  would  be  strictly  classical.  This  is,  of  course,  the  so-called 
Type  II  Quantum  Computer. 

The  first  step  in  our  project  was  to  design  the  PBG  cavity  structure,  using  FTDT 
simulation.  Figure  2  illustrates  the  typical  process  we  followed  in  performing  the 
simulation. 

Typically,  NV-Diamond  is  realized  by  electron  bombardment  and  annealing  in  a 
nitrogen  atmosphere.  Several  different  types  of  diamond  in  principle  can  be  used  in  this 
process.  These  include: 
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•  Type  IA  Diamond:  This  contains  nitrogen  as  an  impurity  in  fairly  substantial 
amounts  (of  the  order  of  0.1%),  and  which  appears  to  have  segregated  into  small 
aggregates.  Also  contains  platelets,  associated  with  the  nitrogen  impurity,  the 
exact  structure  of  which  is  not  known.  It  has  an  absorption  edge  at  ~300nm  and  a 
broad  absorption  band  between  7  and  10pm.  Most  natural  diamonds  are  of  this 
type. 

•  Type  IB  diamond:  This  also  contains  nitrogen  as  an  impurity  but  in  dispersed 
substitutional  form.  Almost  all  synthetic  diamonds  are  of  this  type. 

•  Type  IIA  Diamond:  This  is  effectively  free  of  nitrogen  impurity.  Very  rare  in 
nature,  these  diamonds  have  enhanced  optical  and  thermal  properties.  The  high 
purity  results  in  an  ultraviolet  transmission  band  down  to  approximate  230pm  and 
absence  of  infrared  absorption  in  the  7- 10pm  band. 

•  Type  IIB  Diamond:  A  very  pure  type  of  diamond  which  has  semiconducting 
properties:  generally  blue  in  color.  Extremely  rare  in  nature. 

•  CVD  Grown  Diamond:  These  are  diamond  films  grown  on  a  substrate  via 
chemical  vapor  deposition.  These  films  can  be  made  very  thick,  and  are  used  for 
a  wide  range  of  industrial  applications. 

•  UNCD(Ultra  Nano-Crystalline  Diamond):  UNCD  is  the  name  for  thin-film 
diamond  material  grown  through  a  patented  plasma-enhanced  chemical  vapor 
deposition  process.  UNCD  is  not  a  diamond-like  carbon  material,  but  crystalline 
diamond  consisting  of  nano-sized  grains. 

We  obtained  several  different  types  of  diamond  for  our  project.  These  include: 

•  CVD  Diamond  Wafer:  3pm  thick  diamond  on  a  3”  dia  x  0.060”  thick  Si 
substrate.  Vendor:  Delaware  Diamond  Knives  (DDK),  Inc. 

•  Bulk  Diamond  Ciystal:  2  mm  diameter  X  0.25  mm  thick  [110  cut]  diamond  disk 
cut  from  bulk  diamond  crystal.  Vendor:  Drukker  International 

•  Crsytalline  Optical  Grade  Diamond :  0. 15mm  thick  x5.0  mm  diameter  disc. 
Vendor:  Harris  Int. 

•  Polycrystalline  Optical  Grade  Diamond:  0.25mm  thick  x5.0  mm  diameter  disc. 
Vendor:  Harris  Int. 

•  UNCD  Diamond  Films:  UNCD  film  on  Si  wafer  of  several  different  thicknesses 
(from  0.5  pm  to  5  pm)  and  diameters.  Supplier:  the  group  of  Dr.  J.  Carlisle’s  at 
Argonne  National  Lab. 

Our  initial  approach  was  to  perform  lithography  on  a  single  crystal  diamond  (the  bulk 
diamond  crystal  from  Drukker  International).  However,  the  orientation  of  the  crystal 
made  it  essentially  inert  against  the  range  of  etching  chemicals  we  tried.  In  principle,  it 
would  be  possible  to  make  use  of  a  single  crystal  diamond  with  a  different  orientation 
that  would  be  etchable.  However,  another  constraint  of  the  type  of  PBG  structure  needed 
for  this  project  is  that  the  substrate  must  be  very  thin,  ideally  close  to  2  microns,  backed 
by  a  silicon  substrate.  Single  crystal  diamonds  this  thin  are  not  available.  A  possible 
alternative  is  to  polish  a  thicker  crystal  down  to  this  thickness.  We  were  unable  to  find 
any  vendor  that  would  be  able  to  perfonn  this  polishing.  Harris  offered  to  make  a  film 
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with  a  thickness  down  to  18  microns,  without  guaranteeing  good  surface  quality,  and 
stated  that  they  do  not  know  how  to  make  it  any  thinner. 

Given  this  impasse,  we  proceeded  to  explore  the  other  types  of  diamond  that  were 
available  as  films  deposited  on  a  substrate.  First,  we  tried  the  CVD  films  obtained  from 


Figure  2:  A  cantilever  structure  produced  in  a  UNCD  film. 

DDK  Inc.  However,  the  surface  of  this  crystal  was  very  rough,  and  clearly  unsuited  for 
lithography.  At  this  point,  we  were  made  aware  of  the  use  of  UNCD  diamond  crystals 
that  were  used  by  our  colleagues  in  the  Mechanical  Engineering  department  for  creating 
cantilevers,  as  shown  in  figure  2.  As  such,  we  decided  to  pursue  the  use  of  this  material 
for  creating  PBG  structure. 

Figure  3  summarizes  the  first  result  we  obtained  with  UNCD.  Figure  3a  shows  the 
photolithographic  steps  used  for  transferring  a  computer  generated  mask  to  the  diamond 
substrate  covered  by  Silicon  Nitride.  Figure  3b  shows  the  actual  wafer,  which  consists  of 
a  silicon  dioxide  base  under  the  diamond,  and  a  Silicon  Nitride  film  above  it,  produced 
by  CVD.  This  choice  wa  made  after  several  different  trials  with  other  types  of  layers. 
Figure  3c  shows  the  details  of  the  patterns  etched  in  diamond.  Of  course,  the  feature  size 
here  is  rather  large,  due  to  the  fact  that  we  used  photolithography  for  the  pattern  transfer. 

The  next  step  was  to  repeat  this  process,  but  with  much  smaller  feature  size, 
employing  E-Beam  techniques  instead  of  photolithography.  Before  doing  this,  we 
embarked  first  on  designing  an  actual  band-gap  structure,  tailored  specifically  to  the 
index  of  this  diamond  film.  Figure  4  summarizes  our  two-dimensional  modeling  of  PBG 


6 


SAMPLE  TYPICAL  PATTERNS  MICROSCOPE  IMAGE  DPN  IMAGE  OF  EDGES 


Figure  3:  Illustration  of  the  lithographic  processes  performed  for  producing  pattern  in  diamond.  [A] 
photolithography  used  to  transfer  mask  pattern.  [B]  the  layers  used  for  performing  the  lithography 
sequentially.  [C]  observation  of  patterns  produced  in  diamond  film. 


Figure  4:  FDTD  simulation  for  photonic  bandgap  structure,  performed  for  GHz-scale  signals.  This  is 
then  used,  through  scaling  laws,  to  develop  the  corresponding  structure  for  diamond  [A]  Two-dimensional 
structure  with  plane  wave  input  [B]  The  calculated  bandgap.  [C]  typical  profiles  of  input  fields  and  field 
propagating  into  the  structure. 


structure.  It  is  well  known  that  a  PBG  structure  can  be  designed  first  at  any  wavelength, 
and  for  any  index  variation,  and  then  scaling  rules  can  be  used  to  transfer  the  design  to 
any  other  wavelength  range  and  indices.  As  such,  we  chose  to  perform  the  simulation  in 
the  GHZ  range.  This  makes  it  possible  to  compare  easily  the  predictions  of  our  model  to 
existing  calculation,  without  concerning  ourselves  with  the  scaling  issue  up  front.  The 
result  shown  in  figure  4c  is  in  close  agreement  with  known  calculations  of  this  type. 
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Figure  5  shows  the  parameters  for  the  corresponding  PBG  structure  in  NVD.  The 
larger  hole  added  in  the  center  leads  to  a  cavity  structure.  We  have  also  calculated  the 
mode  pattern  of  several  cavity  designs,  as  summarized  in  figure  6. 
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Figure  5:  Design  of  a  photonic  bandgap  structure  for  diamond,  by  sealing  the  parameters  from  figure  6. 
The  larger  whole  inserted  in  the  middle  leads  to  a  cavity.  We  have  also  calculated  (not  shown)  the  cavity 
field  mode  pattern. 


3D  View  of  Field 
Over  the  Center  and 
Its  Neighborhood 


3D  View  of  Field  Over  the  Center 
and  Its  Neighborhood 


Field  Concentration  at  the 
Defect  Location 


Figure  6:  The  field  concentrations  produced  in  several  alternate  designs  of  PBG  cavities  in  diamond. 


Figure  7  illustrates  schematically  the  E-beam  lithographic  process  we  used  to 
produce  very  small  patterns  in  an  UNCD  film.  The  substrate  produced  here  is  similar  to 
the  one  shown  in  figure  3b.,  except  that  the  photoresist  is  replaced  by  a  layer  of  PMMA, 
and  the  SiCF  layer  is  replaced  by  Al,  in  order  to  allow  the  proper  etch  stop. 
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1.  Deposit  UNCD  film 


5.  Development  on  PMMA 


Diamond  film  :  2  um 

Si-substrate 

2.  Deposit  Al  metal  6  4/  Wet  Chemical  Etching 

A I  :  0.2  um 


3.  Deposit  PMMA  on  Al  ?  ft/E  Etching  on  Diamond  film 


PMMA  :  0.5ii  in 

4.  E-beam  Lithography  8.  Removal  of  Al  mask 


Figure  7:  Schematic  illustration  of  the  E-beam  lithographic  process  we  used  to 
produce  very  small  patterns  in  an  UNCD  film.  The  substrate  produced  here  is 
similar  to  the  one  shown  in  figure  3b.,  except  that  the  photoresist  is  replaced  by  a 
layer  of  PMMA,  and  the  Si02  layer  is  replaced  by  Al,  in  order  to  allow  the  proper 
etch  stop. 


Figure  8  shows  a  typical  pattern  produced  on  a  PMMA  substrate.  Briefly,  we 
placed  a  thin  film  of  PMMA,  using  spin  coating,  on  top  of  the  type  of  wafer  shown  in 
figure  7,  step  2.  We  then  used  an  E-beam  (with  a  resolution  of  20  nm)  to  transfer  the 
desired  pattern  to  the  PMMA.  Here,  we  have  shown  an  AFM  picture  of  the  PMMA 
layer. 
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Figure  8:  20  nm  resolution  E-beam  lithography  was  used  to  produce  this  pattern  (AFM  picture)  in  a 
PMMA  film  produced  on  top  of  a  Silicon  Nitride  film,  which  in  turn  is  above  a  diamond  film.  The  next  step 
is  to  use  this  PMMA  pattern  to  etch  the  nitride  film,  and  then  etch  the  diamond  film. 


Figure  9  shows  in  greater  details  the  AFM  data  obtained  from  this  structure.  As  can  be 
seen,  the  patterns  are  not  very  smooth  at  the  bottom.  This  is  attributable  primarily  to  the 
surface  quality  of  the  underlying  UNCD  film.  When  we  attempted  to  transfer  this  pattern 
to  the  Aluminum  and  then  to  the  diamond,  this  non-uniformity  produced  very  uneven 
patterns,  , making  it  obvious  that  it  would  be  very  difficult  using  this  approach  to  produce 
the  high-quality  PBG  cavity  in  the  UNCD  film. 


Figure  9:  More  details  of  the  AFM  data  obtained  from  this  structure.  As  can  be  seen,  the  patterns  are  not 
very  smooth  at  the  bottom.  This  is  attributable  primarily  to  the  surface  quality  of  the  underlying  UNCD 
film. 


At  this  point,  we  decided  that  we  must  go  back  to  the  idea  of  etching  single 
crystal  diamond.  However,  as  described  above,  we  were  unable  to  find  any  vendor 
willing  or  able  to  polish  a  crystal  down  to  few  microns  while  maintaining  the  optical 
quality  surface.  As  such,  we  decided  to  collaborate  with  the  group  of  Steven  Prawer  at 
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University  of  Melbourne  in  Australia.  His  group  was  pursuing  the  use  of  focused  ion 
beam  (FIB)  etching  to  produce  micro-structure  in  a  single  crystal  diamond.  Figure  10 
shows  some  typical  structures  produced  by  his  group  using  the  FIB  process. 


Cantilever  Waveguide 


Figure  10.  Cantilever  and  waveguide  structures  produced  by  FIB  etching  of  a  single  crystal  diamond 
by  the  group  of  Steven  Prawer  at  Melbourne  University. 


They  also  produced  a  PBG  pattern;  however,  the  FIB  process  produced 
microscopic  damages  that  severly  degraded  the  quality  of  the  pattern,  making  this  process 
unsuitable  for  a  high-Q  PBG  cavity  of  the  type  needed  for  our  project.  Annealing  after 
the  FIB  etching  process  improved  the  quality  of  the  PBG  pattern  somewhat,  but  not 
enough.  It  is  plausible  that  further  refinement  of  the  FIB  process  as  well  as  annealing 
may  suppress  the  damages  to  a  level  that  would  be  adequate  for  a  high  quality  PBG 
cavity. 


Given  our  findings,  we  envision  only  one  possible  route  for  realizing  high-Q  PBG 
cavities  suitable  for  NV-Diamond  quantum  computing.  This  idea  is  summarized  in 
figure  1 1 .  Briefly,  we  will  start  with  a  polished,  1 8  micron  thick  single  crystal  diamond 
film  on  a  silicon  substrate,  to  be  obtained  from  Harris  International,  with  an  etchable 
crystal  orientation.  We  will  then  follow  the  process  of  E-beam  lithoography  summarized 
above  to  produce  arrays  of  PBG  cavity  structures  on  the  top  surface  of  this  film  (steps  1 
through  3).  The  process  will  then  be  repeated  on  the  other  side,  thus  ending  up  with  a 
thin  (2  microns)  layer  of  diamond  with  the  PBG  cavities  (steps  4  through  6),  backed  by  a 
silicon  substrate  on  the  top.  Finally,  electron  bombardment  would  be  used  to  produce  the 
color  centers  in  order  to  convert  the  diamond  into  NVD. 

This  approach  is  promising.  However,  the  quality  of  the  structure  that  can  be 
realized  this  way  can  only  be  detennined  through  experimental  studies.  We  had 
exhausted  the  funds  and  the  time  alloted  under  this  project  by  the  time  we  reached  this 
conclusion.  As  such,  this  task  can  only  be  carried  out  with  additional  support  from 
AFOSR  or  other  funding  agencies. 
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STEP  1 


STEP  4 


PMMA 

AL 

18  micron  Single  Crystal 
Diamond  (SCD) 

Si-Substrate 


STEP  2 


After  Etching  the  PBG 
patternin  SCD 


STEP  6 


After  Etching  away  a  large 
area  of  the  Si-Substrate 


Figure  11:  Schematic  illustration  for  creating  PBG  cavities  with  E-beam  lithography  and  a  single 
crystal  diamond  film.  Briefly,  we  will  start  with  a  polished,  18  micron  thick  single  crystal  diamond  film 
on  a  silicon  substrate,  to  be  obtained  from  Harris  International,  with  an  etchable  crystal  orientation. 
We  will  then  follow  the  process  of  E-beam  lithoography  summarized  above  to  produce  arrays  of  PBG 
cavity  structures  on  the  top  surface  of  this  film  (steps  1  through  3).  The  process  will  then  be  repeated 
on  the  other  side,  thus  ending  up  with  a  thin  (2  microns)  layer  of  diamond  with  the  PBG  cavities  (steps 
4  through  6),  backed  by  a  silicon  substrate  on  the  top. 
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B.2.  Arbitrary  Pattern  Lithography  Using  Bose-Condensate  Based  Atomic 
Interferometry 


Summary: 

As  a  possible  alternative  to  the  E-beam  lithography,  we  have  developed  a  scheme  for  the 
lithography  of  arbitrary,  two-dimensional  nanostructures  via  matter-wave  interference. 
The  required  quantum  control  is  provided  by  a  rc/2-  n-  n/2  atom  interferometer  with  an 
integrated  atom  lens  system.  The  lens  system  is  developed  such  that  it  allows 
simultaneous  control  over  the  atomic  wave-packet  spatial  extent,  trajectory,  and  phase 

87 

signature.  We  demonstrate  arbitrary  pattern  formations  with  two-dimensional  Rb  wave 
packets  through  numerical  simulations  of  the  scheme  in  a  practical  parameter  space. 
Prospects  for  experimental  realizations  of  the  lithography  scheme  are  also  discussed. 
This  technique  may  be  able  to  produce  patterns  as  small  as  2  nm,  and  could  be  of 
widespread  utility  for  many  applications.  The  high  resolution  of  this  technique  could  also 
make  it  possible  to  make  the  PBG  patterns  with  higher  precision,  which  would  ultimately 
be  required  for  precision  quantum  computing 
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We  propose  a  scheme  for  the  lithography  of  arbitrary,  two-dimensional  nanostructures  via  matter-wave 
interference.  The  required  quantum  control  is  provided  by  a  tt/2-tt-tt/2  atom  interferometer  with  an  integrated 
atom  lens  system.  The  lens  system  is  developed  such  that  it  allows  simultaneous  control  over  the  atomic 
wave-packet  spatial  extent,  trajectory,  and  phase  signature.  We  demonstrate  arbitrary  pattern  formations  with 
two-dimensional  °'Rb  wave  packets  through  numerical  simulations  of  the  scheme  in  a  practical  parameter 
space.  Prospects  for  experimental  realizations  of  the  lithography  scheme  are  also  discussed. 
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I.  INTRODUCTION 

The  last  few  decades  have  seen  a  great  deal  of  increased 
activity  toward  the  development  of  a  broad  array  of  litho¬ 
graphic  techniques  [1,2].  This  is  because  of  their  fundamen¬ 
tal  relevance  across  all  technological  platforms.  These  tech¬ 
niques  can  be  divided  into  two  categories:  parallel 
techniques  using  light  and  serial  techniques  using  matter. 
The  optical  lithography  techniques  have  the  advantage  of 
being  fast  because  they  can  expose  the  entire  pattern  in  par¬ 
allel.  However,  these  techniques  are  beginning  to  reach  the 
limits  imposed  upon  them  by  the  laws  of  optics,  namely,  the 
diffraction  limit  [3].  The  current  state  of  the  art  in  optical 
lithography  that  is  used  in  industry  can  achieve  feature  sizes 
on  the  order  of  hundreds  of  nanometers.  Efforts  are  being 
made  to  push  these  limits  back  by  using  shorter-wavelength 
light  such  as  x  rays  [2],  but  this  presents  problems  of  its  own. 
The  serial  lithography  techniques,  such  as  electron  beam  li¬ 
thography  [1],  can  readily  attain  a  resolution  on  the  order  of 
tens  of  nanometers.  However,  because  of  their  serial  nature 
these  methods  are  very  slow  and  do  not  provide  a  feasible 
platform  for  the  industrial  mass  fabrication  of  nanodevices. 

A  different  avenue  for  lithography  presents  itself  out  of 
recent  developments  in  the  fields  of  atomic  physics  and  atom 
optics,  namely,  the  experimental  realization  of  a  Bose- 
Einstein  condensate  (BEC)  [4,5]  and  the  demonstration  of 
the  atom  interferometer  [6-12].  In  essence,  these  develop¬ 
ments  provide  us  with  the  tools  needed  in  order  to  harness 
the  wave  nature  of  matter.  This  is  advantageous  for  lithogra¬ 
phy  because  the  comparatively  smaller  de  Broglie  wave¬ 
length  of  atoms  readily  allows  for  a  lithographic  resolution 
on  the  nanometer  scale.  The  atom  interferometer  provides  a 
means  of  interfering  matter  waves  in  order  to  achieve  lithog¬ 
raphy  on  such  a  scale.  The  BEC,  on  the  other  hand,  provides 
a  highly  coherent  and  populous  source  with  which  to  per¬ 
form  this  lithography  in  a  parallel  fashion.  The  opportunity 
thus  presents  itself  to  combine  the  enhanced  resolution  of 
matter  interferometry  with  the  high  throughput  of  traditional 
optical  lithography. 

It  should  be  noted  that,  although  there  has  been  research 
activity  on  atom  lithography  [13-15]  for  a  number  of  years, 
most  of  the  work  has  involved  using  standing  waves  of  light 
as  optical  masks  for  the  controlled  deposition  of  atoms  on  a 
substrate.  The  primary  limitations  of  using  such  optical 


masks  are  that  the  lithographic  pattern  cannot  be  arbitrary 
and  that  the  resolution  of  the  pattern  is  limited  to  the  100  nm 
scale.  Since  our  scheme  uses  the  atom  interferometer,  how¬ 
ever,  it  allows  for  pattern  formation  by  self-interference  of  a 
matter  wave,  and  is  thus  unhampered  by  the  inherent  limita¬ 
tions  of  the  optical  mask  technique. 

In  this  paper  we  seek  to  demonstrate  theoretically  the  use 
of  the  atom  interferometer  as  a  platform  for  nanolithography 
by  proposing  a  technique  that  allows  for  the  manipulation  of 
a  single-atom  wave  packet  so  as  to  achieve  two-dimensional 
lithography  of  an  arbitrary  pattern  on  the  single-nanometer 
scale.  To  do  this  our  scheme  employs  a  lens  system  along 
one  arm  of  the  interferometer  that  performs  Fourier  imaging 
[3]  of  the  wave-packet  component  that  travels  along  that 
arm.  By  investigating  such  a  technique  for  a  single  atom 
wave  packet,  we  hope  to  establish  the  viability  of  using  a 
similar  technique  for  a  single  BEC  wave  packet,  which 
would  allow  for  truly  high-throughput  lithography. 

The  paper  is  organized  as  follows.  Section  II  presents  an 
overview  of  the  proposed  technique.  Sections  III  and  IV  pro¬ 
vide  a  theoretical  analysis  of  the  atom  interferometer  itself 
and  our  proposed  imaging  system,  respectively.  Section  V  is 
devoted  to  some  practical  considerations  of  the  setup  and  its 
parameter  space,  and  Sec.  VI  gives  the  results  of  numerical 
simulations.  Finally,  we  touch  upon  the  issue  of  replacing  the 
single-atom  wave  packet  with  the  macroscopic  wave  func¬ 
tion  of  a  BEC  in  Sec.  VII.  Appendixes  A  and  B  show  some 
of  the  steps  in  the  derivations. 

II.  PROPOSED  INTERFEROMETER 
A.  Principles  of  operation 

In  a  7r/2-7r-7r/2  atom  interferometer  (AI),  which  was  first 
theoretically  proposed  by  Borde  [6]  and  experimentally  dem¬ 
onstrated  by  Kasevich  and  Chu  [7],  an  atom  beam  is  released 
from  a  trap  and  propagates  in  free  space  until  it  encounters  a 
77/2  pulse,  which  acts  as  a  50-50  beam  splitter  [16-22].  The 
split  components  then  further  propagate  in  free  space  until 
they  encounter  a  77  pulse,  which  acts  as  a  mirror  so  that  the 
trajectories  of  the  split  beam  components  now  intersect.  The 
beams  propagate  in  free  space  again  until  they  encounter 
another  ir/2  pulse  at  their  point  of  intersection,  which  now 
acts  as  a  beam  mixer.  Because  of  this  beam  mixing,  any 
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phase  shift  <b  introduced  between  the  beams  before  they  are 
mixed  will  cause  an  interference  to  occur  such  that  the  ob¬ 
served  intensity  of  one  of  the  mixed  beams  at  a  substrate  will 
be  proportional  to  1+cos  cf>,  much  as  in  the  Mach-Zehnder 
interferometer  [23]  from  classical  optics.  For  our  scheme  we 
propose  the  same  type  of  interferometer,  but  with  a  single 
atom  released  from  the  trap  instead  of  a  whole  beam. 

Now,  if  we  introduce  an  arbitrary,  spatially  varying  phase 
shift  4>{x,y)  between  the  two  arms  of  the  interferometer  be¬ 
fore  they  mix,  the  intensity  of  their  interference  pattern  as 
observed  on  a  substrate  will  be  proportional  to  1 
4-  cos  4>{x,y).  Thus,  in  our  system,  we  use  an  appropriate 
choice  of  4>(x,y)  in  order  to  form  an  arbitrary,  two- 
dimensional  pattern.  This  quantum  phase  engineering  (al¬ 
ready  demonstrated  for  BECs  [24,25])  is  achieved  by  using 
the  ac-Stark  effect  so  that  I(x,y)  °c  <p(x,y),  where  I{x,y)  is 
the  intensity  of  an  incident  light  pulse. 

Also,  in  order  to  achieve  interference  patterns  on  the 
nanoscale,  4>{x,y)  must  itself  be  at  nanometer  resolution. 
However,  reliable  intensity  modulation  of  a  light  pulse  is 
limited  to  the  submicrometer  range  due  to  diffraction  effects. 
One  way  to  address  this  is  by  focusing  the  wave  packet  after 
it  is  exposed  to  the  submicrometer  resolution  phase  shift 
4>{x,y),  thereby  further  scaling  down  4>{x,y)  to  nanometer 
resolution  after  it  is  applied  to  the  wave  packet.  Our  scheme 
achieves  this  scaling  via  an  atom  lens  system. 

Additionally,  just  as  with  a  Gaussian  laser  beam,  exposing 
a  single  Gaussian  wave  packet  to  a  spatially  varying  phase 
shift  4>{x,y)  will  cause  it  to  scatter.  In  order  for  both  the 
phase-shifted  and  non-phase-shifted  components  of  the  wave 
packet  to  properly  interfere,  our  lens  system  is  also  used  to 
perform  Fourier  imaging  [30]  such  that,  at  the  substrate,  the 
phase-shifted  component  of  the  wave  packet  is  an  unscat¬ 
tered  Gaussian  that  is  properly  aligned  with  its  non-phase- 
shifted  counterpart  and  has  the  phase  information  4>{x,y)  still 
intact.  Indeed,  the  lens  system,  which  is  created  using  the 
ac-Stark  effect,  serves  the  double  purpose  of  scaling  down 
the  phase  information  <fi(x,y)  from  submicrometer  resolution 
to  single-nanometer  resolution  and  neutralizing  the  wave- 
packet  scattering  caused  by  the  same  phase  shift  4>{x,y). 

B.  Schematic 

In  our  overall  scheme,  represented  by  Fig.  1,  the  atoms 
are  treated  as  A  systems  [26-33]  (inset  B)  and  are  prepared 
in  the  ground  state  |l).  A  single-atom  trap  [34-36]  is  used  to 
release  just  one  atomic  wave  packet  along  the  z  axis.  After 
traveling  a  short  distance,  the  wave  packet  is  split  by  a  tt/2 
pulse  into  internal  states  |l)  and  |3).  The  state-|3)  component 
gains  additional  momentum  along  the  y  axis  and  separates 
from  the  state-|l)  component  after  they  both  travel  further 
along  the  z  axis.  Next,  a  tt  pulse  causes  the  two  components 
to  transition  their  internal  states  and  thereby  reflect  their  tra¬ 
jectories.  The  component  along  the  top  arm  is  now  in  the 
original  ground  state  |l)  and  proceeds  to  be  exposed  to  the 
lens  system.  The  lenses  of  the  lens  system  are  pulses  of  light 
that  intercept  the  state-|l)  component  of  the  wave  packet  at 
different  times.  By  modulating  their  spatial  intensity  in  the 
x-y  plane,  these  pulses  of  light  are  tailored  to  impart  a  par¬ 


ticular  phase  pattern  in  the  x-y  plane  to  the  wave-packet 
component  that  they  interact  with  via  the  ac-Stark  effect.  As 
shown  in  inset  B,  Fig.  1,  the  detuning  of  the  light  that  the 
lenses  are  composed  of  is  several  times  larger  for  state  |3) 
than  for  state  |l).  The  lenses  can  therefore  be  considered  to 
have  a  negligible  ac-Stark  effect  on  the  state-]  3)  wave-packet 
component  as  compared  to  the  state-|l)  component.  This  is 
important,  because  in  a  practical  situation  the  separation  be¬ 
tween  the  wave  packets  for  |l)  and  |3)  may  be  small  enough 
so  that  the  transverse  extent  of  the  lens  pulses  could  overlap 
both  wave  packets. 

The  first  light  pulse  is  intensity  modulated  to  carry  the 
phase  information  of  the  first  lens  of  the  lens  system.  It  then 
intercepts  the  state-|l)  wave-packet  component  and  adds  the 
phase  <pi(x,y).  After  some  time  the  state-|l)  component  has 
evolved  due  to  the  first  lens  such  that  it  is  an  appropriate  size 
for  exposure  to  the  phase  information  corresponding  to  the 
arbitrary  pattern  image  (inset  A).  Another  light  pulse  is  in¬ 
tensity  modulated  to  carry  the  phase  information  of  both  the 
second  lens  and  the  inverse  cosine  of  the  arbitrary  pattern. 
The  pulse  intercepts  the  state-]  1)  component  and  adds  the 
additional  phase  <f>2{. x,y).  After  some  time  a  third  light  pulse 
is  prepared  and  applied  to  the  state-]  1)  component  to  add  a 
phase  of  <^3(x,y),  which  act  as  the  third  lens  of  the  lens 
system.  Soon  after,  the  final  7 t/ 2  pulse  mixes  the  trajectories 
of  the  wave-packet  components.  A  chemically  treated  wafer 
is  set  to  intercept  the  state-|l)  component  in  the  x-y  plane. 
Due  to  the  mixing  caused  by  the  last  tt/2  pulse,  only  a  part 
of  what  is  now  the  state-]  1)  component  has  gone  through  the 
lens  system.  Because  of  the  lens  system,  it  arrives  at  the 
wafer  with  a  phase  that  is  a  scaled-down  version  of  the  im¬ 
age  phase  </>p(x,y)  =  arccos  P{x,y).  The  other  part  of  what  is 
now  the  state-]  1)  component  did  not  go  through  the  lens 
system.  There  is  therefore  a  phase  difference  of  4>P{x,y )  be¬ 
tween  the  two  parts  of  the  state-]  1)  component  and  the  wave 
packet  strikes  the  wafer  in  an  interference  pattern  propor¬ 
tional  to  1  +cos[arccos  P(x,y)]  =  1  +  P(x,y).  The  impact  with 
the  wafer  alters  the  chemically  treated  surface,  and  the  pat¬ 
tern  is  developed  through  chemical  etching. 

As  a  note,  one  preparation  for  the  wafer  is  to  coat  it  with 
a  self-assembled  monolayer  [37].  However,  Hill  et  al.  [38] 
demonstrate  an  alternate  approach  using  hydrogen  passiva¬ 
tion,  which  may  be  better  suited  for  lithography  at  the  single¬ 
nanometer  scale  due  to  its  inherent  atomic-scale  granularity. 

Finally,  note  that  the  coated  wafer  may  reflect  as  well  as 
scatter  the  pulses  of  the  lens  system.  The  phase  fronts  of  the 
wave  packets  may  potentially  be  distorted  if  exposed  to  these 
reflections  and  scatterings.  However,  this  problem  can  be 
overcome  easily  as  follows.  During  the  time  window  over 
which  the  lens  pulses  are  applied,  a  small  mirror  is  placed  at 
an  angle  in  front  of  the  wafer,  so  as  to  deflect  the  lens  pulses 
in  a  harmless  direction.  This  will  also  have  the  added  benefit 
of  not  exposing  the  wafer  to  the  lens  pulses  at  all.  Right  after 
the  last  lens  pulse  has  been  applied  and  deflected,  the  mirror 
will  be  moved  out  of  the  way,  thus  allowing  the  atomic 
waves  to  hit  the  wafer  surface. 

III.  ANALYSIS  OF  THE  INTERFEROMETER  (ttII-tt-ttII) 
A.  Formalism 

As  explained  in  the  previous  section,  we  consider  the  be¬ 
havior  of  a  single-atomic  wave  packet  in  our  formulation  of 
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FIG.  1.  A  single-atomic  wave  packet  is  released  from  the  atom  trap.  (2)  The  wave  packet  is  split  using  a  ir/2  pulse.  (3)  The  split 
components  are  reflected  by  a  tt  pulse.  (4a)  The  spatial  light  modulator  (SLM)  modulates  a  light  pulse  such  that  it  will  act  as  the  first  lens 
of  the  atom  lens  system.  (4b)  The  light  pulse  intercepts  the  wave-packet  component  that  is  in  state  |l)  and  imparts  a  phase  signature  0](x,y) 
via  the  ac-Stark  effect.  (5a)  Now  the  SLM  modulates  a  second  light  pulse  such  that  it  will  impart  both  the  phase  information  corresponding 
to  the  arbitrary  image  {arccos[P(v,y)]}  and  the  phase  information  of  the  second  lens  of  the  lens  system.  (5b)  The  second  light  pulse 
intercepts  the  same  wave-packet  component  as  the  first  one  and  imparts  the  phase  signature  (f>2{x,y).  (6a)  The  SLM  modulates  a  third  light 
pulse,  preparing  it  to  act  as  the  third  lens  of  the  lens  system.  (6b)  The  third  light  pulse  intercepts  the  same  wave-packet  component  as  the 
other  two  pulses  and  imparts  a  phase  (f>2(x,y).  (7)  Both  wave-packet  components  are  mixed  along  the  two  trajectories  by  a  ir/2  pulse.  (8) 
A  chemically  treated  wafer  intercepts  the  state-|l)  component  so  that  an  interference  pattern  forms  on  the  wafer  proportional  to  1 
+  cos{arccos[/>(x,y)]}=  1  +P(x,y).  Inset  A:  The  image  P(x,y)  that  is  to  be  transferred  ultimately  to  the  wafer.  Inset  B:  The  internal  energy 
states  of  the  wave  packet  modeled  as  a  A  system.  The  light  pulses  used  for  the  atom  lenses  have  a  much  larger  detuning  for  ground  state  |3) 
than  they  do  for  ground  state-|l)  so  that  they  effectively  only  interact  with  the  state-|l)  component  of  the  wave  packet.  The  ir/ 2  pulses  and 
the  tt  pulse  use  light  that  is  largely  detuned  for  both  ground  states. 

the  problem.  Also,  in  order  to  understand  and  simulate  the  AI 
[6-12]  properly,  the  atom  must  be  modeled  both  internally 
and  externally.  It  is  the  internal  evolution  of  the  atom  while 
in  a  laser  field  that  allows  for  the  splitting  and  redirecting  of 
the  beam  to  occur  in  the  AI.  However,  the  internal  evolution 
is  also  dependent  on  the  external  state.  Also,  while  the  exter¬ 
nal  state  of  the  atom  accounts  for  most  of  the  interference 
effects  which  result  in  the  arbitrary  pattern  formation,  the 
internal  state  is  responsible  for  some  nuances  here  as  well. 

In  following  the  coordinate  system  as  shown  in  Fig.  1,  we 
write  the  initial  external  wave  function  as 

|%(r,f  =  0)>  =  -4=exp|^|4-j  (1) 

where  f=xi+yj. 


Internally,  the  atom  is  modeled  as  a  three-level  A  system 
[26-33]  (as  shown  in  Fig.  1,  inset  B)  and  is  assumed  to  be 
initially  in  state  |l): 


IW>  =  c1(f)|l)  +  c2(f)|2)  +  c3(0|3),  (2) 


where  we  consider  c^O^l,  c2(0)  =  0,  c3(0)  =  0.  States  |l) 
and  |3)  are  metastable  states,  while  state  \2)  is  an  excited 
state. 

As  will  become  evident  later,  in  some  cases  it  is  more 
expedient  to  express  the  atom’s  wave  function  in  k  space 
[39].  To  express  our  wave  function,  then,  in  terms  of  momen¬ 
tum,  we  first  use  Fourier  theory  to  reexpress  the  external 
wave  function  as 
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\Ve(x,y,t))=  —  J  J  \^e{px,py,t))\px)\py)dplldpy,  (3) 

where  we  let  \px)=e,l-Pxl^x  and  \py)  =  el{-pylh>y .  The  complete 
wave  function  is  simply  the  outer  product  of  the  internal  and 
external  states  [Eqs.  (2)  and  (3)]: 

|^(X,y,f))  =  —  jj  [C\(Px,Py,t)\\,Px,Py)  +  C2{px,Py,t) 

X|2  ,PX,Py)  +  C3(px,py,t)\3,px,py)]dpxdpy,  (4) 

where  Cn(px,py,t)=cn(t)\<l>e(px,py,t)).  In  position  space,  the 
outer  product  gives 

|^(r,f)>  =  C\{t)\\,^e{r,t))  +  c2{t)\2pV  (,{r,t)) 

+  c3(f)|3,'Pe(r,f)).  (5) 


B.  State  evolution  in  free  space 

The  free-space  evolution  of  a  wave  function  is  fully 
derived  in  Appendix  A.  Presented  here  are  simply  the 
results  cast  in  our  particular  formalism.  For  the  free- 
space  Hamiltonian  H=J  Jll=l[{pl+p2)l2m+h(on\\n,px,py) 
X(n,px,py\dpxdpy,  if  the  wave  function  is  known  at  time  t 
=  0,  then  after  a  duration  of  time  T  in  free  space,  the  wave 
function  becomes 

\xHrJ  =  T)) 

=  ^j  j  [Ci (px,py, Q)e-^l+P^l2mh+^T\ l,p„py) 

+  C2(px,py ,  0  )e<tf+r>2mh+^T\2,px,py) 

+  C3(j>x,py,0)e-,[^2M)'2mh+^T\3,px,py)]dpxdpy, 

(6a) 


or 


\'V(r,t  =  T))  =  e-i“izc1(0)|l,''I,e(r,7))  +  e~,01l-rc2(Q) 

X|2  +  e_!"37c3(0)|3,'vIfe(r, ’,  T)) .  (6b) 


C.  State  evolution  in  tt  and  7t/2  pulse  laser  fields 

The  electromagnetic  fields  encountered  by  the  atom  at 
points  2,  3,  and  7  in  Fig.  1  that  act  as  the  tt/2,  tt,  and  77/ 2 
pulses  are  each  formed  by  two  lasers  that  are  counterpropa- 
gating  in  the  y-z,  plane  parallel  to  the  y  axis.  We  use  the 
electric  dipole  approximation  to  write  the  Hamiltonian  in 
these  fields  as 


H : 


■//£ 

J  J  n=i 


p 2  +  pi 

\  ■'  +  )\n,px,Py){n,px,py\dpxdpy 


[g«(“A*-*Ay+0A) 


+  e-i 


r 


(7) 


where  Eao  and  Eb0  are  vectors  denoting  the  magnitude  and 
polarization  of  the  fields  traveling  in  the  +  and  —y  directions, 
respectively,  e  is  the  position  vector  of  the  electron,  and  e0  is 
the  electron  charge.  Refer  to  Appendix  B  for  the  complete 
derivation  of  the  wave  function  evolution  in  these  fields. 
Only  the  results  are  presented  here. 

If  the  atom  begins  completely  in  state  1 1  ,xVe(r,t))  then 
after  a  time  T  of  evolving  in  the  above  described  fields,  the 
result  is 

|^'(r,t=r))  =  cos|y7’)|l,3Pe(r,0)) 

_  iei(a,B-coA)T+i(<t>B-<f>A)  sin|^rj|3,'4ff,(r,0))e"'(^+^, 

(8) 

where  we  have  used  the  definitions  given  in  Sec.  Ill  A.  We 
see  that  for  a  7r  pulse  (T=tt/CI),  Eq.  (8)  becomes 

\V(x,y,t  =  tt/O))  =  - 

X 1 3 ,  ^(x,  y,  0))e~^kA+k^y ,  (9) 

while  for  a  tt/2  pulse  [T=  77/(20)],  Eq.  (8)  yields 

|  ty{x,y,t=  77/(20))) 

l,3Pf(x,  V,0))  -  je1'("B-'»A)TO'2n+i(<fe-0A) 

X^=\3,Ve(x,y,0))e-KkA+kB».  (10) 

Similarly,  if  the  atom  begins  completely  in  state  |3,'vPc(r, f)), 
the  wave  function  after  a  time  T  becomes 

\^{x,y,t  =  T))  =  -  ieK°>A-°>B)T+K<l>A-<l>B)  sin^T 

X\\,We(x,yMei(kA+kB)y  +  CO 

X|3,<Pe(x,y,0)),  (11) 

so  that  for  a  77  pulse,  Eq.  (11)  gives 

| '¥(x,y,t=  tt/O))  =  -  ieiteA-naWn+HfA-tB) 

X|l,3Pe(x,y,0))e'(^+^)y,  (12) 

and  for  a  tt/2  pulse,  Eq.  (11)  becomes 

\V{x,y,t  =  77/(20)))  =  - 

x|l,^'e(x,y,0))e‘(^+^)3’ 

+  i|3,3Pf(x,y,0)).  (13) 

D.  State  evolution  through  the  whole  interferometer 

To  see  the  effects  of  phase  explicitly,  we  make  use  of  the 
analysis  that  we  have  done  for  the  state  evolution  of  the 


043606-4 


TWO-DIMENSIONAL  NANOLITHOGRAPHY  USING  ATOM... 


PHYSICAL  REVIEW  A  71.  043606  (2005) 


wave  packet.  Take  our  initial  wave  packet  I'T')  to  have  initial 
conditions  as  discussed  in  Sec.  Ill  A.  At  time  1=0  the  first 
7t/2  pulse  equally  splits  ['T)  into  two  components  |'Ta)  and 
I'Pj,)  such  that 


|^fl)  =  -  iei(^-MA)^2n+i(0fil-0A1)_L|3iqf^X;3;jO^e-,aA+/:i;h; 

V2 

(14a) 


\Vb)=-^\l,'9e{x,y,0)), 


(14b) 


=  _  —  (g'(MB-wA)7r/2ft+/(0B2_<Afl3_0A2+<fcl3) 


4.  e'(“A-"fi)'n'/2R+!(0A2-0Al+0Bl-0B2)) 

X|l,^e(x,y-y0,27’0+  r1))e“i(“>1+“3)7'0,  (17a) 


=  /  — (e'(0A2-0Al_^,A3+</lBl_0B2+fe)  _  g4“B-“A)'ITATH'(0B2-0A2)) 


X  |3,^c(x,y  -y0  -  yx,2 T0  +  Tl))e-^A+kB)y-i(a,l+lo3)T0 

(17b) 


These  have  populations 


where  we  used  Eq.  (8).  After  a  time  t=T0  of  free  space  [Eq. 
(6b)]  and  then  a  77  pulse,  Eqs.  (8)  and  (11)  yield 

•Vp  ^  —  _  e'(<uA_"s)'n'/2n+/(0A2-0Ai+0Bi-0B2)_">'32b^= 

V2 

x|l,Tfc(u>’->o,7'0)),  (15a) 

=  —  ieil'a>B~o>A^'K/^+'^B2~^A2^~io>iTa—^\3,'ire(x,y,T0)) 

v2 

Xe-^A+bib.  (15b) 

The  | ''P„)  component  becomes  shifted  in  space  by  y0  due 
to  the  momentum  it  gained  in  the  +y  direction  from  the  77 
pulse.  Now  another  zone  of  free  space  for  a  time  T0  [Eq.  (7)] 
followed  by  the  final  rr/2  pulse  [using  Eqs.  (8)  and  (11)] 
forms 


(ipi|^i)=  i[l  +cos(^0)],  <^3^3)=  ^[1  -cos(<£0)], 

(18) 

where  <f>0  =  ( 77/  O)  ( wA  -  coB)  -  cf>A  x  +  <j)B  x  +  2  4>A1 -  2  (f)B2  -  (f>A3 
+  <bHX.  We  see  that  the  state  populations  are  functions  of  the 
phase  differences  of  the  laser  fields.  Since  we  can  choose 
these  phase  differences  arbitrarily,  we  can  populate  the  states 
arbitrarily.  If  we  choose  the  phases,  for  example,  such  that 
(f>()  is  some  multiple  of  277,  then  the  wave-packet  population 
will  end  up  entirely  in  internal  state  |l). 

IV.  ARBITRARY  IMAGE  FORMATION 

If,  however,  between  the  77  pulse  and  the  second  tt/2 
pulse  we  apply  a  spatially  varying  phase  shift  (f>P(r)  to  |,'Ta), 
but  keep  4>{]  as  a  multiple  of  277,  then  the  populations  in  Eqs. 
(20)  become  instead 


\|/^  _  _  e;(<DA-«B)W2n+z'(0A2-0Ai+0Bi-0B2)~,(6'l+“,3)2b 

X-lb^Uy-To^To)) 

+  ;e'(0A2_<AAl-^’A3+<Asl_0B2+^B3)_'("'l+M3)2b 

X^\3,%(x,y-y0,2T0))e-***+k>to,  (16a) 

_  _  e;(MB-«A)i7/2n+;(0B2-0B3-0A2+0A3)~!(<«i+“3)7b 

x-lb^Uy-To^To)) 

—  /£7(w5-wAW^+,(^52-^A2)_,(wl+a,3)7b 

X  i|3,^e(x,y  -  y0,2T0))e~l(k*+k^.  (16b) 

Now  the  \'irh)  component  is  spatially  aligned  with  the 
component.  However,  another  split  occurs  because  both 
of  these  components  are  partially  in  internal  state  |3).  After 
some  further  time  7]  in  free  space,  state  |3)  has  drifted  fur¬ 
ther  in  the  +v  direction.  The  substrate  can  now  intercept  the 
two  internal  states  of  the  total  wave  function  in  separate  lo¬ 
cations.  We  write  the  state-|l)  and  -|3)  wave  functions  as 


(3pi|3pi)=  i{l  +cos[^P(7r)]}. 

<^3^3)=  “{1 -cos[0P(r)]}.  (19) 

Therefore,  if  we  let  0F(r)  =  arccos[P(r)],  where  P{r)  is  an 
arbitrary  pattern  normalized  to  1,  the  state  |l)  population  will 
be 

<¥il*i>=  ^[l  +  /m  (20) 

If  the  substrate  at  8  in  Fig.  1  intercepts  just  this  state,  the 
population  distribution  will  be  in  the  form  of  the  arbitrary 
image.  Over  time,  depositions  on  the  substrate  will  follow 
the  population  distribution,  and  thereby  physically  form  the 
image  on  the  substrate. 

A.  Imparting  an  arbitrary,  spatially  varying  phase  shift  for 
arbitrary  image  formation 

We  now  review  how  to  do  such  phase  imprinting  [24,25] 
to  a  single  wave  packet  using  the  ac-Stark  effect.  First,  con¬ 
sider  the  Schrodinger  equation  (SE)  for  the  wave  packet  ex¬ 
pressed  in  position  space: 
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r)W(r  t))  -h2 

ih'  y  =  —  V2|^(;,f))  +  V(f)VV(r,t)).  (21) 

dt  2m 

If  we  consider  a  very  short  interaction  time  t  with  the  poten¬ 
tial  V(r),  we  find 

dV^irj  +  t)) 

ih-L~^ - -«V(r)|  V(?,t+T))  (22a) 

dt 

=>|¥(r,f  +  t)>  «  | ¥(r,*)>e“(w)vWr.  (22b) 

Thus,  we  see  that  an  arbitrary  phase  shift  <fiP(r)  is  im¬ 
parted  on  the  wave  packet  in  position  space  by  choosing 
V(r)  =  (h/  T)(f>P{r).  Although  this  would  give  the  negative  of 
the  desired  phase,  it  makes  no  difference  because  it  is  the 
cosine  of  the  phase  that  gives  the  interference  pattern. 

In  order  to  create  the  arbitrary  potential  needed  to  impart 
the  arbitrary  phase  shift,  we  use  the  ac-Stark  effect  (light 
shift).  As  illustrated  in  Fig.  1  at  4b,  5b,  and  6b,  the  atom  will 
be  in  the  internal  state  |l).  If  exposed  to  a  detuned  laser  field 
that  only  excites  the  1 1)  — >■  |2)  transition,  the  eigenstates  be¬ 
come  perturbed  such  that  their  energies  shift  in  proportion  to 
the  intensity  of  the  laser  field.  A  spatially  varying  intensity 
will  yield  a  spatially  varying  potential  energy.  Specifically,  in 
the  limit  that  g/  8-  (),  where  g  is  proportional  to  the  square 
root  of  the  laser  intensity  and  S  is  the  detuning,  it  is  found 
that  the  energy  of  the  ground  state  is  approximately 
hg2/(4S).  To  impart  the  pattern  phase,  then,  we  subject  the 
atomic  wave  packet  at  4b,  5b,  and  6b  in  Fig.  1  to  a  laser  field 
that  has  an  intensity  variation  in  the  x-y  plane  such  that 

g2(r)  =  (4  S/r)(f)p(r)  =  (4<5/r)arccos[P(r)],  (23) 

where  P(r)  is  the  arbitrary  pattern  normalized  to  1  and  t  is 
the  interaction  time. 


B.  The  need  for  a  lens  system 

The  need  for  a  lens  system  for  the  atomic  wave  packet 
arises  due  to  two  separate  considerations.  First,  there  is  a 
need  for  expanding  and  focusing  the  wave  packet  in  order  to 
shrink  down  the  phase  pattern  imparted  at  5b  in  Fig.  1.  We 
have  shown  above  how  the  phase  pattern  is  imparted  using 
an  intensity  variation  on  an  impinging  light  pulse.  However, 
due  to  the  diffraction  limit  of  light,  the  scale  limit  of  this 
variation  will  be  on  the  order  of  100  nm.  This  will  cause  the 
interference  at  8  to  occur  on  that  scale.  To  reach  a  smaller 
scale,  we  require  a  lens  system  that  allows  expansion  and 
focusing  of  the  wave  packet  to  occur  in  the  transverse  plane. 
Using  such  a  system,  we  could,  for  example,  expand  the 
wave  packet  by  two  orders  of  magnitude  prior  to  5b,  impart 
the  phase  pattern  at  5b,  and  then  focus  it  back  to  its  original 
size  by  the  time  it  reaches  8.  The  interference  would  then 
occur  on  the  scale  of  1  nm. 

The  second  consideration  that  must  be  made  is  that  an 
arbitrary  phase  shift  (f>(x,y)  introduced  at  5b,  if  it  has  any 
variation  at  all  in  the  transverse  plane,  will  cause  the  wave 
packet  traveling  along  that  arm  of  the  AI  to  alter  its  momen¬ 
tum  state.  Any  free-space  evolution  after  this  point  will  make 
the  wave  packet  distort  or  go  off  trajectrory,  causing  a  noisy 


interference  or  even  eliminating  interference  at  8  all  together. 

Our  lens  system,  then,  must  accomplish  two  objectives 
simultaneously:  (1)  allow  for  an  expansion  and  focusing  of 
the  wave  packet  to  occur  and  (2)  have  the  wave  packet  prop¬ 
erly  aligned  and  undistorted  when  it  reaches  8.  To  do  this,  we 
employ  techniques  similar  to  those  developed  in  classical 
Fourier  optics  [3].  First  we  develop  a  diffraction  theory  for 
the  two-dimensional  (2D)  quantum-mechanical  wave  packet; 
then  we  use  the  theory  to  set  up  a  lens  system  that  performs 
spatial  Fourier  transforms  on  the  wave  packet  in  order  to 
achieve  the  two  above  stated  objectives. 


C.  Development  of  the  quantum-mechanical  wave-function 
diffraction  theory 

Consider  the  2D  SE  in  freespace 

<?|^(r,f))  -h2l  d2  d2\.  „ 

ih - : - =  ~+ TT  W^))-  (24) 


dt 


2m  \  dxl  dy 1 


By  inspection,  we  see  that  it  is  linear  and  shift  independent. 
If  we  can  then  find  the  impulse  response  of  this  “system”  and 
convolve  it  with  an  arbitrary  input,  we  can  get  an  exact  ana¬ 
lytical  expression  for  the  output.  To  proceed,  we  first  try  to 
find  the  transfer  function  of  the  system. 

Using  the  method  of  separation  of  variables,  it  is  readily 
shown  that  all  solutions  of  the  system  (the  2D  SE  in  free 
space)  can  be  expressed  as  linear  superpositions  of  the  fol¬ 
lowing  function: 

^(r,  t)  =  Ae^-(hl2m)  l*|2,]  (25) 

where  A  is  some  constant  and  k=kxi+kyj  can  take  on  any 
values.  Now  let  us  take  some  arbitrary  input  to  our  system  at 
time  t= 0  and  express  it  in  terms  of  its  Fourier  components: 

l'I'in(r)}  =  ~  f  Y\\Jk))e,irdk.  (26) 

2tt  J 


We  can  then  evolve  each  Fourier  component  for  a  time  T 
by  using  Eq.  (25)  to  get  the  output 

l^outW)=^  J 

=  —  f  (|<Fin(k))e-,(ft/2m)l ^e^dk 

2i t  J 

=  ^~  [  \<l>Uk)yk?dk.  (27) 

2  7 T  J 

It  follows  that 

(28) 

Our  transfer  function,  then,  for  a  free-space  system  of 
time  duration  T  is 

H(k)  =  (29) 

After  taking  the  inverse  Fourier  transform,  we  find  the 
impulse  response  to  be 

h(r)  =  -i[—  W'^kl2.  (30) 
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Finally,  convolving  this  with  some  input  to  the  system  at 
time  t= 0,  |,'Ifin(r)),  gives  the  output  at  time  t=T,  |’'P,out(r))>  to 
be 

!*-(*>  -  -  i(^T^)ei{m,2hm2  J  l^in(^)> 

Xei(ml2hT)\?'\2e-i(mlhT)f  r' dp  _  (31) 

This  expression  is  analogous  to  the  Fresnel  diffraction 
integral  from  classical  optics. 


D.  Fourier  transform  lens  scheme 

Consider  now  the  following. 

(1)  Take  as  input  some  wave  function  |^(r)),  and  use  the 
light  shift  to  apply  a  “lens”  (in  much  the  same  way  as  we 
show  above  how  to  apply  the  arbitrary  pattern  phase)  such 
that  it  becomes 


(2)  Pass  it  through  the  free-space  system  for  a  time  T 
using  the  above  derived  integral  to  get 


(  —  —  \j(m/2hT)\Pr 

\2i tUT) 


(3)  Now  use  the  light  shift  again  to  create  another  “lens” 
where  the  phase  shift  is  e-'i(mllhr)\A  -ir/2l  so  that  we  are  left 
with 


—  — )  [  \V{?))e-i{mim”' dp . 

2irhT j  J  1 


We  see  that  this  is  simply  a  scaled  version  of  the  Fourier 
transform  (FT)  of  the  input.  This  lens  system,  then,  is  such 
that 


"-«>■ -libs 


where  |<l>in)  is  the  FT  of  | 


(32) 


E.  Using  the  FT  lens  scheme  to  create  a  distortion-free 
expansion  and  focusing  system  for  applying  the  pattern  phase 

In  order  to  achieve  our  desired  goals  of  doing  expansion 
and  focusing  and  preventing  distortion,  we  propose  the  sys¬ 
tem  illustrated  in  Fig.  2(a).  We  first  input  our  Gaussian  wave 
packet  into  a  FT  scheme  with  a  characteristic  time  parameter 
T=Ta.  We  will  then  get  the  Fourier  transform  of  the  input 
(also  a  Gaussian)  scaled  by  m/{hTA).  Then,  we  give  the 
wave  packet  a  phase  shift  that  corresponds  to  the  desired 
interference  pattern  (pattern  phase)  and  put  it  through  an¬ 
other  FT  scheme  with  the  same  time  parameter  TA.  The  wave 
function  is  now  the  convolution  of  the  original  input  with  the 
pattern  phase.  Finally,  a  third  FT  scheme  is  used  with  T 
=  Tb  so  that  the  output  is  the  same  as  the  wave  function  just 
before  the  second  FT  scheme,  but  is  now  reflected  about  the 
origin  and  scaled  by  m/{hTB)  instead  of  m/(hTA).  The  pat¬ 
tern  phase,  therefore,  has  been  scaled  down  by  a  factor  of 
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Ta/Tb.  Since  both  TA  and  TB  can  be  chosen  arbitrarily,  we 
can,  in  principle,  scale  down  the  pattern  phase  by  orders  of 
magnitude.  If,  for  example,  the  pattern  phase  is  first  imparted 
on  a  scale  of  —100  nm,  we  can  choose  TA/TB  to  be  100  so 
that  at  the  output  of  our  lens  system,  it  is  on  a  scale  of 
—  1  nm.  By  scaling  down  the  pattern  phase,  we  can  scale 
down  the  interference  pattern  at  point  8  in  Fig.  1. 

Within  the  context  of  the  interferometer,  our  lens  system 
is  placed  at  4b,  5b,  and  7  in  Fig.  1.  Now,  since  the  system 
provides  us  with  the  desired  output  immediately  in  time  after 
the  final  lens  [lens  3b  in  Fig.  2(a)],  this  final  lens,  the  final 
7t/2  pulse,  and  the  substrate  6  all  need  to  be  adjacent.  If  they 
are  not,  the  wave  packet  will  undergo  extra  free-space  evo¬ 
lution  and  may  distort.  However,  such  a  geometry  is  difficult 
to  achieve  experimentally  so  we  propose  a  modification  to 
the  lens  system  [Fig.  2(b)].  Specifically,  we  can  move  the 
lens  3b  in  Fig.  2(a)  to  occur  immediately  before  lens  2a,  as 
long  as  we  rescale  it  to  account  for  the  different  wave-packet 
size  at  that  location.  We  call  the  rescaled  version  3b',  which 
is  the  same  as  3b  except  for  the  parameter  TA  in  place  of  T B. 
We  can  then  place  the  substrate  at  8  in  Fig.  2(a)  to  be  where 
the  lens  3b  previously  was;  that  is,  a  time  TB  away  from  lens 
3a.  The  final  tt/2  pulse  can  occur  anywhere  between  lens  3a 
and  the  substrate,  as  long  as  it  is  far  enough  away  from  the 
substrate  to  allow  sufficient  time  for  the  state-|3)  component 
to  separate  from  the  state-|l)  component.  To  avoid  disturbing 
the  requisite  symmetry  of  the  AI,  we  accomplish  this  by 
choosing  TB  to  be  sufficiently  large  while  leaving  the  final 
7t/2  pulse  itself  in  its  original  location.  This  geometry  will 
allow  the  substrate  to  intercept  the  state-|l)  component  ex¬ 
clusively  and  at  precisely  the  right  moment  such  that  it  does 
not  undergo  too  little  or  too  much  free-space  evolution  with¬ 
out  having  any  of  the  final  7r/2  pulse,  final  lens,  or  substrate 
adjacent.  Finally,  we  can  simplify  the  lens  system’s  imple¬ 
mentation  if  we  combine  the  lenses  that  are  adjacent.  Lenses 
lb,  2a,  3b',  and  4>P(r)  can  be  combined  into  lens  a\  lenses 
2b  and  3a  can  be  combined  into  lens  f3.  Explicitly,  lens  a  has 
phase  shift 


0«W=-(j^)H2+'n‘-<Mr)  (33) 

and  lens  /3  has  phase  shift 


<t>p(r)  =  ~ 


m  m  77 

— —  +  — —  rP+  — . 
2  fiTA  2  hTBF  2 


(34) 


Figure  2(c)  shows  the  implementation  of  the  lens  system 
within  the  context  of  the  whole  AI. 

A  cause  for  concern  may  arise  in  the  fact  that  with  the 
lens  system  in  place,  the  part  of  the  wave  packet  that  travels 
along  the  arm  without  the  lens  will  be  interfering  not  with  a 
phase-modified  version  of  itself,  but  with  a  phase-modified 
Fourier  transform  of  itself.  That  is,  the  output  of  the  lens 
system  is  a  phase-modified  Fourier  transform  of  its  input.  As 
such,  the  effective  width  of  the  wave  packet  coming  from  the 
lens  system  may  be  significantly  larger  than  the  effective 
width  of  that  coming  from  the  arm  without  lenses,  thus  caus¬ 
ing  a  truncation  of  the  pattern  formation  around  the  edges. 
This  problem  is  addressed  by  selecting  T B  such  that  the  wave 
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FIG.  2.  (A)  The  lens  system.  Each  lens  is  actually  a  pulse  of  light  with  a  transverse  intensity  modulation.  Between  lenses  la  and  lb  and 
2a  and  2b  are  free-space  regions  of  time  duration  TA,  while  between  lenses  3a  and  3b  there  is  a  free-space  region  of  duration  TB.  Lenses  la 
and  2a  give  the  wave  function  a  phase  </>ia=02a  =  _(m/2/i7,/j)|r|2,  lenses  lb  and  2b  impart  a  phase  (f>ib=(f>2b  =  -(ni/2hTA)\r\2+Tr/2,  lens  3a 
gives  a  phase  (/)3.d=-(m/2fiTB)\f\2,  and  lens  3b  gives  a  phase  (f>T,b=-(m/2hTB)\r\2  +  TT/2.  (B)  The  lens  system  from  (A)  rearranged.  The  input 
and  output  are  still  the  same,  but  the  output  is  no  longer  immediately  preceded  by  a  lens.  Lens  3b  is  the  same  as  lens  3b  from  (A)  except 
for  Ta  in  place  of  TB  so  that  it  gives  a  phase  shift  of  ^>3£/  =  -(m/2ftr4)|f|2+'n-/2.  (C)  The  modified  lens  system  in  context.  Lenses  a  and  /3 
are  composites  of  the  lenses  from  the  system  of  (B).  Between  lenses  la  and  a:  is  a  free-space  region  of  time  length  TA ,  as  well  as  between 
lenses  a  and  /?.  Between  lens  (3  and  the  substrate  is  a  free-space  region  of  time  duration  TB.  Values  of  4>a(r)  and  4>p(r)  are  as  in  Eqs.  (33) 
and  (34),  respectively. 
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packet  from  the  lens  system  is  scaled  to  have  an  effective 
width  equivalent  to  or  smaller  than  the  wave  packet  from  the 
other  arm.  Also,  because  of  the  Fourier  transform,  the  wave 
packet  coming  from  the  lens  system,  even  without  an  added 
pattern  phase,  may  have  a  different  phase  signature  from  the 
wave  packet  coming  from  the  other  arm.  Regarding  this  is¬ 
sue,  our  numerical  experiments  show  that  after  free-space 
propagation  for  a  time  on  the  order  of  the  time  scale  deter¬ 
mined  as  practical  (see  Sec.  V),  the  phase  difference  between 
the  original  wave  packet  and  its  Fourier  transform  is  very 
small  over  the  span  of  the  effective  width  of  the  wave  packet. 
Thus,  the  effect  of  this  phase  noise  on  the  interference  pat¬ 
tern  is  negligible. 

V.  SOME  PRACTICAL  CONSIDERATIONS 
A.  Wave-packet  behavior 

The  behavior  of  the  wave  packet  primarily  has  implica¬ 
tions  for  the  time  and  wave-packet  effective  width  param¬ 
eters  of  the  lithography  scheme.  As  mentioned  earlier,  the 
scale  limit  of  the  intensity  variation  that  creates  the  pattern 
phase  when  it  is  first  applied  is  — 10-7  m.  The  lens  system 
then  further  reduces  the  scale  of  the  pattern  phase  by  a  factor 
of  TaITb.  To  achieve  lithography  features  on  the  scale  of 
~1  nm,  this  ratio  needs  to  be  —100.  However,  we  must  also 
take  into  consideration  the  extent  of  the  entire  intensity 
variation.  In  other  words,  referring  to  Fig.  2(c),  the  effective 
width  of  the  wave  packet  at  lens  a  must  be  large  enough  to 
accommodate  the  entire  pattern  on  the  light  pulse  bearing  the 
phase  pattern  information.  We  assume  that  this  dimension 
will  be  on  the  order  of  a  millimeter.  We  know  that  the  wave 
packet  at  lens  a  is  a  scaled  Fourier  transform  of  the  wave 
packet  immediately  before  lens  la,  so  that  its  effective  width 
at  lens  a  is  hTA/mcrin.  This  must  be  on  the  order  of  10-3  m. 
Also,  another  way  in  which  the  time  parameters  are  restricted 
is  by  the  total  amount  of  time  that  the  atom  spends  in  the  AI. 

Now,  as  shown  earlier,  it  is  the  state-]  1)  component  in  our 
scheme  that  will  form  the  desired  interference  pattern.  The 
substrate  must  therefore  intercept  this  component  exclusive 
of  the  state-|3)  component.  Fortunately,  the  state-|3)  compo¬ 
nent  will  have  an  additional  velocity  in  the  y  direction  due  to 
photon  recoil  so  that  the  two  states  will  separate  if  given 
enough  time.  Also  recall  that  each  wave-packet  state  after  the 
final  tt/2  pulse  is  composed  of  two  elements,  one  that  went 
through  the  lens  system  and  one  that  did  not,  such  that  the 
elements  that  traveled  along  the  arm  without  the  lens  system 
will  have  larger  effective  widths  (since  the  output  of  the  lens 
system  is  smaller  than  its  input).  The  two  states  will  be  suf¬ 
ficiently  separated,  then,  when  the  state-|3)  component  has 
traveled  far  enough  in  the  +y  direction  after  the  final  tt/2 
pulse  such  that  there  is  no  overlap  of  the  larger  effective 
widths.  Since  we  know  that  photon  recoil  gives  the  state-|3) 
component  an  additional  momentum  of  2hk  in  the  +y  direc¬ 
tion,  we  have  mv=2hk.  Also,  it  can  be  shown  that  the  effec¬ 
tive  width  of  a  wave  packet  after  passing  through  free  space 
for  a  time  T  is  <r\  1  +  (77  t),  where  T=mcr/h  and  cr  is  the 
original  effective  width.  Therefore,  for  sufficient  spatial 
separation  of  the  states  (assuming  that  the  time  between  the 
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FIG.  3.  The  poposed  transition  scheme.  We  consider  5  Sl/2(F 
=  1  ,m=  1)  =  1 1),  52P1/2(F'  =  2 ,  m = 0)  =  |2),  and  52Sm(F=2, 
m=-l)  =  |l),  where  the  quantization  axis  corresponds  to  the  y  di¬ 
rection  in  Fig.  1.  cr_-polarized  light  excites  the  |l)<-^|2)  transition 
and  cr+-polarized  light  excites  the  |2)<->|3)  transition.  Both  lasers 
are  detuned  by  680  MHz.  For  the  tt/2,  t t,  and  tt/2  pulses,  the 
above  mentioned  two  transitions  are  simultaneously  excited.  For  the 
light-shift-based  lens  system,  as  shown  in  the  picture,  a  different 
transition  from  the  state  |l)  to  the  sublevel  52P1/2(F'  =2,m=  1)  is 
used,  and  Tr-polarized  light  is  applied.  Since  this  field  is  highly 
detuned  from  the  F=2<->F'=2  transitions,  the  light  shift  for  |3)  can 
be  neglected. 

final  tt/2  pulse  and  the  substrate  is  on  the  order  of  TB )  we 
need  v  X  triny/l  +{Tb/t). 

To  summarize,  our  restrictions  are 

Ta  <  10  s,  Tb  <  lO-2^, 

and 


After  using  some  simple  algebra,  we  find  that  the  first 
three  restrictions  are  satisfied  if  we  apply  the  following 

trin  10-5  m,  —  a-JT A  S  10-6  m/s  TB  ;£  lO-2^. 

We  can,  for  example,  choose  crin=10“5  m,  TA  — 10  s,  TB 
— 10”1  s.  A  simple  check  shows  that  these  choices  also  sat¬ 
isfy  the  fourth  restriction. 

Finally,  since  our  proposed  lithography  scheme  involves 
the  use  of  a  single  atom  at  a  time,  it  entails  the  drawback  of 
being  very  slow.  To  make  this  type  of  lithography  truly  prac¬ 
tical,  a  Bose-Einstein  condensate  [4,5]  would  have  to  be  used 
instead  of  a  single-atomic  wave  packet. 

B.  Proposed  87Ru  levels  and  transitions  for  the  nanolithography 

scheme 

For  practical  implementation  of  our  three-level  atom,  we 
use  the  D 1  transitions  in  :;7Rb  [40].  Figure  3  illustrates.  One 
of  the  restrictions  is  that,  in  order  to  be  able  to  neglect  spon¬ 
taneous  emission,  we  need  for  each  single  transition 
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xrx  T<§1,  (35) 

where  g0  is  the  Rabi  frequency,  S  is  the  detuning,  T  is  the 
decay  rate,  and  t  is  the  interaction  time.  Both  the  Raman 
pulse  scheme  and  the  light  shift  scheme  also  require 

go  <  8.  (36) 

We  have  the  following  relation: 

*u=(Hr2'  <3?) 

'  •'sat  ' 

If  we  assume  7sat=3  mW/cm2,  7max=2  mW/mm2,  and  T 
=  3.33  X  107  s_1,  we  find  that  go,max~8-6X  109  Hz. 

We  choose  the  relevant  Raman  A  transition  levels  to  be 
1 1)  =  5  2Sl/2(F=  1  ,ra  =  1),  |2)  =  5  2Pl/2(F=2,m  =  0),  and 
5  28ii2(F=2,m  =  -l)  =  |3),  with  the  quantization  axis  being 
in  the  y  direction  of  Fig.  1.  The  energy  difference  between 
the  levels  |l)  and  |2)  is  6.8  GHz.  The  1 1 )  — >  |2)  and  |2) 
— >|3)  transitions  are  performed  by  simultaneously  applying 
cr~ -  and  cr+-polarized  fields.  The  two  ground  states  1 1)  and  |3) 
have  equal  and  opposite  g  factors,  so  that  they  will  experi¬ 
ence  the  same  force  for  a  given  magnetic  field  gradient  used 
for  slowing  them.  For  the  ac-Stark  shift,  we  apply 
7r-polarized  light  that  couples  1 1 )  to  the  F'  =  2, m=\  excited 
state.  Because  this  field  is  highly  detuned  from  the  F 
=  2<->7?'  =  2  transition,  the  corresponding  light  shift  of  level 
|3)  can  be  neglected. 

In  order  to  satisfy  the  constraint  that  the  Rabi  frequency 
be  much  less  than  the  detuning,  we  choose  g0 =68  MHz.  This 
is  well  below  the  maximum  limit  calculated  above. 

As  far  as  the  interaction  time  for  the  tt/2  and  it  pulse 
scheme,  it  is  the  Raman  Rabi  frequency  that  is  of  interest: 

2 

n  =  — .  (38) 

28 

Using  this  in  Eq.  (36),  we  get 

a 

2-XTXHl 

8 

8 

=>CIt<  (39) 

2r 

Plugging  in  the  chosen  value  for  8  and  the  typical  value  of 
33.33  MHz  for  T,  we  find  that  12 r<S  10.2.  We  can  satisfy  this 
restraint  by  choosing  12  r=  7r  for  the  7r  pulse  and  half  as  much 
for  the  tt/2  pulse,  giving  a  pulse  duration  of  r=  77-/11 
~924  ns  for  a  t t  pulse  and  t~462  ns  for  a  tt/2  pulse. 

For  the  light  shift  we  use  the  same  7r-polarized  excitation 
of  state  |l)— >5  2P  m(F' -2,m  =  1)  as  shown  in  Fig.  3.  The 
time  constraint  in  this  case  is 

2 

—  t=2tt.  (40) 

4  8 

This  gives  an  interaction  time  of  r~3.7  /xs.  Ideally,  the 
light  shift  pulse  will  only  interact  with  the  wave  packet  in 
state  |l).  This  may  actually  be  possible  if  we  choose  TA  to  be 


large  enough  such  that  the  two  states  gain  enough  of  a  trans¬ 
verse  separation.  If,  as  by  example  above,  we  choose  TA 
~  10  s,  then  the  separation  between  the  two  states  will  be  on 
the  order  of  a  centimeter  and  there  will  be  virtually  no  over¬ 
lap  between  the  two  components  of  the  wave  packet  in  the 
separate  arms.  The  light  pulse  could  then  simply  intercept 
only  state  |3).  If,  however,  the  situation  is  such  that  the  states 
are  overlapping,  then  state  |l)  will  also  experience  the  light 
shift,  but  it  will  be  about  a  factor  of  ten  less  because  of  the 
detuning  being  approximately  ten  times  larger  for  it  than  for 
the  state-|3)  transition. 


VI.  NUMERICAL  EXPERIMENTS 


The  numerical  implementation  of  our  lithography  scheme 
was  done  by  distributing  the  wave  packets  across  finite 
meshes  and  then  evolving  them  according  to  the  Schrodinger 
equation.  This  evolution  was  done  in  both  position  and  mo¬ 
mentum  space  according  to  expediency.  To  go  between  the 
two  domains,  we  used  two-dimensional  Fourier  transform 
and  inverse  Fourier  transform  algorithms. 

The  initial  wave  packet  was  taken  in  momentum  space 
and  completely  in  internal  state  |l).  Specifically,  the  wave 
packet  was  given  by  the  Fourier  transform  of  Eq.  (1): 


|<J>e(M  =  0)) 


(41) 


The  evolution  of  the  wave  packets  in  the  tt  and  tt/2 
pulses  was  done  in  momentum  space  in  order  to  be  able  to 
account  for  the  different  detunings  that  result  for  each  mo¬ 
mentum  component  due  to  the  Doppler  shift.  Specifically,  we 
numerically  solved  Eq.  (B15)  for  the  different  components  of 
the  7-space  wave  packet  mesh,  then  applied  the  inverse  of 
the  transformation  matrix  given  by  Eq.  (B9)  to  go  to  the 
original  basis. 

Outside  of  the  lens  system,  the  free-space  evolution  of  the 
wave  packets  was  also  done  in  momentum  space.  This  was 
achieved  easily  by  using  Eqs.  (A4).  Within  the  lens  system, 
however,  it  was  more  computationally  efficient  to  use  Eq. 
(31)  for  the  free-space  evolution  because  of  the  need  to  apply 
the  lenses  in  position  space.  The  results  of  using  Eq.  (32) 
were  initially  cross-checked  with  the  results  of  using  Eqs. 
(A4)  and  were  found  to  agree. 

Figure  4(a)  is  a  targeted  (arbitrary)  pattern.  Figures  4(b) 
and  4(c)  demonstrate  the  formation  of  the  arbitrary  pattern 
by  interference  of  the  state-|l)  wave  packets  at  the  output  of 
the  interferometer.  Both  figures  were  the  result  of  applying 
the  same  arbitrary  pattern  phase,  but  Fig.  4(b)  was  formed 
without  any  shrinking  implemented  (i.e.,  TA-TB).  Figure 
4(c),  however,  demonstrates  the  shrinking  ability  of  the  lens 
system  by  yielding  a  version  of  Fig.  4(b)  that  is  scaled  by  a 
factor  of  2{Ta/Tb=2).  The  length  scales  are  in  arbitrary  units 
due  to  the  use  of  naturalized  units  for  the  sake  of  computa¬ 
tional  viability. 


VII.  SUGGESTIONS  FOR  EXTENSION  TO  BEC 

As  mentioned  above,  in  order  to  make  the  lithography 
scheme  truly  practical,  a  Bose-Einstein  condensate  is  re- 
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FIG.  4.  (a)  A  targeted  arbitrary  (e.g.,  a  tic-tac-toe  board  here) 
image,  (b)  The  arbitrary  image  is  now  formed  with  the  lens  system 
in  place,  but  without  any  scaling.  We  see  that  it  is  a  more  complex 
pattern  than  just  a  simple  periodic  structure  such  as  sinusoidal 
fringes,  (c)  The  same  image  as  in  (b)  is  formed  with  the  lens  system 
still  in  place,  but  a  scaling  factor  of  2  has  been  used  to  shrink  the 
pattern. 

quired  in  place  of  the  single  atom.  Indeed,  the  self¬ 
interference  of  a  BEC  has  already  been  demonstrated 
[41,42].  The  difficulty  in  using  the  BEC  for  controlled  imag¬ 
ing,  however,  arises  from  the  nonlinear  term  in  the  Gross- 
Pitaevskii  equation  (GPE).  Our  lens  system,  for  example, 
would  not  be  valid  as  it  was  developed  from  the  linear  SE. 

One  approach  to  getting  around  this  problem  is  to  try  to 
eliminate  the  nonlinear  term  in  the  GPE.  Specifically,  the 
GPE  for  the  BEC  takes  the  form 

(AP  l  -ft2  \ 

=  [  — V'+V+GoW2]*,  (42) 

where  the  nonlinear  term  coefficient  is  U0=4Trh2a/m  and  a 
is  the  scattering  length  for  the  atom.  It  has  been  demon¬ 
strated  for  s7Rb  that  the  scattering  length  can  be  tuned  over  a 
broad  range  by  exposing  the  BEC  to  magnetic  fields  of  vary¬ 
ing  strength  near  Feshbach  resonances  [43,44].  The  relation¬ 


ship  between  the  scattering  length  and  the  applied  magnetic 
field  B  when  near  a  Feshbach  resonance  can  be  written  as 


A 


^  ^peak 


(43) 


where  ah„  is  the  background  scattering  length,  Bpeak  is  the 
resonance  position,  and  A=Bzer0-Bpeak.  Setting  B=Bzei0 
would  therefore  set  the  scattering  length  to  zero  and  elimi¬ 
nate  the  nonlinear  term  in  the  GPE.  While  the  atom-atom 
interaction  may  not  be  completely  eliminated  in  reality  due 
to  the  fluctuation  in  density  that  we  wish  to  effect  through 
the  lens  system,  it  is  worth  investigating  if  it  could  be  made 
to  be  negligible  over  an  acceptable  range.  We  could  then  use 
our  previously  developed  lens  system  to  perform  the  imaging 
and  thereby  interfere  a  large  number  of  atoms  simulta¬ 
neously.  Alternatively,  one  must  redevelop  the  design  of  the 
lenses  and  the  imaging  optics  as  applied  to  the  equation  of 
motion  for  a  BEC  [Eq.  (42)]  for  a  nonzero  value  of  U0.  This 
effort  is  in  progress. 
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APPENDIX  A:  STATE  EVOLUTION  IN  FREE  SPACE 


In  free  space,  the  Hamiltonian  can  be  expressed  in  the 
momentum  domain  as 


H 


+  hco„  I n,px,py)(n,px,py\dpjpy. 


(Al) 

where  u>n  is  the  frequency  corresponding  to  the  eigenenergy 
of  internal  state  | n).  For  a  single  momentum  component  (px 
=px o  and  Py=pyo),  the  Hamiltonian  for  the  total  evolution  in 
momentum  space  is  given  by 


H-- 


P.xO  +  PyO 

2m 


+  ha> i  0 


0 


Pi  o  +  Plo 

2m 


+  ha>2  0 


Pm  +  pIo 

2m 


+  hu>2 
(A2) 

Using  this  in  the  SE,  we  get  the  equations  of  the  ampli¬ 
tude  evolution  in  momentum  space: 


Ci(pM,Py0,t)  =  -  ~ 


'  I  P2x0  +  P2vq 
2m 

2  ,  2 


+  )C\{pxQ,py0,t), 


Cl(px0,PyoJ)  =  ~  +  h(o2)C2(Px0’Py0’t)’ 


(A3) 
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C3(Px0’Py0’t)  =  -  ^(PV0^/;V°  +  haj^Jc3{pxQ,py0,t). 
These  yield  the  solutions 
Ctip^PyoJ)  = 

C2(PxO’PyO’ t)  =  C2(px0,py0, ( A4) 


C3(px0,py0,t)  =  C^Pyo^e-^Wyo)'2"^'. 

We  see  that  if  the  wave  function  is  known  at  time  t= 0, 
then  after  a  duration  of  time  T  in  free  space,  the  wave  func¬ 
tion  becomes 

|¥(r,f  =  7)> 

=  2^11  VCx(px,py,  0)e<^l2m^T\l,px,py) 

+  C2(px,py,  0  )e^+P>mh^T\2,px,py) 

+  Q(/Jx>/7>'>0)e“'L^()/2mft+“3]7'|3,px,py)]Jp^p>,. 

(A5) 

We  can  also  write  it  as 


Eb  =  Eb  o  cos(«Br  +  kBy  +  <f>B ) 


— —[  pi(",B,+kBy+'t’B)  +  e-'(u>B,+kBy- 
2 


(B2) 


where  Eao  and  Eb0  are  vectors  denoting  the  magnitude  and 
polarization  of  their  respective  fields.  Keeping  in  mind  that 
our  wave  function  is  expressed  in  the  momentum  domain, 
we  take  position  as  an  operator. 

The  Hamiltonian  here  is  expressed  as  the  sum  of  two 
parts:  H=H0+H1.  The  first  part  corresponds  to  the  noninter¬ 
action  energy: 


fin 


■JJs 

J  J  n=  1 


3/2  2 

Px+Py 


2m 


+  fuo„  \n,px,p){n,px,py\dPxdp 


(B3) 


The  second  part  accounts  for  the  interaction  energy,  for 
which  we  use  assumption  (3)  from  above  to  make  the  elec¬ 
tric  dipole  approximation  and  get 


£ 

Hl  =  -  eQe  ■  — +  e-!'(&>A«-£Ay+0A)j 


I^(r,r  =  T))  =  e-i“lTc1(0)|l,'?re(r,7))  +  e_,'“2rc2(0) 

X|2,¥e(r,r)>  +  e~i,oiTc3{0)\?>,'Ve{r,Tj). 

(A6) 


APPENDIX  B:  STATE  EVOLUTION  IN  n  AND  rr/2  PULSE 
LASER  FIELDS 


^e'('ufl,+*si'+<As) 


(B4) 


where  e  is  the  position  vector  of  the  electron,  and  eQ  is  the 
electron  charge.  Now,  seeing  that  expressions  of  the  form 
( n\s-EA0\n )  and  ( n\e-EB0\n )  are  zero,  and  using  assumption 
(4),  we  can  express  Eq.  (B4)  as 


The  electromagnetic  fields  encountered  by  the  atom  at 
points  2,  3,  and  7  in  Fig.  1  that  act  as  the  i t  and  tt/2  pulses 
are  each  formed  by  two  lasers  that  are  counterpropagating  in 
the  y-z  plane  parallel  to  the  y  axis.  We  will  refer  to  the  laser 
propagating  in  the  +y  direction  as  EA,  and  the  one  propagat¬ 
ing  in  the  —y  direction  as  EB.  In  deriving  the  equations  of 
motion  under  this  excitation,  we  make  the  following  assump¬ 
tions:  (1)  the  laser  fields  can  be  treated  semiclassically  [45], 
(2)  the  intensity  profiles  of  the  laser  fields  forming  the  7r  and 
7t/2  pulses  remain  constant  over  the  extent  of  the  atomic 
wave  packet,  (3)  the  wavelengths  of  the  lasers  are  signifi¬ 
cantly  larger  than  the  separation  distance  between  the 
nucleus  and  electron  of  the  atom,  (4)  EA  excites  only  the 
1 1 ) •«— >  J 2)  transition  and  EB  only  the  |3)<->|2)  transition,  (5) 
Ea  and  Eb  are  far  detuned  from  the  transitions  that  they 
excite,  and  (6)  EA  and  EB  are  of  the  same  intensity. 

Using  assumptions  (1)  and  (2),  we  write  the  laser  fields  as 


^(|l>P*.PyX2>P;c,PVl  +  |2,^,Py)<l,^,Py|) 


EgB 


X[e'K'-Uv+0A>  +  +  '2M(\3,Px,Pv) 

X(2,Px,Py\  +  \2,px,py)(3,px,py\)[e^>°B'+kBy+^ 


.  e_'('uB,+&gy+0fl)j 


dPxdPy, 


(B5) 


where  we  let  gA=(l|e-£A0|2)  =  (2|8-£A0|l)  and  gB 
=  <3|e-£*)|2)=<2|  e-Ego^)-  Finally,  we  can  use  the  identities 
[39] 


0iky  . 


■2//  K 


PX’Py)(n,Px,Py  -  hk\dPxdp  ,  (B6a) 


Ea  =  Ea o  cos («Af  -  kAy  +  <f>A) 


Ea  q 
2 


| -  i(aAt-kAy+$A) 


+  e-i(<oAf-I:Ay+0A)] 


(Bl) 


-iky , 


■s/J  (»• 


px,Py)(n,px,py  +  hk\dPxdPy,  (B6b) 


and 


and  the  rotating  wave  approximation  [45]  in  Eq.  (B5)  to  give 


043606-12 


TWO-DIMENSIONAL  NANOLITHOGRAPHY  USING  ATOM... 


PHYSICAL  REVIEW  A  71.  043606  (2005) 


l  =  J  J  1  ,px,py)(2,px,py  +  hkA\ 

+  ^e-*“*"+J\2,p„py  +  hkA)(l,px,py\  +  b > 

X  |3  ,px,py  +  hkA  +  hkB)(2,px,py  +  hkA\  + 


X  \2,px,p  +  hkA){3,px,p  +  hkB  +  hkA 


dpxdpy.  (B7) 


We  note  that  the  full  interaction  between  the  internal 
states  1 1),  |2),  and  |3)  occurs  across  groups  of  three  different 
momentum  components:  \px,py),  \ px,py+hkA),  and  \px,py 
+  hkA  +  hkB).  This  can  be  understood  physically  in  terms  of 
photon  absorption  and  emission  and  conservation  of  momen¬ 
tum.  Keeping  in  mind  assumption  (4),  if  an  atom  begins  in 


state  |l  ,pxo,pYo)  and  absorbs  a  photon  from  field  EA,  it  will 
transition  to  internal  state  \2)  because  it  has  become  excited, 
but  it  will  also  gain  the  momentum  of  the  photon  ( hkA )  trav¬ 
eling  in  the  +y  direction.  It  will  therefore  end  up  in  state 
\2,px0,px()+hkA).  Now  the  atom  is  able  to  interact  with  field 
EB,  which  can  cause  stimulated  emission  of  a  photon  with 
momentum  hkB  in  the  —y  direction.  If  such  a  photon  is  emit¬ 
ted,  the  atom  itself  will  gain  an  equal  momentum  in  the 
opposite  direction,  bringing  it  into  external  state  \pXQ,pxo 
+  hkA+hkB).  The  atom  will  also  make  an  internal  transition 
to  state  |3)  because  of  the  deexcitation.  The  total  state  will 
now  be  \3  ,pxQ,pyi)+hkA  +  hkg}.  We  thereby  see  that  our 
mathematics  is  corroborated  by  physical  intuition. 

Getting  back  to  the  Hamiltonian,  we  look  at  the  general 
case  of  one  momentum  grouping  so  that  we  get  in  matrix 
form  H=H0+HX  from  Eqs.  (B3)  and  (B7): 


H-- 


2  .  2 


Px  +  P}  +  h(Ol  ^V("A'+« 
2m  2 


P.v  +  (Pv  +  **a)-  > 

2  2m  "  2 


0 


p\  +  (Py  +  hkA  +  hkg )2 


2  2m 

In  order  to  remove  the  time  dependence  we  apply  some  transformation  Q  [39]  of  the  form 


+  ha>2 


(B8) 


so  that  the  SE  becomes 


i(6»p+0l)  q 

0 

)  £*(  Hyt+tb'p 

0 

)  0 

ei(03t+(/>3) 

ih 


dl'P) 

dt 


where  |'4f)  =  2|1P)  and  H=QHQ  1  +  ih(0QI dt)Q  1 .  The  matrix  representation  is 


(B9) 


(BIO) 


H  = 


2  ,  2 
Px  +  Px 

2m 


+  ha I,  -  h0l 


^i(oiA+a]-d2)i+i(4>A+tl>]  ~<hi  q 

2 


■' (Py  '  Yl'  +  jita  —h6  — ^  g-|(Mg+^-g2)f-i(0B+^-02) 

2  2m  -  -  2 


^Bei(aB+e3-e2),+i(M-4>2)  P:<+(Py  +  *kA  +  fikB)-  _ 

2  2m  3  3 


(Bll) 


Choosing  0l  =  -a>A,  02= 0,  02  =  -a>B,  <t>l  =  -(f>A,  <f>2  =  0,  and  <f)3  =  -(f>B,  Eq.  (Bll)  becomes 
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H  = 


E\ (Px’Pv)  +  hv  l  +  ~ 


Ha 

2 

0 


E2(px,py)  +  a  co2 

Hb 

2 


0 

Hb 

2 

EyiPx’Py)  +  ?l(03  + 


(B12) 


where  we  have  taken 


E\(px,py)  ■ 


2  ,  2 

Px  +  Pv 

2m 


E2(PxiPy)  = 


P2X  +  (Py  +  hkA)2 


E2(Px,Py)  = 


2m 

p;  +  jpy  +  hkA  +  hkB )2 
2m 


(B13a) 


(B13b) 


(B13c) 


In  order  to  further  simplify  the  analysis,  we  set  the  zero 
energy  at  El(pxQ,py0)  +  hu>l  +  hcoA  for  some  specific  momen¬ 
tum  group  with  px=px o  and  py=py o-  Also,  since  a>A  and  coB 
can  be  chosen  independently,  we  can  let  E3(px 0,py0)+h(o3 
+  hcoB=El(px0,pv0)  +  hojl  +  hcoA.  With  the  energies  thus  set, 
Eq.  (B12)  becomes 


H  = 


0 

Ha 

2 

0 


Ha 

2 

-  S 
Hb 


0 

Hb 

2 

0 


(B14) 


where  S=[El(px0,py0)  +  hajl+h(oA\-[E2{px0,py0)  +  hM2].  Us¬ 
ing  this  Hamiltonian  in  Eq.  (BIO),  we  get  the  equations  of 
motion  as 

CiiPxO’PyO’1 )  =  -  i—C2(px0,Pyo+hkA,t),  (B15a) 


~  8  A _ 

C2(Pxo,pyo  +  hkA,  t )  =  -  k—Ci(pxQ,py0,t)  +  iSC2(pxo,py o 


.gB- 


+  hkA,t)  -  i—C3(px0,py0  +  hkA 


+  fi  kBl  f) , 


(B15b) 


. gB ; 


C3(pxQ,pyo  +  hkA  +  hkB,t)  =  -  i—C2(px0,Py  +  hkA,t). 

(B15c) 

Assumption  (5)  allows  us  to  make  the  adiabatic  approxi¬ 
mation  so  that  we  can  set  C2(px0,py0,t)~0,  and  assumption 


(6)  gives  us  gA=gB=go-  Equations  (B15)  then  simplify 


to 


.go 


Ci(Px0,Py0,t)  =  -  i^Ct(px0,py0,t) 


^2 

-  i  —fPiip&Pyo  +  hkA  +  hkB,  t) , 


(B16a) 


g2  ~ 

C3{px0,Pyo  +  hkA  +  hkB,t )  =  -  ifsCt(px0,py0,t) 


g2  ~ 

-  i^C3[px0,py0  +  hkA  +  hkB,  t) , 

(B16b) 

where  we  have  chosen  to  neglect  state  C2  from  here  on  due 
to  the  adiabatic  approximation.  We  can  now  use  another 
transformation  on  this  system  to  make  it  more  tractable.  Let 

c  1  (p^Pyo,  t)  =  C,  (px0,py0,  t)el(s^4S)' ,  (B 1 7  a) 

C3(px0,py0  +  hkA  +  hkB,t)  =  C3(px0,pyQ  +  hkA  +  hkB,t)e'igol4S>' . 

(B17b) 

The  system  in  Eqs.  (B16)  then  becomes 

2 

.go 


»  60  ~ 

Ci(px0,py0,t)  =  ~  i^C3(px0,py0  +  hkA  +  hkB,t), 


(B18a) 


-go. 


C3(px0,pyo  +  hkA  +  hkB,  t)  =  -i  ,  (px0,py  0,  t) . 

(B18b) 

Solving  this  and  reversing  the  transformations  of  Eqs. 
(B17)  and  (B9),  we  arrive  at 

/ft 

Ci(p*o>PyO’t)  =  Ci(px0,p3,0,0)cosl  —t 


C3(px0,Pyo  +  hkA 


+  hkBA))e-,(ol''-nn),-,rl,B  sin ( 


(B19a) 


043606-14 


TWO-DIMENSIONAL  NANOLITHOGRAPHY  USING  ATOM... 


PHYSICAL  REVIEW  A  71,  043606  (2005) 


CliPxO’PyO  +  hkA  +  hkB,t) 


X 


+  C3(px0,py0  +  hkA  +  Mb,0)cos^  —  t J ,  (B  19b) 

where  we  let  Ll=g^l2S.  It  should  be  noted,  however,  that 
these  solutions  were  arrived  at  only  for  the  specific  momen¬ 


tum  group  where  px=px -o  and  px=py o-  This  was  the  case 
where  both  laser  fields  were  equally  far  detuned.  Other  mo¬ 
mentum  groups  will  have  slightly  different  solutions  due  to 
the  Doppler  shift,  which  causes  the  detunings  to  be  per¬ 
turbed.  For  a  more  accurate  description,  we  need  to  numeri¬ 
cally  solve  each  momentum  group’s  original  three  equations 
of  motion  without  making  any  approximations.  This  is  what 
we  do  in  our  computational  model.  For  a  basic  phenomeno¬ 
logical  understanding  of  the  interferometer,  however,  it  is 
sufficient  to  assume  that  the  above  analytical  solution  is  ac¬ 
curate  for  all  momentum  components. 


[1]  Handbook  of  Micro-lithography,  Micromachining  and  Micro¬ 
fabrication,  edited  by  P.  Rai-Choudhury,  SPIE  Press,  Belling¬ 
ham,  WA  (1979). 

[2]  Microlithography  Science  and  Technology,  edited  by  I.  R. 
Sheats  and  Bruce  W.  Smith  (Marcel  Dekker,  New  York,  1997) 
Chap.  7. 

[3]  Joseph  W.  Goodman,  Introduction  To  Fourier  Optics 
(McGraw-Hill,  New  York,  1996). 

[4]  M.  H.  Anderson  et  al.,  Science  269,  198  (1995). 

[5]  K.  B.  Davis  et  al,  Phys.  Rev.  Lett.  75,  3969  (1995). 

[6]  C.  J.  Borde,  Phys.  Lett.  A  140,  10  (1989). 

[7]  M.  Kasevich  and  S.  Chu,  Phys.  Rev.  Lett.  67,  181  (1991). 

[8]  L.  Gustavson,  P.  Bouyer,  and  M.  A.  Kasevich,  Phys.  Rev.  Lett. 
78,  2046  (1997). 

[9]  M.  J.  Snadden,  J.  M.  McGuirk,  P.  Bouyer,  K.  G.  Haritos,  and 
M.  A.  Kasevich,  Phys.  Rev.  Lett.  81,  971  (1998). 

[10]  J.  M.  McGuirk,  M.  J.  Snadden,  and  M.  A.  Kasevich,  Phys. 
Rev.  Lett.  85,  4498  (2000). 

[11]  Y.  Tan,  J.  Morzinski,  A.  V.  Turukhin,  P.  S.  Bhatia,  and  M.  S. 
Shahriar,  Opt.  Commun.  206,  141  (2002). 

[12]  D.  Keith,  C.  Ekstrom,  Q.  Turchette,  and  D.  E.  Pritchard,  Phys. 
Rev.  Lett.  66,  2693  (1991). 

[13]  J.  H.  Thywissen  and  M.  Prentiss,  e-print  physics/0209084. 

[14]  B.  Brezger  et  al.,  Europhys.  Lett.  46,  148  (1999). 

[15]  R.  Stiitzle  et  al,  J.  Opt.  B:  Quantum  Semiclassical  Opt.  5, 
S164  (2003). 

[16]  D.  S.  Weiss,  B.  C.  Young,  and  S.  Chu,  Phys.  Rev.  Lett.  70, 
2706  (1993). 

[17]  T.  Pfau  et  al.,  Phys.  Rev.  Lett.  71,  3427  (1993). 

[18]  U.  Janicke  and  M.  Wilkens,  Phys.  Rev.  A  50,  3265  (1994). 

[19]  R.  Grimm,  J.  Soding  and  Yu.  B.  Ovchinnikov,  Opt.  Lett.  19, 
658  (1994). 

[20]  T.  Pfau,  C.  S.  Adams,  and  J.  Mlynek,  Europhys.  Lett.  21,  439 
(1993). 

[21]  K.  Johnson,  A.  Chu,  T.  W.  Lynn,  K.  Berggren,  M.  S.  Shahriar, 
and  M.  G.  Prentiss,  Opt.  Lett.  20,  1310  (1995). 

[22]  M.  S.  Shahriar,  Y.  Tan,  M.  Jheeta,  J.  Morzinsky,  P.  R.  Hem- 
mer,  and  P.  Pradhan,  e-print  quant-ph/0311163;  J.  Opt.  Soc. 


Am.  B  (2005)  (in  press). 

[23]  B.  E.  A.  Saleh  and  M.  C.  Teich,  Fundamentals  of  Photonics 
(Wiley-Interscience,  New  York,  1991). 

[24]  L.  Dobrek  et  al.,  Phys.  Rev.  A  60,  R3381  (1999). 

[25]  J.  Denschlag  et  al..  Science  287,  97  (2000). 

[26]  S.  Swain,  J.  Phys.  B  15,  3405  (1982). 

[27]  M.  Prentiss,  N.  Bigelow,  M.  S.  Shahriar,  and  P.  Hemmer,  Opt. 
Lett.  16,  1695  (1991). 

[28]  P.  R.  Hemmer,  M.  S.  Shahriar,  M.  Prentiss,  D.  Katz,  K.  Berg¬ 
gren,  J.  Mervis,  and  N.  Bigelow,  Phys.  Rev.  Lett.  68,  3148 
(1992). 

[29]  M.  S.  Shahriar  and  P.  Hemmer,  Phys.  Rev.  Lett.  65,  1865 
(1990). 

[30]  M.  S.  Shahriar,  P.  Hemmer,  D.  P.  Katz,  A.  Lee,  and  M.  Prent¬ 
iss,  Phys.  Rev.  A  55,  2272  (1997). 

[31]  P.  Meystre,  Atom  Optics  (Springer  Verlag,  Berlin,  2001). 

[32]  Atom  Interferometry,  edited  by  P.  Berman  (Academic  Press, 
New  York,  1997). 

[33]  Fiber  Optic  Rotation  Sensors,  edited  by  S.  Ezekiel  and  A. 
Arditty  (Springer- Verlag,  Berlin,  1982). 

[34]  J.  McKeever  et  al.,  Phys.  Rev.  Lett.  90,  133602  (2003). 

[35]  S.  J.  van  Enk  et  al.,  Phys.  Rev.  A  64,  013407  (2001). 

[36]  A.  C.  Doherty  et  al.,  Phys.  Rev.  A  63,  013401  (2001). 

[37]  K.  K.  Berggren  et  al..  Science  269,  1255  (1995). 

[38]  S.  B.  Hill  et  al.,  Appl.  Phys.  Lett.  74,  2239  (1999). 

[39]  M.  S.  Shahriar,  M.  Jheeta,  Y.  Tan,  P.  Pradhan,  and  A.  Gangat, 
Opt.  Commun.  243,  183  (2004). 

[40]  A.  Corney,  Atomic  and  Laser  Spectroscopy  (Oxford  Science, 
Oxford,  1976). 

[41]  M.  R.  Andrews  et  al..  Science  275,  637  (1997). 

[42]  Y.  Torii  et  al.,  Phys.  Rev.  A  61,  041602  (2000). 

[43]  A.  Marte,  T.  Volz,  J.  Schuster,  S.  Durr,  G.  Rempe,  E.  G.  M. 
van  Kempen,  and  B.  J.  Verhaa,  Phys.  Rev.  Lett.  89,  283202 
(2002). 

[44]  T.  Volz,  S.  Durr,  S.  Ernst,  A.  Marte,  and  G.  Rempe,  Phys.  Rev. 
A  68,  010702(R)  (2003). 

[45]  L.  Allen  and  J.  H.  Eberly,  Optical  Resonance  and  Two-Level 
Atoms  (Dover,  New  York,  1974). 


043606-15 


B.3.  Development  of  a  Quantum  Computing  Architecture  Using  Atomic  Ensembles 
as  Quantum  Bits 

Summary: 

Given  the  technical  difficulties  encountered  in  realizing  a  quantum  computer 
based  on  NV-diamond,  coupled  with  the  prospect  that  it  may  be  very  difficult  to  achieve 
a  very  small  decoherence  rate  in  such  a  system  (a  problem  endemic  to  all  solid  state  QC 
ideas),  we  worked  on  developing  alternative  ideas  for  robust  and  high-fidelity  quantum 
computing,  based  on  atomic  ensembles  as  quantum  bits. 

Current  proposals  focusing  on  neutral  atoms  for  quantum  computing  are  mostly 
based  on  using  single  atoms  as  quantum  bits  (qubits),  while  using  cavity  induced 
coupling  or  dipole-dipole  interaction  for  two-qubit  operations.  An  alternative  approach  is 
to  use  atomic  ensembles  as  quantum  bits.  However,  when  an  atomic  ensemble  is  excited, 
by  a  laser  beam  matched  to  a  two-level  transition  (or  a  Raman  transition)  for  example,  it 
leads  to  a  cascade  of  many  states  as  more  and  more  photons  are  absorbed.  In  order  to 
make  use  of  an  ensemble  as  a  qubit,  it  is  necessary  to  disrupt  this  cascade,  and  restrict  the 
excitation  to  the  absorption  (and  emission)  of  a  single  photon  only.  We  have  shown  how 
this  can  be  achieved  by  using  a  new  type  of  blockade  mechanism,  based  on  the  light-shift 
imbalance  (LSI)  in  a  Raman  transition.  We  developed  first  a  simple  example  illustrating 
the  concept  of  light-shift  imbalance  induced  blockade  (LS11B)  using  a  multi-level 
structure  in  a  single  atom,  and  show  verifications  of  the  analytic  prediction  using 
numerical  simulations.  We  then  extended  this  model  to  show  how  a  blockade  can  be 
realized  by  using  LSI  in  the  excitation  of  an  ensemble.  Specifically,  we  showed  how  the 
LSI1B  process  enables  one  to  treat  the  ensemble  as  a  two-level  atom  that  undergoes  fully 
deterministic  Rabi  oscillations  between  two  collective  quantum  states,  while  suppressing 
excitations  of  higher  order  collective  states. 

Finally,  we  developed  an  explicit  scheme  that  shows  how  such  a  system  can  be 
used  as  a  quantum  bit  (qubit)  for  quantum  communication  and  quantum  computing. 
Specifically,  we  showed  how  to  realize  a  controlled-NOT  (CNOT)  gate  using  the 
collective  qubit  and  an  easily  accessible  ring  cavity  or  a  tapered  nano-fiber,  via  an 
extension  of  the  so-called  Pellizzari  scheme.  We  also  describe  how  multiple,  small-scale 
quantum  computers  realized  using  these  qubits  can  be  linked  effectively  for 
implementing  a  quantum  internet.  We  developed  an  explicit  description  of  the  details  of 
the  energy  levels  and  transitions  in  an  Rb  atom  that  could  be  used  for  implementing 
these  schemes. 
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Abstract 

Current  proposals  focusing  on  neutral  atoms  for  quantum  computing  are  mostly  based  on  using  single  atoms  as  quantum  bits 
(qubits),  while  using  cavity  induced  coupling  or  dipole-dipole  interaction  for  two-qubit  operations.  An  alternative  approach  is  to 
use  atomic  ensembles  as  quantum  bits.  However,  when  an  atomic  ensemble  is  excited,  by  a  laser  beam  matched  to  a  two-level  transition 
(or  a  Raman  transition)  for  example,  it  leads  to  a  cascade  of  many  states  as  more  and  more  photons  are  absorbed  [R.H.  Dicke,  Phys. 
Rev.  93  (1954)  99].  In  order  to  make  use  of  an  ensemble  as  a  qubit,  it  is  necessary  to  disrupt  this  cascade,  and  restrict  the  excitation  to 
the  absorption  (and  emission)  of  a  single  photon  only.  Here,  we  show  how  this  can  be  achieved  by  using  a  new  type  of  blockade  mech¬ 
anism,  based  on  the  light-shift  imbalance  (LSI)  in  a  Raman  transition.  We  describe  first  a  simple  example  illustrating  the  concept  of 
light-shift  imbalance  induced  blockade  (LSIIB)  using  a  multi-level  structure  in  a  single  atom,  and  show  verifications  of  the  analytic 
prediction  using  numerical  simulations.  We  then  extend  this  model  to  show  how  a  blockade  can  be  realized  by  using  LSI  in  the  exci¬ 
tation  of  an  ensemble.  Specifically,  we  show  how  the  LSIIB  process  enables  one  to  treat  the  ensemble  as  a  two-level  atom  that  under¬ 
goes  fully  deterministic  Rabi  oscillations  between  two  collective  quantum  states,  while  suppressing  excitations  of  higher  order  collective 
states. 

©  2007  Elsevier  B.V.  All  rights  reserved. 

OCIS:  020.1670;  270.0270;  270.1670 


Many  different  technologies  are  currently  being  pursued 
for  realizing  a  quantum  computer  (QC).  One  of  the  most 
promising  approaches  involve  the  use  of  neutral  atoms. 
This  approach  is  particularly  attractive  because,  in  princi¬ 
ple,  it  is  possible  to  achieve  very  long  decoherence  times 
and  very  high  fidelities  when  using  neutral  atoms.  Current 
proposals  for  quantum  computing  focusing  on  neutral 
atoms  are  based  on  using  single  atoms  as  quantum  bits, 
often  while  using  cavity  induced  coupling  or  dipole-dipole 
interaction  for  two-qubit  operations.  However,  given  the 
degree  of  difficulties  encountered  in  isolating  and  control¬ 
ling  single  atoms,  this  process  has  proven  very  difficult  to 
realize,  especially  on  a  large  scale.  An  alternative  approach 
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is  to  use  atomic  ensembles  as  quantum  bits.  However, 
when  an  atomic  ensemble  is  excited,  by  a  laser  beam 
matched  to  a  two-level  transition  (or  a  Raman  transition) 
for  example,  it  leads  to  a  cascade  of  many  states  as  more 
and  more  photons  are  absorbed  [1-3].  In  order  to  make 
use  of  an  ensemble  as  a  qubit,  it  is  necessary  to  disrupt  this 
cascade,  and  restrict  the  excitation  to  the  absorption  (and 
emission)  of  a  single  photon  only.  In  principle,  this  can 
be  achieved  through  the  use  of  the  so-called  dipole  block¬ 
ade,  which  can  be  particularly  efficient  if  Rydberg  transi¬ 
tions  are  used  [4,5]. 

Dipole  blockades  generally  occur  between  individual 
atoms  within  an  ensemble.  In  order  to  make  use  of  this 
blockade  mechanism  in  a  manner  that  is  consistent  with 
a  quantum  computing  architecture,  it  is  necessary  to  con¬ 
trol  the  distribution  of  inter-atomic  distances  between  each 
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pair  of  atoms  in  the  ensemble  in  a  precise  manner.  Further¬ 
more,  in  order  to  achieve  long  decoherence  times,  it  is 
necessary  to  make  use  of  dipole  blockades  based  on  spin- 
spin  coupling,  which  is  necessarily  much  weaker  than  the 
optical  dipole-dipole  coupling. 

Here,  we  show  that  a  new  type  of  blockade  mechanism, 
based  on  the  light-shift  imbalance  in  a  Raman  transition, 
can  overcome  these  limitations.  The  resulting  system  does 
not  impose  any  constraint  on  the  distribution  of  inter-atomic 
distance  within  an  ensemble.  Furthermore,  no  dipole-dipole 
coupling  is  necessary,  so  that  a  relatively  low  density  system 
can  be  used. 

In  what  follows,  we  describe  first  a  simple  example  illus¬ 
trating  the  concept  of  light-shift  imbalance  induced  block¬ 
ade  (LSIIB)  using  a  multi-level  structure  in  a  single  atom, 
and  show  verifications  of  the  analytic  prediction  using 
numerical  simulations.  We  then  extend  this  model  to  show 
how  a  blockade  can  be  realized  by  using  LSI  in  the  excita¬ 
tion  of  an  ensemble.  Specifically,  we  show  how  the  LSIIB 
process  enables  one  to  treat  the  ensemble  as  a  two-level 
atom  that  undergoes  fully  deterministic  Rabi  oscillations 
between  two  collective  quantum  states,  while  suppressing 
excitations  of  higher  order  collective  states.  We  then  show 
how  this  transition  can  be  used  to  realize  a  quantum  bit 
(qubit)  embodied  by  the  ensemble.  Using  multiple  energy 
levels  inside  each  atom,  the  LSIIB  enables  the  transfer  of 
quantum  information  between  neighboring  ensembles,  as 
well  as  the  realization  of  a  C-NOT  gate.  In  effect,  this  rep¬ 
resents  a  generalization  of  the  so-called  Pellizari  scheme  for 
quantum  information  processing  [6].  Furthermore,  the 
LSIIB  can  be  used  to  link  two  separate  quantum  computers 
(QC),  by  transferring  the  quantum  state  of  any  ensemble 
qubit  in  one  QC  to  any  ensemble  qubit  in  another  QC. 
In  a  separate  paper,  we  discuss  details  of  these  quantum 
computation  and  communication  protocols,  offer 
practical  ways  to  implement  this  scheme,  and  propose 
experiments  to  demonstrate  the  feasibility  of  these  schemes 
[7]. 

The  significance  of  the  LSIIB  process  can  be  summa¬ 
rized  as  follows:  (a)  It  can  be  used  to  realize  a  determinis¬ 
tic  quantum  bit  encoded  in  the  collective  excitation  states 
of  an  atomic  ensemble,  (b)  Along  with  a  moderate-Q  cav¬ 
ity,  it  can  be  used  to  realize  a  two-qubit  gate  (e.g.,  a  C- 
NOT  gate)  between  two  ensemble-based  qubits,  (c)  It 
can  be  used  to  transport,  deterministically,  the  quantum 
state  of  an  ensemble  qubit  from  one  location  to  another 
separated  by  macroscopic  distances,  and  (d)  It  can  be  used 
to  establish  a  quantum-link  between  two  ensembles-and- 
cavity  based  quantum  computers.  The  scheme  proposed 
here  and  expanded  further  in  Ref.  [7]  therefore  offers  a 
robust  technique  for  realizing  a  quantum  network  without 
using  the  single  atom  and  super-cavity  based  approaches 
[8-10]. 

In  order  to  illustrate  the  basic  mechanism  that  underlies 
the  LSIIB,  it  is  convenient  to  consider  first  a  simple  exam¬ 
ple  of  a  set  of  five  levels  inside  a  single  atom,  as  shown  in 
Fig.  1.  Under  the  rotating  wave  approximation  (RWA), 


Fig.  1.  Schematic  illustration  of  a  five  level  system,  illustrating  the  process 
of  LSIIB. 


and  the  rotating  wave  transformation  [11]  (RWT),  the 
Hamiltonian  (with  h  =  1 )  describing  this  interaction  is 
given,  in  the  bases  of  {|1>,  |2>,  |3>,  |4>,  |5>},  by 
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where  we  have  defined  8  =  +  82)/2  and  A  =(8\  —  h2). 

Under  the  conditions  that  8^>  A,  8  Ul5  and  8  Q2, 

one  can  eliminate  the  optically  excited  states  |2>  and 
|4>  adiabatically.  The  effective  Hamiltonian  in  the  bases 
of  { 1 1>,  |3>,  1 5> }  is  then  given  by  [12-15] 
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where  UR  =  (QiQ2)/28  is  the  Raman-Rabi  frequency,  and 
Ej  =  Q^/48  is  the  light-shift  due  to  Q/j  =  1,  2).  Note  that 
the  levels  |1>,  |3>  and  |5>  are  light-shifted  by  different 
amounts.  In  general,  this  Hamiltonian  describes  a  process 
wherein  populations  can  oscillate  between  the  states  |1>, 
1 3>  and  1 5>,  with  the  maximum  amplitude  in  each  level  being 
determined  by  the  relative  values  of  the  parameter. 

Consider  now  the  case  where  Q2  -C  Ui.  Furthermore, 
assume  that  A  =  e2.  Under  this  condition,  the  Raman  cou¬ 
pling  between  |1>  and  |3>  become  resonant,  while  the 
Raman  coupling  between  |3>  and  |5>  becomes  detuned 
by  -(£1  +  £2).  Explicitly,  this  can  be  seen  by  subtracting 


an  energy  (£1 
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to  give 
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In  the  limit  of  e2  — >  0,  A  =  0,  the  detuning  for  the  |3>  to 
|5>  Raman  coupling  is  simply  —  gj.  For  UR  <C  £i  (which 
is  the  same  condition  as  02<cUi),  the  coupling  to 
level  |5>  can  be  ignored.  The  net  result  is  that  the  sys¬ 
tem  will  oscillate  between  |1>  and  |3>.  Note  that  this 
result  is  due  to  the  fact  that  levels  |1>  and  |3>  get 
light-shifted  by  nearly  the  same  amounts,  thus  remaining 
resonant  for  the  Raman  transition,  while  level  |5>  sees 
virtually  no  light-shift.  Thus,  the  excitation  to  level  |5> 
is  essentially  blockaded  by  the  imbalance  in  the  light- 
shifts. 

In  order  to  verify  this  blockade,  we  used  numerical 
techniques  to  integrate  the  equation  of  motion  of  the  sys¬ 
tem  shown  in  Fig.  1.  Since  each  of  the  optical  detunings 
(<51  and  52)  is  assumed  to  be  much  larger  than  the  Rabi 
frequencies  and  the  spontaneous  decay  rates,  the  effect  of 
the  decay  terms  is  not  significant  and  was  ignored  for 
simplicity.  The  population  of  each  of  the  levels  with  all 
the  atoms  starting  in  level  1,  is  shown  in  Fig.  2.  These 
are  calculated  for  the  following  parameter  values: 
A  =  0,  =0.1,  and  5/Q\  =  1 0.  Analytically,  the 

residual  populations  in  levels  |2>  and  |4>  are  expected 
to  be  of  the  order  of  {Q\/5)2  and  that  of  level  |5>  is 
expected  to  be  of  the  order  of  (Or/e!)2,  which  in  turn 
is  of  the  order  of  {Q2/Q\)2  and  are  consistent  with  the 
values  seen  here.  These  excitations  can  be  suppressed  fur¬ 
ther  by  making  these  ratios  smaller. 

Consider  next  the  excitation  of  an  ensemble  of  N  atoms. 
To  start  with,  we  consider  each  atom  in  the  ensemble  to  be 
a  three-level  system,  as  illustrated  in  Fig.  3.  Using  Dicke’s 
model  of  collective  excitation  [1],  we  can  show  that  the 
ensemble  excitation  can  be  represented  as  shown  in 
Fig.  4.  The  collective  states  in  this  diagram  are  defined  as 
follows: 


1  N 

I  A)  =  |  aua2,...,aN),  \Gi)=-j=Y2\aua2,...,gj,...aN), 

1  N 

\Ci)=-^='Y2\aua2,...,c],...aN), 

j  nc2 

lG2)=-^p=  Yl 

^  2  iW+k) 

N  f-' 

1  Cl 

|c2)  =  ...  /  ~  \d\  ,U2, . . .  .  . .  ,C£, . . .  Uy)  > 

v  CljMm 
j  7NC2 

|Gi4)=-^==  ^2  \a\ia2 


1  z 

|G2,i)  =~v=  ^ 

^  Z  C i+HV) 

1  z 

|Gi,2)ee—  V  \ai,a2,...,gj,...,ck,...,ch...aN)  [Z=  3jVC3] 


'  ( i+k+i) 


(4) 


The  relevant  coupling  rates  between  these  collective  states 
are  also  illustrated  in  Fig.  4.  Note  that  for  large  detunings, 
the  excitations  to  the  intermediate  states  |Gi>  and  |Gj . i > 
are  small,  so  that  higher  order  states  such  as  |G2>  and 
|G2ji>  can  be  ignored.  The  remaining  system  looks  very 
similar  to  the  five  level  system  considered  earlier.  (Paren¬ 
thetically  at  this  point,  note  that  the  coupling  between 
|Gi>  and  |Q>  does  not  scale  with  yjN,  unlike  the  coupling 
between  |A>  and  |Gi>,  which  scales  as  s/N).  We  can  now 
proceed  in  a  manner  similar  to  the  one  described  earlier. 
First,  in  the  rotating  wave  transformation  frame,  the  trun¬ 
cated,  six  level  Hamiltonian,  in  the  bases  of  |A>,  |Gi>, 


o  - 1 - I - - - — - 

0  m  IfflO  14M  2  (ID  2‘M 


Tine  t  [  n  IQ-1T1  inta  | 


Fig.  2.  Numerical  results  showing  populations  of  the  five  level  system,  starting  with  all  the  atoms  in  level  1. 
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Fig.  3.  Schematic  illustration  of  a  three-level  transition  in  each  atom  in  an 
ensemble. 


Fig.  4.  Schematic  illustration  of  the  relevant  collective  states  and  the 
corresponding  transition  rates. 


|Ci>,  |GU>,  |C  2>  and  |Gi;2>  is  given  by  (the  justification 
for  not  including  the  state  |C3>  will  be  made  by  showing 
that  the  excitation  to  |C2>  can  be  suppressed,  thus  in  turn 
making  the  amplitude  of  |C3>  insignificant) 


for  the  three  remaining  states  (|A>,  |Ci>,  and  |C2>)  are 
given  by  (assuming  5  »  A) 


EA  +  d  /  2 


H  = 


Qr0/2 


0 


Qro/2 

eci  —  d/2 

ff[2(N-l)/N]QRo/2 


0 

ff[2(N-l)/N]QKo/2 
£C2  —  3d/2 


(6) 


where  eA,  ec, ,  and  eCi  are  the  light-shifts  of  the  states  |A>, 
|Ci>,  and  |C2>,  respectively  and  URo  =  {ffNQ\Q2)/{28). 
To  first  order,  these  light-shifts  are  balanced,  in  the  sense 
that  (fiCi  —  £a)  =  (Ec2  —  £ci)-  This  means  that  if  the  explicit 
two-photon  detuning,  A,  is  chosen  to  make  the  Raman 
transition  between  |A>  and  |Ci>  resonant  (i.e., 
A  =  (eCl  —  eA)),  then  the  Raman  transition  between  |Ci> 
and  |C  2>  also  becomes  resonant.  This  balance  is  broken 
when  the  light-shifts  are  calculated  to  second  order  in 
intensity  (which  corresponds  to  fourth  order  in  Q)  ,  and 
the  blockade  shift  is  then  given  by 

dB  =  (SC2  —  £ci)  ~  (£ci  —  £a)  =  ~(&2  +  ) / (8<53 ) .  (7) 


With  the  proper  choice  of  two-photon  detuning 
(d  =  (eCi  —  £a))  to  make  the  Raman  transition  between 
|A>  and  |Ci>  resonant,  the  effective  Hamiltonian  (after 
shifting  the  zero  of  energy,  and  assuming  N  »  1 )  is  now  gi¬ 
ven  by 


H  = 


0 

Qr0/2 

0 


QRo/2  0 

0  ^2ro /\[2 

^ro/^/2  Ar 


(8) 


This  form  of  the  Hamiltonian  shows  clearly  that  when 
CdB,  the  coupling  to  the  state  |C2>  can  be  ignored. 
As  such,  the  collective  excitation  process  leads  to  a  Rabi 
oscillation  in  an  effectively  closed  two-level  system  consist¬ 
ing  of  |A>  and  |Ci>.  This  is  the  LSIIB  in  the  context  of 
ensemble  excitation,  and  is  the  key  result  of  this  paper. 

While  it  may  be  rather  obvious  at  this  point,  we  empha¬ 
size  nonetheless  that  we  can  now  represent  a  quantum  bit 
by  this  effectively  closed  two-level  system.  In  the  process, 
we  have  also  shown  how  to  perform  an  arbitrary  single 
qubit  rotation,  an  essential  pre-requisite  for  quantum  com- 


d/2  ffNQ{/2 
ffNQi/2  -S 
0  Q2/2 

0  0 

0  0 
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0 

Q2/2 

-A  12 

{VN^\)Ql/2 

0 

0 
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0 

{y/N-\)Qi/2 

-(<5  +  d) 
V2Q2/2 
0 


0  0 

0  0 

0  0 

V2Q2/2  0 

— 3d/2  {ffN~^2)Qff2 
~2)Ql/2  -(<5  +  2zl)  . 


(5) 


where  the  detunings  are  defined  just  as  before: 
5  =  (^3  +  S2)/2  and  A  =  (Si  —  S2). 

If  the  detunings  are  large  compared  to  the  transition 
rates,  we  can  eliminate  states  |Gi>,  |GU>  and  |GR2>  adi- 
abatically.  Under  this  condition,  the  effective  Hamiltonian 


puting.  The  details  of  how  such  a  qubit  can  be  used  for 
quantum  computation,  quantum  communication,  and  the 
realization  of  a  quantum  network  is  described  in  Ref.  [7], 
To  be  explicit,  let  us  consider  a  specific  numerical  exam¬ 
ple.  Choosing  the  natural  decay  rate,  T,  of  the  excited  state 
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to  be  unity,  assume  Q2  =  100,  kb  =  rj  <C  1,  and  8  =  1000. 
This  corresponds  to  a  value  of  A  =  bi  —  82  = 
£Ci  —  £a  —  2.5,  satisfying  the  condition  that  8)$>A.  The 
value  of  the  blockade  shift  is  AB  ^  —  1/80.  To  be  concrete, 
we  satisfy  the  requirement  that  (>Ro/\/'2  <c  AB  by  demand¬ 
ing  that  (Orq/^/2)  =  dB/10.  This  translates  to  the  condition 
that  rj^/N  —  0.035.  Note  that  these  parameters  satisfy  the 
constraint  that  5  »  Q2,  and  8  )$>  Q\^/N.  The  acceptable 
range  for  N  will  be  dictated  by  the  choice  of  rj,  or  vice- 
versa,  depending  on  the  particular  experiment  at  hand. 
As  an  example,  for  rj  =  10-3,  we  need  N=  1200.  To  see 
whether  such  a  range  of  parameters  are  potentially  realis¬ 
tic,  let  us  consider  an  ensemble  of  cold  S7Rb  atoms  caught 
in  a  FORT  trap,  excited  by  control  beams  with  a  cross-sec¬ 
tional  diameter  of  about  200  pm.  Assuming  that  the  transi¬ 
tions  employ  dipole  matrix  element  amplitudes  that  are 
half  as  strong  as  those  of  the  strongest  transitions,  the 
power  needed  for  the  Q2  beam  is  about  100  mW,  and  that 
for  the  beam  is  about  10  pW.  The  time  for  a  7t-transition 
going  from  |A>  to  |Ci>  is  about  50  ps.  Given  that  the 
decoherence  time  in  a  FORT  can  be  of  the  order  of  min¬ 
utes,  as  many  as  106  qubit  operations  can  be  carried  out 
at  this  rate.  The  number  of  photons  in  the  Qi  is  close  to 
2000,  so  that  its  treatment  as  a  classical  beam  is  valid. 

The  essence  of  the  LSI1B  for  ensemble  excitation  is  sum¬ 
marized  in  Fig.  5.  Briefly,  whenever  we  have  a  three-level 
optically  off-resonant  transition  for  the  individual  atoms, 
this  can  be  translated  into  a  corresponding  off-resonant 
three-level  transition  involving  collective  states,  which  in 
turn  is  reduced  to  an  effective  two-level  transition.  In  order 
for  this  to  hold,  the  primary  constraint  is  that,  for  the  col¬ 
lective  states,  the  Rabi  frequency  on  one  leg  must  be  much 
bigger  than  the  same  for  the  other. 

To  summarize,  we  have  described  a  new  type  of  block¬ 
ade  that  allows  one  to  treat  an  ensemble  excitation  as  a  sin¬ 


gle,  deterministic  quantum  bit  consisting  of  only  two-levels. 
Such  a  system  can  be  used  to  realize  a  two-qubit  gate  (e.g., 
a  C-NOT  gate)  between  two  ensemble-based  qubits.  It  can 
also  be  used  to  transport,  deterministically,  the  quantum 
state  of  an  ensemble  qubit  from  one  location  to  another 
separated  by  macroscopic  distances,  and  it  can  be  used  to 
establish  a  quantum-link  between  two  ensembles-and-cav- 
ity  based  quantum  computers. 
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Quantum  communication  and  computing  with  atomic  ensembles 
using  a  light-shift-imbalance-induced  blockade 
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Recently,  we  have  shown  that  for  conditions  under  which  the  so-called  light-shift  imbalance  induced  block¬ 
ade  occurs,  the  collective  excitation  of  an  ensemble  of  a  multilevel  atom  can  be  treated  as  a  closed  two-level 
system.  In  this  paper,  we  describe  how  such  a  system  can  be  used  as  a  quantum  bit  (qubit)  for  quantum 
communication  and  quantum  computing.  Specifically,  we  show  how  to  realize  a  controlled-NOT  gate  using  the 
collective  qubit  and  an  easily  accessible  ring  cavity,  via  an  extension  of  the  so-called  Pellizzari  scheme.  We 
also  describe  how  multiple,  small-scale  quantum  computers  realized  using  these  qubits  can  be  linked  effec¬ 
tively  for  implementing  a  quantum  internet.  We  describe  the  details  of  the  energy  levels  and  transitions  in  an 
87Rb  atom  that  could  be  used  for  implementing  these  schemes. 
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I.  INTRODUCTION 

When  an  atomic  ensemble  is  excited,  by  a  laser  beam 
matched  to  a  two-level  transition  (or  a  Raman  transition),  for 
example,  it  leads  to  a  cascade  of  many  states  as  more  and 
more  photons  are  absorbed  [1-3],  In  order  to  make  use  of  an 
ensemble  as  a  quantum  bit  (qubit),  it  is  necessary  to  disrupt 
this  cascade,  and  restrict  the  excitation  to  the  absorption  (and 
emission)  of  a  single  photon  only.  In  principle,  this  can  be 
achieved  through  the  use  of  the  so-called  dipole  blockade 
[4,5],  In  order  to  make  use  of  this  blockade  mechanism  in  a 
manner  that  is  consistent  with  a  quantum  computing  archi¬ 
tecture,  it  is  necessary  to  control  the  distribution  of  inter¬ 
atomic  distances  between  each  pair  of  atoms  in  the  ensemble 
in  a  precise  manner  Furthermore,  in  order  to  achieve  long 
decoherence  times,  it  is  necessary  to  make  use  of  dipole 
blockade  based  on  spin-spin  coupling,  which  is  necessarily 
much  weaker  than  the  optical  dipole-dipole  coupling.  Re¬ 
cently,  we  have  shown  that  a  different  type  of  blockade 
mechanism,  based  on  the  light-shift  imbalance  in  a  Raman 
transition,  can  overcome  these  constraints.  The  resulting  sys¬ 
tem  does  not  impose  any  constraint  on  the  distribution  of 
interatomic  distance  within  an  ensemble.  Furthermore,  no 
dipole-dipole  coupling  is  necessary,  so  that  a  relatively  low 
density  system  can  be  used. 

In  Ref.  [6],  we  have  shown  how  this  light-shift- 
imbalance-induced  blockade  (LSIIB)  process  enables  one  to 
treat  the  ensemble  as  a  two-level  atom  that  undergoes  fully 
deterministic  Rabi  oscillations  between  two  collective  quan¬ 
tum  states,  while  suppressing  excitations  of  higher-order  col¬ 
lective  states.  In  this  paper,  we  show  how  this  transition  can 
be  used  to  realize  a  qubit  embodied  by  the  ensemble.  Using 
multiple  energy  levels  inside  each  atom,  we  show  how  the 
LSIIB  enables  the  transfer  of  quantum  information  between 
neighboring  ensembles,  as  well  as  the  realization  of  a 
controlled-NOT  gate.  In  effect,  this  represents  a  generaliza¬ 
tion  of  the  so-called  Pellizzari  scheme  for  quantum  informa¬ 
tion  processing  [7].  Furthermore,  we  show  that  the  LSIIB 
can  be  used  to  link  two  separate  quantum  computers  (QCs), 
by  transferring  the  quantum  state  of  any  ensemble  qubit  in 
one  QC  to  any  ensemble  qubit  in  another  QC.  We  discuss 


practical  ways  to  implement  this  scheme  using  specific  en¬ 
ergy  levels  and  transitions  in  87Rb  atoms,  and  propose  ex¬ 
periments  to  demonstrate  the  feasibility  of  these  protocols. 

The  significance  of  the  light-shift-imbalance-induced 
blockade  (LSIIB)  process  can  be  summarized  as  follows:  (i) 
It  can  be  used  to  realize  a  deterministic  quantum  bit  encoded 
in  the  collective-excitation  states  of  an  atomic  ensemble,  (ii) 
Along  with  an  easily  accessible  ring  cavity,  it  can  be  used  to 
realize  a  two-qubit  gate  (e.g.,  a  controlled-NOT  gate)  between 
two  ensemble-based  qubits,  (iii)  It  can  be  used  to  transport, 
deterministically,  the  quantum  state  of  an  ensemble  qubit 
from  one  location  to  another  separated  by  macroscopic  dis¬ 
tances.  (iv)  It  can  be  used  to  establish  a  quantum  link  be¬ 
tween  two  ensembles-and-cavity  based  quantum  computers. 
The  scheme  proposed  here  thus  offers  a  robust  technique  for 
realizing  a  quantum  internet  without  using  the  single-atom 
and  ultrashort  cavity  based  approaches  [8-10]. 


II.  FUNDAMENTALS  OF  THE  LSIIB  PROCESS 

We  summarize  first  the  basic  notations  for  describing  the 
LSIIB  process  [6].  In  the  simplest  version,  we  consider  each 
atom  in  the  ensemble  to  be  a  A-type  three-level  system,  with 
two  metastable  states  (|a)  and  |c)),  and  an  optically  excited 
state  (|e)),  as  shown  in  Fig.  1.  The  collective  states  of  interest 
are  defined  as  follows  [6]: 


|A)  =  | al,a2,  ...  ,aN), 


|Gi)  =  ~t=E  \ai,a2,  ■  ■  ■  ,gj,  •••<%>, 
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FIG.  1.  Schematic  illustration  of  a  three-level  transition  in  each 
atom  in  an  ensemble. 
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where  nCm=(^i)  =  NI/\M\(M-N)\ ].  The  relevant  coupling 
rates  between  these  collective  states  are  also  illustrated  in 
Fig.  2.  Note  that  for  large  detunings,  the  excitations  to  the 
intermediate  states  |G[)  and  | G u )  are  small,  so  that  higher- 
order  states  such  as  |G2)  and  | G2  1)  can  be  ignored.  (Paren¬ 
thetically  at  this  point,  note  that  the  coupling  between  |  G 1 ) 
and  |C3)  does  not  scale  with  JN,  unlike  the  coupling  between 
| A)  and  |GX),  which  scales  as  JN.) 

In  the  rotating  wave  transformation  frame,  the  truncated, 
six-level  Hamiltonian,  in  the  bases  of  |A),  |Gi),  |Ci),  |G[ j), 
|C2),  and  |G1i2)  is  given  by  (the  justification  for  not  including 
the  state  |C3)  will  be  made  by  showing  that  the  excitation  to 
|C2)  can  be  suppressed,  thus  in  turn  making  the  amplitude  of 
|C3)  insignificant) 


H  = 


A/2 

\,YH/2 

0 

0 

0 

0 


\7Vfl  |  /2  0 

-  5  H2/2 

H2/2  -  A/2 

0  (V/v-  m.n 

0  0 

0  0 


0 

0 

(\W-  l)fl,/2 
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V2Hz/2 
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0 

0 

0 

V2H2/2 

-3A/2 

(V/v-  2)n1/2 


0 

0 

0 

0 

(\.Y- 2)12/2 
-(<5+ 2  A) 


(2a) 


where  the  detunings  are  defined  as  8=  (<5,  +  82)!2  and  A 
=  ( -  <%)  ■ 

If  the  detunings  are  large  compared  to  the  transition  rates, 
we  can  eliminate  states  | G3),  |G3 j)  and  |G3  2)  adiabatically. 
Under  this  condition,  the  effective  Hamiltonian  for  the  three 
remaining  states  (|A),  |Cj),  and  |Co))  are  given  by  (assuming 
S>  A) 


H  = 


ea  +  A/2  f  Iav,/2 

^Ro/2  ecl  -  A/2 

0  \[2(N-  l)/lV]fift/2 


0 

V[2(/V-  l)/N\tlgJ2 

eci  ~  3A/2 


(2b) 


first  order,  these  light  shifts  are  given  by  £A=ND,2l/4S,  ec 
=  [03  +  (N-l)nj]/45,  and  EC=[m\  +  (N-2)ti\\l AS,  and 
are  balanced,  in  the  sense  that  (£cl-£A)  =  (£C2-ecl).  This 
means  that  if  the  explicit  two-photon  detuning  A  is  chosen  to 
make  the  Raman  transition  between  |A)  and  |Cj)  resonant 
[i.e.,  A  =  (sc1-£A)],  then  the  Raman  transition  between  | Cj) 
and  |C2)  also  becomes  resonant.  This  balance  is  broken  when 
the  light  shifts  are  calculated  to  second  order,  and  the  block¬ 
ade  shift  is  then  given  by 


A b  —  (Bc2  ~  eci)  -  (eci  -  sa)  -  ~  (^2  +  n{)/(8<53) .  (3) 


where  sA,  ec1,  and  sc2  are  the  light  shifts_of  the  states  |A),  With  the  proper  choice  of  two-photon  detuning  [A  =  (scl 
| C?!),  and  | Co),  respectively,  and  12 Ko  =  ( y  Ml |  ll2)  /  (2  <5) .  To  —  ea)]  to  make  the  Raman  transition  between  |A)  and  IC^ 
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FIG.  2.  Schematic  illustration  of  the  relevant 
collective  states  and  the  corresponding  transition 
rates.  See  text  for  details. 


resonant,  the  effective  Hamiltonian  (after  shifting  the  zero  of 
energy,  and  assuming  N9>  1)  is  now  given  by 


0 


H  = 


nRo/2 


o 


2  0 
0  _  fwVS 
nSo/V 2  ab 


(4) 


This  form  of  the  Hamiltonian  shows  clearly  that  when 
[fl*,/  \  2]  Aft,  the  coupling  to  the  state  |C2)  can  be  ignored. 
As  such,  the  collective  excitation  process  leads  to  a  Rabi 
oscillation  in  an  effectively  closed  two-level  system  consist¬ 
ing  of  | A)  and  |Cj).  This  is  the  LSIIB  in  the  context  of  en¬ 
semble  excitation,  and  is  the  key  result  upon  which  most  of 
the  protocols  we  propose  here  are  based. 

Thus  we  can  now  represent  a  quantum  bit  by  this  effec¬ 
tively  closed  two-level  system.  In  the  process,  we  have  also 
shown  how  to  perform  an  arbitrary  single  qubit  rotation,  an 
essential  prerequisite  for  quantum  computing.  The  essence  of 
the  LSIIB  for  ensemble  excitation  can  be  summarized  as 
follows.  Whenever  we  have  a  three-level  optically  off- 
resonant  transition  for  the  individual  atoms,  this  can  be  trans¬ 
lated  into  a  corresponding  off-resonant  three-level  transition 
involving  collective  states,  which  in  turn  is  reduced  to  an 
effective  two-level  transition,  as  illustrated  in  Fig.  3.  In  order 
for  this  to  hold,  the  primary  constraint  is  that,  for  the  collec¬ 
tive  states,  the  Rabi  frequency  on  one  leg  must  be  much 
bigger  than  the  same  for  the  other. 

At  this  point,  it  is  instructive  to  consider  a  specific  nu¬ 
merical  example.  Choosing  the  natural  decay  rate  T  of  the 
excited  state  to  be  unity,  assume  02  =  1 00,  11,  =  77  1 ,  and 

<5=1000.  This  corresponds  to  a  value  of  A  =  A ,  —  S2  =  cCi  -cA 
=  2.5,  satisfying  the  condition  that  A.  The  value  of  the 
blockade  shift  is  ABss-J/80.  To  be  concrete,  we  satisfy  the 
requirement  that  ClRo/y2<  AB  by  demanding  that  (£lRo/  J 2) 
=  Ab/10.  This  translates  to  the  condition  that  r],  /V=0.035. 


Note  that  these  parameters  satisfy  the  constraint  that  S>  fl2> 
and  The  acceptable  range  for  N  will  be  dictated 

by  the  choice  of  77,  or  vice  versa,  depending  on  the  particular 
experiment  at  hand.  As  an  example,  for  y=  10-3,  we  need 
N=  1200.  To  see  whether  such  a  range  of  parameters  is  po¬ 
tentially  realistic,  let  us  consider  an  ensemble  of  cold  87Rb 
atoms  caught  in  a  far  off-resonant  trap  (FORT)  trap,  excited 
by  control  beams  with  a  cross-sectional  diameter  of  about 
200  ptm.  Assuming  that  the  transitions  employed  dipole  ma¬ 
trix  element  amplitudes  that  are  half  as  strong  as  those  of  the 
strongest  transitions,  the  power  needed  for  the  fl2  beam  is 
about  100  mW,  and  that  for  the  11,  beam  is  about  10  pW. 
The  time  for  a  tt  transition  going  from  |A)  to  |Ci)  is  about 
50  /lls.  Given  that  the  decoherence  time  in  a  FORT  can  be  of 
the  order  of  minutes,  as  many  as  106  qubit  operations  can  be 
carried  out  at  this  rate.  The  number  of  photons  in  the  11,  is 
close  to  2000,  so  that  its  treatment  as  a  classical  beam  is 
valid. 

HI.  REALIZATION  OF  A  CONTROLLED-not  GATE 
BETWEEN  TWO  ENSEMBLES 


The  LSIIB  process  can  be  used  to  perform  two  qubit  op¬ 
erations  between  two  ensembles.  In  particular,  we  show  here 


FIG.  3.  Summary  of  the  LSIIB  process  in  an  ensemble.  See  text 
for  details. 
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FIG.  4.  Illustration  of  the  basic  configuration  for  coupling  two 
ensembles  (top),  and  the  requisite  energy  levels  for  each  atom  (bot¬ 
tom).  See  text  for  details. 


FIG.  5.  (a)  Schematic  illustration  of  the  single  qubit  state  prepa¬ 
ration  in  ensemble  I.  (b)  Schematic  illustration  for  the  single  qubit 
operation  for  ensemble  II.  See  text  for  details. 


how  a  controlled-NOT  gate  can  be  realized  in  such  a  system. 
The  process  we  describe  is  essentially  the  ensemble  version 
of  the  well-known  Pellizzari  scheme  [7].  Figure  4  (top) 
shows  the  basic  configuration  involving  only  two  ensembles 
and  a  cavity.  A  more  general  scheme  that  can  couple  the 
nearest-neighbor  ensembles  in  a  large  array  of  ensembles  is 
described  schematically  later  on.  In  what  follows,  the  two 
ensembles  will  be  denoted  as  ensembles  I  and  II,  respec¬ 
tively.  Figure  4  (bottom)  shows  the  necessary  energy  levels 
in  each  of  the  atoms.  The  transition  corresponding  to  the 
cavity  mode  is  also  illustrated,  with  a  vacuum  Rabi  fre¬ 
quency  of  gc.  Later  on,  we  will  discuss  how  an  alkali-metal 
atom  such  as  87Rb  may  be  used  to  implement  these  energy 
levels. 

We  assume  that  both  ensembles  are  prepared  in  the  |A) 
state.  For  the  other  states,  the  only  ones  that  will  be  relevant 
are  the  ones  that  are  of  the  following  form: 


IGi) 


\In 


N 

2  \ai,a2,  ...aN), 

j=  i 


q  =  s,b,c,d,e,f,g  or  h  (excluding  a) .  (5) 

We  now  describe  the  controlled-NOT  gate  operation  in  mul¬ 
tiple  steps. 

Step  1:  Initialize  each  ensemble  in  an  arbitrary  quantum 
superposition .  At  1=0,  we  assume  that  each  ensemble  is  in 
state  |A):|i/f)/=|A)/;  |  «/r)7/=  | A)n.  This  can  be  achieved,  in 
principle,  by  optical  pumping  in  a  real  system.  Operations 
for  ensemble  1:  The  transitions  to  be  used  for  performing 
the  single  qubit  rotation  is  shown  in  Fig.  5.  With 
we  apply  these  two  pulses  simultaneously  for  duration 
7j.  The  state  of  ensemble  I  (E-I)  is  then  given  by 
|  if/\ = cosLQ^r,  /  2]|A)I+ isintn^r,  /  2]|  Cj)!  -  a\A\  +  0\  Cx\. 


Here,  0A>„  is  as  defined  in  Eq.  (2)  above.  This  achieves  the 
goal  of  producing  an  arbitrary  quantum  state  in  E-I.  Opera¬ 
tions  for  ensemble  II:  For  the  second  ensemble,  we  apply  the 
same  transitions  as  in  Fig.  5(a),  but  for  a  duration  T2  such 
that  HRoT2=tt.  This  corresponds  to  a  pi  pulse,  and  ensemble 
II  (E-II)  is  now  in  state  |Ci).  This  is  now  followed  by  the 
excitation  pulses  shown  in  Fig.  5(b).  Note  in  particular  that 
the  coupling  rates  for  the  collective  states  no  longer  include 
the  v/V  factor.  The  pulses  at  (l  j  and  Oj  are  applied  for  a 
duration  T3  (satisfying  the  condition  that  il'2<  il\),  so  that 
the  state  for  E-II  is  now  given  by  |(/f)II=cos[n,7’3/2]|C)n 
+  fsin[n,7’1/2]|Z))II=  £|C)n+  rj\D)n.  Thus  we  now  have 
prepared  E-II  in  an  arbitrary  qubit  state.  Finally,  note  that 
the  cavity  so  far  is  in  the  ground  state  with  no  photons: 
l'A)c=|0)c- 

Step  2:  Transfer  state  of  E-I  to  the  cavity.  The  transitions 
to  be  employed  for  this  step  is  illustrated  in  Fig.  6(a).  The 
pulse  at  fl7  is  applied  for  a  duration  74  such  that 
UN£llgcl2S)  74  =tt.  As  a  result  of  this  pi  pulse,  the  states  of 
E-II  remain  unchanged,  but  the  states  of  E-I  and  the  cavity 
become  as  follows:  |</f)c=aj0)c+/?|l)c;  \if/)I=\Cl)I.  This  ac¬ 
complishes  the  goal  of  transferring  the  quantum  state  of  E-I 
to  the  cavity. 

Step  3:  Transfer  state  of  the  cavity  to  E-II.  The  transi¬ 
tions  used  for  this  step,  in  E-II,  is  illustrated  in  Fig.  6(b). 
Note  again  that  the  transition  rates  for  coupling  the  collective 
states  now  do  not  include  the  -JN  factor.  The  pulse  at  ilu  is 
applied  for  a  duration  T5  such  that  (ilngc/2S)  T5=tt.  As  a 
result  of  this  pi  pulse,  the  states  of  E-I  remains  unchanged, 
but  the  states  of  E-II  and  the  cavity  become  as  follows: 
l'A)c=|0)c;  |<A)//=a£|Si)//+£/3|C1)//+  yP\D\)n.  This 

accomplishes  the  goal  of  transferring  the  quantum  informa¬ 
tion  from  the  cavity  to  E-II.  Of  course,  steps  2  and  3  together 
achieve  the  objective  of  transferring  all  the  quantum  infor- 
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FIG.  6.  (a)  Illustration  of  the  excitation  step  used  to  transfer  the 
state  of  E-I  to  the  cavity,  (b)  Schematic  illustration  of  the  pulses 
inside  E-II  to  transform  the  cavity  state  to  E-II.  See  text  for  details. 

mation  from  both  ensembles  into  ensemble  II  only. 

Step  4:  Perform  “effective  controlled-NOT  transition”  in¬ 
side  E-II.  Using  the  transitions  shown  in  Fig.  5(b),  we  apply 
a  pi  pulse  to  E-II,  which  exchanges  the  amplitudes  of  |C|)n 
and  |ZD1)II.  As  a  result,  while  E-I  and  the  cavity  stays  un¬ 
changed,  the  state  of  E-II  now  becomes  \f)n=at;\Sl)n 
+  &\D\ )//+  va\Bx)u+  rj(3\Cfu. 

Step  5:  Transfer  E-II  quantum  state  partly  to  the  cavity. 
Using  the  transitions  shown  in  Fig.  6(b),  we  now  apply  an¬ 
other  pi  pulse  in  E-II.  As  a  result,  while  E-I  stays  unchanged, 
the  cavity  and  the  E-II  now  become  entangled:  |  t//)n.c 
=  a;^|Ci)//|l)c+^/3|E)i)//|0)c+  rja\Dx)n\l)c+  ?7^|C1)//|0)c. 

Step  6:  Transfer  cavity  state  back  to  E-I.  Using  the  tran¬ 
sitions  shown  in  Fig.  3(a),  we  now  apply  another  pi  pulse  in 
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E-I.  As  a  result,  the  cavity  returns  to  the  zero-photon  state, 
and  E-I  and  E-II  are  now  entangled:  |(/f)c=|0)c>  l*A)/-// 
=  a^\A)I\Cl)II+  ?7a:|A)/|E)1)//+^|C1)/|Z)1)//+^77|C1)/|C1)//. 
This  is  to  be  compared  with  the  direct-product  state  at  the 
beginning  of  the  protocol:  \4di-  /Xinit)  =  |  if)  n 

=  a^\A)j\Cl)n+ay\A)I\D1)n+f3^\Cl)I\Cl)II+f3r]\Cl)I\Dl)n. 
Thus  the  quantum  state  of  E-II  has  been  flipped  only  for  the 
case  where  the  quantum  state  of  E-I  is  in  |Q),  which,  of 
course,  corresponds  to  the  controlled-NOT  operation. 

IV.  QUANTUM  COMMUNICATION  BETWEEN 
LSIIB-BASED  QUANTUM  COMPUTERS 

The  whole  operation  for  realizing  interprocessor  commu¬ 
nication  is  illustrated  in  Fig.  7.  Briefly,  each  QC  is  modeled 
as  an  array  of  ensembles  that  can  be  moved  linearly,  using 
several  different  possibilities,  including  quantum  motors  or 
sliding  standing  waves  through  a  cavity  [11-15].  It  is  also 
possible  to  consider  an  implementation  where  race-track- 
shaped  microcavities  are  used  on  alternating  sides,  with  the 
atom  seeing  the  cavity  mode  through  evanescent  modes. 
Note  also  that  while  we  are  concentrating  primarily  on  en¬ 
sembles  of  neutral  atoms  here,  there  are  solid  material  and 
systems  that  display  the  kind  of  narrow  Raman  transitions 
necessary  for  implementing  such  a  scheme  [16-18],  although 
we  have  not  yet  worked  out  the  details  for  an  explicit  solid- 
state  system. 

The  objective  of  the  communication  link  is  to  transfer  the 
quantum  state  of  the  ensemble  (Ql)  on  the  edge  of  one  QC 
(QC1)  to  the  ensemble  (Q2)  on  the  edge  of  the  other  QC 
(QC2).  To  start  with,  assume  that  Ql  is  in  an  arbitrary  quan¬ 
tum  state  (the  protocol  would  work  just  as  well  in  the  more 
useful  case  when  this  atom  is  entangled  with  the  rest  of 
QC1),  and  Q2  is  in  its  ground  state:  \4>)qi  =  oc\A)qi 
+P\C\)q\,  \>!>)q2  =  \A)q2.  The  transitions  to  be  used  for  the 
transfer  process  are  illustrated  in  the  bottom  of  Fig.  7.  Note 
that  here  we  are  using  the  same  set  of  transition  as  in  Fig. 
5(a),  except  that  instead  of  using  a  semiclassical  mode  to 


FIG.  7.  Schematic  illustration  of  quantum 
communication  between  two  quantum  computer. 
See  text  for  details. 
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couple  the  |g)  to  |c)  transition,  we  use  a  free-space  single¬ 
photon  mode.  The  parameters  of  this  mode  (such  as  the 
vacuum  Rabi  frequency  gF)  is  determined  essentially  by  the 
pump  pulse _[3, 18].  We  now  apply  a  pulse  for  a  duration  tr 
such  that  ( iN£l\CgFl28 )  tr=tt.  The  result  would  be  a  trans¬ 
fer  of  the  quantum  state  of  Q1  to  the  free-space  photon  mode 
(in  essentially  the  same  way  as  in  step  2  of  the  controlled- 
NOT  protocol).  A  fiber  delay  line  will  allow  one  to  let  this 
photon  state  leave  QC1  completely.  The  pump  pulse  on  Q2 
will  be  timed  with  the  anticipated  arrival  of  this  photon  state 
(this  can  be  done,  since  the  delay  between  the  photon  state 
and  the  pump  in  Q1  can  be  precalibrated  as  well  as  calcu¬ 
lated).  This  pump  will  then  be  applied  for  a  duration  rw  such 
that  (\M22g?f/2<5)  tw=tt.  Note  that  the  presence  of  the  JN 
factor  makes  both  of  these  processes  fast  enough  even  for  a 
free-space  photon  mode.  After  the  second  pulse,  the  free- 
space  photon  mode  will  be  in  the  zero  photon  state,  and  the 
quantum  states  of  Q1  and  Q2  would  be  i/y)y2=  ^\^)q2 
+/8|C1>fi2;  |(/f)Gi  =  |A)gi.  It  is  obvious  that  if  Q1  were  en¬ 
tangled  with  the  rest  of  QC 1 .  then  Q2  will  now  be  entangled 
with  the  rest  of  QC1.  This  is,  of  course,  all  that  is  needed  to 
link  the  two  quantum  computers. 

V.  EXPERIMENTAL  CONSTRAINTS  ON  CAVITY 
AND  ENSEMBLE  PARAMETERS 

In  order  to  see  the  feasibility  of  these  operations,  it  is 
important  to  consider  the  relevant  cavity  parameters.  When 
the  ensemble  qubit  is  interacting  with  the  cavity  mode  (e.g., 
in  step  2  above),  the  vacuum  coupling  rate  is  enhanced  by  a 
factor  of  \N,  so  that  the  cavity  parameters  for  reaching  the 
strong-coupling  regime  is  significantly  relaxed.  To  illustrate, 
consider  the  situation  where  a  cavity  holds  an  ensemble  ru¬ 
bidium  qubit  of  N  atoms,  with  a  length  of  L,  and  an  effective 
mode  diameter  of  d  at  the  center.  The  single-photon  electric- 
field  amplitude  is  then  given  by  yJ(2fuo)/(s0V)  where  V 
—  {tt! A)d2L  is  the  effective  mode  volume  [19],  and  co  is  the 
frequency  of  the  photon.  For  the  rubidium  atom,  with  N=  1 , 
the  corresponding  vacuum  rabi  frequency  is  proportional  to 
this  field  amplitude,  given  by  ga  —  54.25T  (where  T  —  2 tt6 
X  106  s_1  is  the  natural  line  width)  for  Lo=40  /jura  and  da 
=  5  /Jim  [20].  If  we  assume  a  matched  pair  of  essentially 
lossless  mirrors,  each  with  an  intensity  transmittivity  of  1.2 
X  10-6,  the  corresponding  cavity  finesse  is  about  2.6  X  106, 
with  a  free  spectral  range  (FSR)  of  about  3.8  X  1012  Hz  [21]. 
This  yields  a  cavity  decay  rate,  y0  (half  width  at  half  maxi¬ 
mum),  of  about  0.12r,  so  that  the  cavity  lifetime,  r0=  y~\  is 
about  222  ns.  From  the  discussion  above,  it  is  easy  to  see 
that  these  parameters  scale  as  follows:  g=g0(d0/ d) 
X(L0/L)mN112,  and  t=t0(T0IT)(LI L0). 

Consider  next  a  situation  where  the  number  of  atoms  in 
the  collective  qubit  is  3000.  Using  the  scaling  laws  stated 
above,  we  see  that  the  vacuum  rabi  frequency  will  remain 
unchanged  if  the  length  is  increased  by  a  factor  of  3000  ( L 
=  12  cm).  Here,  we  consider  a  cavity  length  of  L  =  5  cm,  so 
that  the  vacuum  rabi  frequency  is  ga—  84.04r,  and  the  FSR 
is  3.0  X  109  Hz.  For  the  same  finesse  as  above,  the  cavity 
decay  rate  is  now  only  9.55  X  10-5  F,  with  the  corresponding 
cavity  lifetime  being  close  to  0.3  ms.  With  expected  im¬ 


provements  in  the  mirror  qualities,  the  cavity  storage  time 
may  possibly  become  even  longer.  If  we  choose  the  strength 
of  to  be  7  X  Kr4r,  with  a  common-mode  detuning  of  S 
=  103r,  the  constraints  for  the  LSIIB  are  easily  satisfied.  The 
time  needed  for  a  tt  pulse  is  about  50  /U.s,  which  is  consid¬ 
erably  shorter  than  the  cavity  lifetime.  In  this  context,  we 
point  out  that  a  variant  of  the  Pellizari  scheme  [7] — which 
makes  use  of  the  so-called  Cavity  Dark  State  [18] — can  also 
be  used  to  implement  a  controlled-NOT  gate.  When  this  ap¬ 
proach  is  used,  the  cavity  remains  virtually  unpopulated  dur¬ 
ing  the  operation,  thereby  reducing  the  constraint  on  the  cav¬ 
ity  lifetime.  Another  important  consequence  of  such  a  long 
cavity  is  that  it  allows  one  to  use  a  ring-configuration  rather 
than  a  conventional  Fabry-Perot  configuration,  thereby 
avoiding  the  formation  of  standing  waves  in  the  cavity.  The 
long  cavity  also  makes  it  potentially  easier  to  access  the  en¬ 
semble  with  the  external  control  beams. 

VI.  PROPOSED  EXPERIMENTS  USING  RUBIDIUM 
ATOMIC  ENSEMBLES 

The  energy  levels  and  selection  rules  necessary  for  imple¬ 
menting  these  schemes  can  be  accommodated  by  using,  for 
example,  Zeeman  sublevels  in  the  D1  transition  manifold  in 
87Rb  atoms.  One  explicit  construct  is  shown  in  Fig.  8.  A 
nonzero  magnetic  field  is  used  to  lift  the  degeneracy  between 
the  state  |a)  and  states  {|j)  and  |b)}.  This  is  necessary  in  order 
to  ensure,  for  example,  that  during  the  process  outlined  in 
Fig.  5(a),  there  is  no  Raman  coupling  between  |,v)  and  |c)  (or 
between  | b)  and  \d)).  Similar  constraints  also  apply  for  pro¬ 
cesses  illustrated  in  Fig.  6(a).  For  the  process  in  Fig.  6(b),  the 
opposite  constraint  applies  (that  is,  |s)  to  |c)  as  well  as  \b)  to 
| d)  are  allowed,  but  |a)  is  not  coupled  due  to  Raman  detun¬ 
ing).  The  amount  of  magnetic  field  to  be  applied  can  be 
rather  small,  (typically  less  than  a  gauss),  with  the  constraint 
that  the  Zeeman  shift  should  be  larger  than  the  Raman  tran¬ 
sition  linewidths.  The  second-order  Zeeman  shift  for  the  |o) 
state  would  be  very  small  in  this  case,  and  can  be  taken  into 
account  in  choosing  the  relevant  Raman  resonance  condi¬ 
tions. 

Note  that  the  g  factors  for  the  F=  1  and  F=  2  levels  are 
equal  and  opposite,  which  ensures  that  the  s-c  transition  and 
the  b-d  transitions  can  be  simultaneously  Raman  resonant. 
Note  also  that  the  g  factors  for  the  F'  =  1  and  F'  =  2  are  also 
equal  and  opposite  to  each  other,  and  are  a  factor  of  a  3 
smaller  than  those  for  the  F=  1  and  F= 2  levels.  More  impor¬ 
tantly,  because  the  Raman  transitions  are  far  detuned  from 
the  F'  sublevels  optically,  the  Zeeman  shifts  of  the  sublevels 
for  F'  =  1  and  F'=  2  are  of  virtually  no  significance. 

The  energy  difference  between  the  F' - 1  and  F' -2  is 
more  than  816  MHz.  Thus  for  processes  that  are  detuned 
below  resonance  with  respect  to  the  F'  =  1  [e.g.,  for  the  pro¬ 
cess  in  Fig.  6(b)]  sublevel  will  not  be  influenced  significantly 
by  the  F'  =  2  level.  By  the  same  token,  for  processes  that  are 
detuned  above  resonance  with  respect  to  the  F'=  2  [e.g.,  for 
the  process  in  Fig.  5(a)]  sublevel  will  not  be  influenced  sig¬ 
nificantly  by  the  F'  - 1  level.  This  consideration  was  taken 
into  account  in  designing  the  protocols,  so  that  for  some 
operations  the  detuning  is  negative  [e.g.,  Fig.  2(a)],  while  for 
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FIG.  8.  (a)  Schematic  illustration  of  the  Zeeman  sublevels  in  the  D1  manifold  of  87Rb  atoms  that  can  be  used  to  implement  the  proposed 
schemes,  (b)  Schematic  illustration  of  the  87Rb  transitions  that  can  be  used  for  demonstrating  the  LSIIB  effect  in  a  single  cluster  without  a 
surrounding  cavity.  See  text  for  details. 


others  the  detuning  is  positive  [e.g..  Fig.  3(b)]. 

Finally,  note  that  in  order  to  allow  for  the  selection  rules 
assumed  in  the  protocols,  it  is  necessary  to  use  a  combination 
of  tr+,  <x_,  and  tt  polarized  light.  This  can  be  easily  accom¬ 
modated,  since  two  orthogonal  directions  are  free  for  apply¬ 
ing  classical  field  even  in  the  presence  of  a  cavity. 

Using  a  configuration  where  a  cluster  of  atoms  are 
trapped  using  a  FORT  beam  it  should  be  possible  to  demon¬ 
strate  the  ensemble-based  LSIIB  process,  without  using  a 
cavity.  Note  that,  for  an  ensemble,  this  process  occurs  in  a 
pure  A  system.  This  can  be  realized  by  using  the  states  |l), 
|2),  and  |3)  of  Fig.  8(b).  The  atoms  will  first  be  optically 
pumped  into  state  |l).  The  two-opposite  circularly  polarized 
beams,  detuned  from  the  F'  =  1  level,  will  then  be  used  to 
excite  the  Raman  transition  between  |l),  |2),  and  |3).  Since 
the  off-resonant  excitation  produces  negligible  population  in 
|2),  no  repump  beam  would  be  necessary. 

The  evidence  of  LSIIB  in  this  case,  of  course,  is  not  as 
simple  as  in  the  case  of  single  atoms  [6],  What  is  necessary 
in  this  case  is  to  demonstrate  that  one  and  only  one  atom  (on 
average)  can  be  excited  to  level  |3).  In  order  to  test  this,  one 
can  operate  the  Raman  transition  process  in  a  pulsed  mode, 
with  a  variable  pulse  duration.  After  applying  a  Raman  pulse, 
the  beams  exciting  the  1-2-3  transition  would  be  turned  off, 
and  another  Raman  transition  would  be  applied  to  transfer 
the  atom(s)  in  state  |l)  to  the  F= 2,  mF= 2  level.  A  cycling 
transition  that  couples  this  state  to  the  F=  3,  mF=  3  level  in 
the  D2  manifold  would  then  be  used  to  collect  fluorescence 
produced  from  this  atom.  The  fluorescence  could  be  col¬ 
lected  by  a  high  numerical  aperture  imaging  system,  and 
then  transported  through  a  fiber  to  a  single-photon  counter. 
By  taking  into  account  all  the  experimental  parameters,  one 
should  thus  be  able  to  determine  the  number  of  atoms  that 
were  excited  to  state  |3).  The  proof  of  ensemble-LSIIB 
would  be  established  by  showing  that  this  never  exceeds  1 
for  any  pump-pulse  duration. 

The  specific  geometry  of  the  experimental  setup  that 
could  be  used  to  demonstrate  the  ensemble-LSIIB  based 
controlled-NOT  gate  is  illustrated  schematically  in  Fig.  9(a). 


The  FORT  beams  shown  could  be  produced  from  a  single 
beam  with  imaging  optics.  Using  dichroic  mirrors,  additional 
beams  (to  be  used  for  producing  the  control  pulses  for  each 
ensemble)  could  be  brought  in,  parallel  to  the  FORT  beams. 
Additional  control  beams  would  be  brought  in  a  direction 
perpendicular  to  both  the  cavity  axis  and  the  FORT  axis.  This 
degree  of  freedom  would  enable  us  to  satisfy  the  polarization 
selection  rules. 

In  order  to  keep  the  process  as  clean  as  possible,  it  is 
advisable  to  consider  first  the  simplest  nontrivial  version  of 
the  controlled-NOT  gate.  Recall  the  most  general  form  of  the 
input  (unentangled)  state  of  the  two  ensembles,  and  the  cor¬ 
responding  entangled  state  after  the  controlled-NOT  opera¬ 
tion,  one  could  choose  the  parameters  such  that  a=/3 
=  1  /  J2,  £=  1 .  and  77=0.  The  entangled  state  produced  by  the 
controlled-NOT  gate  would  therefore  be  given  by  |  «/f) 7  /7 
=  [|^)/l  C\)  //  +  lQ)/|£>i  )//]/  \  2.  Here,  Ens-I  is  the  ensemble 
trapped  by  FORT1,  and  Ens-II  is  the  one  trapped  by  FORT2, 
for  example. 

In  order  to  demonstrate  that  such  a  state  has  indeed  been 
produced,  one  can  proceed  in  several  different  ways.  One 
option  is  to  apply  the  pulse  described  in  Fig.  7,  with  the 
pump  beam  parallel  to  the  direction  of  FORT1  (which  holds 
Ens-I)._If  Ens-I  is  in  state  |A),  then  application  of  a  pulse  so 
that  ( y  All , cgF/ 2 S)  tr=tt  will  produce  a  single  photon  in  the 
free-space  mode  (characterized  by  the  single-photon  Rabi 
frequency  of  gF).  As  described  in  detail  previously,  the  field 
at  nlc  is  coupled  to  the  |a)-|g)  transition  with  a  large  posi¬ 
tive  detuning,  and  the  field  for  gF  is  coupled  to  the  |c)-|g) 
transition  with  the  same  detuning.  The  frequency  of  the  free- 
space  mode  is  selected  automatically  in  order  to  satisfy  the 
two-photon  transition  condition.  If  Ens-I  is  in  state  IQ),  then 
application  of  the  same  pulse  will  not  produce  a  photon  in 
this  free-space  mode.  The  photon  in  this  mode  could  be  de¬ 
tected  by  using  an  additional  fiber-coupled  single-photon  de¬ 
tector,  preceded  by  polarization  and  spectral  filters. 

In  parallel  to  this  process,  a  pump  beam  (with  Rabi- 
frequency  fl'2C)  would  be  applied  along  FORT2  (which 


022323-7 


SHAHRIAR,  PATI,  AND  SALIT 


PHYSICAL  REVIEW  A  75,  022323  (2007) 


FIG.  9.  (a)  Schematic  illustra¬ 
tion  of  the  configuration  that  can 
be  employed  for  testing  the 
controlled-NOT  operation  using 
the  ensemble  LSIIB.  (b)  Illustra¬ 
tion  of  the  excitations  that  could 
be  employed  for  confirming  the 
controlled-NOT  gate  operation.  See 
text  for  details. 


houses  Ens-II),  coupled  to  the  |c)-|/t)  transition,  with  a  large, 
positive  detuning.  The  corresponding  free-space  mode  (with 
vacuum  Rabi  frequency  of  g' F),  coupled  to  the  \d)-\h)  tran¬ 
sition,  would  be  automatically  chosen  to  have  the  same  de¬ 
tuning,  in  order  to  satisfy  the  two-photon  resonance  condi¬ 
tion.  This  photon  would  be  passed  through  the  same  set  of 
three  filters,  and  then  detected  by  a  fiber-coupled  photon 
counter.  As  before,  application  of  a  pi  pulse  such  that 


(O' 2Cg ' Fl 28)  t  r  =  tt  would  produce  a  single  photon  in  the 
g' F  only  if  Ens-II  is  in  state  C , ) ;  otherwise,  there  would  be 
no  photons  in  this  mode.  Note,  again,  that  there  is  no  yTV 
factor  involved  in  this  expression.  A  coincidence  between  the 
photon  counts  found  in  the  gF  and  the  g' F  modes  will  prove 
that  the  controlled-NOT  gate  was  applied  correctly,  and  con¬ 
firm  the  entanglement  between  the  ensembles.  The  transi¬ 
tions  are  illustrated  schematically  in  Fig.  9(b). 
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B.4.  Ultra-Low  Power  Non-Linear  Optics  in  a  Tapered  Nanofiber  as  an  Enabling 
Technlogy  for  Optically  Linking  Ensemble-Based  Quantum  Bits 

Summary: 

As  described  in  section  B.3  above,  we  have  developed  a  new  model  of  optical 
quantum  computing  using  atomic  ensembles  as  quantum  bits.  In  our  original  model, 
these  would  make  use  of  moderate-Q  cavities  in  order  to  couple  nearest  neighbor 
quantum  bits  for  CNOT  operations.  Alternatively,  the  coupling  can  be  mediated  by  the 
evanescent  field  of  a  tapered  nano-fiber  (TNF).  A  key  advantage  of  using  the  TNF  for 
this  purpose  is  that  it  can  be  effortlessly  connected  to  regular  optical  fibers,  thereby 
linking  multiple  quantum  computers  in  the  form  of  a  quantum  internet,  or  for  Type-II 
quantum  computing. 

Under  this  grant,  we  have  studied  the  extremely  strong  interaction  between  the 
evanescent  field  of  a  TNF  and  an  atomic  vapor.  Specifically,  we  have  observed  ultra-low 
light  level  non-linear  optical  interactions  in  a  TNF  embedded  in  a  hot  rubidium  vapor  as 
well  as  in  a  cloud  of  trapped  atoms.  The  small  optical  mode  area  plays  a  significant  role 
in  the  optical  properties  of  the  hot  vapor  Rb-TNF  system,  allowing  nonlinear  optical 
interactions  with  nW  level  powers  in  a  hot  vapor  even  in  the  presence  of  transit-time 
dephasing  rates  much  larger  than  the  intrinsic  linewidth.  We  demonstrate  nonlinear 
absorption  and  V-type  electromagnetically-induced  transparency  with  cw  powers  below 
10  nW,  comparable  to  the  best  results  in  any  Rb-optical  waveguide  system.  Using  a  cloud 
of  trapped  atoms,  we  have  also  been  able  to  see  such  an  interaction  at  100  pW.  The  good 
performance  and  flexibility  of  the  Rb-TNF  system  makes  it  a  very  promising  candidate 
for  ultra-low  power  resonant  nonlinear  optical  applications  in  general,  and  ensemble 
quantum  computing  in  particular. 
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We  report  the  observation  of  low-light  level  optical  interactions  in  a  tapered  optical  nanofiber 
(TNF)  embedded  in  a  hot  rubidium  vapor.  The  small  optical  mode  area  plays  a  significant  role  in 
the  optical  properties  of  the  hot  vapor  Rb-TNF  system,  allowing  nonlinear  optical  interactions  with 
nW  level  powers  even  in  the  presence  of  transit-time  dephasing  rates  much  larger  than  the  intrinsic 
linewidth.  We  demonstrate  nonlinear  absorption  and  V-type  electromagnetically-induced  trans¬ 
parency  with  cw  powers  below  10  nW,  comparable  to  the  best  results  in  any  Rb-optical  waveguide 
system.  The  good  performance  and  flexibility  of  the  Rb-TNF  system  makes  it  a  very  promising 
candidate  for  ultra-low  power  resonant  nonlinear  optical  applications. 

PACS  numbers:  42. 50. Gy,  42.81.Qb 


The  interaction  of  light  with  atomic  ensembles  [1-3]  has 
created  a  number  of  advancements  in  generating  and  con¬ 
trolling  classical  and  quantum  states  of  light  [4-6].  In  or¬ 
der  to  enable  optical  interactions  with  few  photons  (as  de¬ 
sired  in  some  quantum  information  processing  studies[7]), 
much  effort  has  been  invested  into  increasing  the  interac¬ 
tion  between  atoms  and  photons.  This  is  often  done  by 
employing  high  finesse  optical  microcavities,  where  quan¬ 
tum  effects  are  readily  observed  between  a  single  atom 
and  a  single  photon  [8-11].  However,  these  systems  are 
extremely  complicated.  For  applications  where  single¬ 
photon  interactions  are  not  needed,  atomic  ensembles 
in  a  simple  nonresonant  geometry  are  used  to  observe 
quantum  optic  effects.  Previously,  most  work  was  per¬ 
formed  with  optical  beams  focused  into  an  alkali  vapor, 
where  the  use  of  a  bulk,  free-space  geometry  fundamen¬ 
tally  limited  the  interaction  strength  and/or  length  [6]. 
Recent  experiments  have  resulted  in  significant  increases 
in  both  interaction  strength  and  interaction  length  by  en¬ 
casing  an  atomic  vapor  into  a  waveguide[12,  13],  with  V- 
system  electromagnetically-induced  transparency  (EIT) 
demonstrated  in  a  photonic  crystal  fiber  with  approxi¬ 
mately  3  nW  of  pulsed  power[12].  In  both  of  these  cases, 
the  enhanced  interaction  comes  from  reducing  the  optical 
mode  area.  Here  we  explore  the  use  of  tapered  optical 
nanofibers  (TNF)  with  optical  mode  areas  smaller  than 
in  previous  experiments,  and  demonstrate  low  light  level 
nonlinear  interactions  with  a  rubidium  atomic  vapor. 

A  TNF  consists  of  a  single-mode  silica  optical  fiber 
where  the  original  fiber  diameter  is  reduced  to  the  sub¬ 
micron  range.  In  this  system,  the  optical  mode  which 
was  initially  guided  by  the  core/cladding  interface  is 
adiabatically  converted  with  negligible  loss  [14]  into  a 
mode  guided  by  the  cladding/environment  interface.  The 
evanescent  field  extends  into  the  surrounding  environ¬ 
ment,  which  for  the  purposes  of  this  work  consists  of  a  ru¬ 
bidium  atom  vapor.  The  use  of  an  evanescent  interaction 
maintains  the  benefit  of  independently  controlling  the 


interaction  strength  (through  cross-sectional  geometry) 
and  length  (due  to  the  invariance  of  the  evanescent  field 
along  the  waveguide),  and  allows  an  increase  of  interac¬ 
tion  strength  over  previous  waveguide  approaches[12,  13]. 
This  system  has  been  previously  used  to  study  evanes¬ 
cent  interactions  with  Cs  atoms,  such  as  probing  atom 
fluorescence[15,  16].  In  principle  the  interaction  length 
can  be  arbitrary  long,  just  like  in  any  waveguide,  with  a 
limit  due  to  intrinsic  optical  waveguide  loss.  For  a  TNF 
of  diameter  400  nm,  approximately  10%  of  the  optical 
field  interacts  with  the  rubidium  vapor.  For  this  diam¬ 
eter,  we  estimate  that  the  total  number  of  atoms  inter¬ 
acting  with  the  TNF  as  1400/cm,  which  for  the  typical 
taper  length  of  3  mm  in  this  work  gives  a  total  number  of 
atoms  of  ss  450  (assuming  a  Rb  atomic  density  of  6  x  1012 
atoms/cm3  for  a  vapor  at  100°C). 

We  have  investigated  numerically  the  performance  of 
the  Rb-TNF  system  for  different  TNF  geometries.  Fig¬ 
ure  1  shows  a  plot  of  the  fraction  of  the  optical  energy 
located  in  the  evanescent  field  and  the  optical  mode  area 
for  a  TNF  versus  diameter  for  an  optical  wavelength  of 
780  nm,  corresponding  to  the  Rb  D2  manifold.  We  see 
that  as  the  diameter  decreases  the  fraction  of  the  optical 
energy  that  extends  into  the  external  Rb  vapor  increases 
monotonically.  The  optical  mode  area  decreases  as  the 
diameter  decreases  to  about  0.45  /iin  due  to  the  added 
dielectric  confinement,  corresponding  to  a  minimum  op¬ 
tical  mode  area  of  0.18  (im2.  As  the  diameter  is  reduced 
further,  the  mode  quickly  delocalizes,  causing  a  rise  in 
mode  area.  This  minimum  mode  area  is  significantly 
smaller  than  that  found  in  alternate  waveguiding  geome¬ 
tries  (>  10  x  smaller  than  ARROW  waveguides  [13]  and 
>  lOOx  smaller  than  PCF  waveguides[12]),  and  is  com¬ 
parable  to  the  atomic  scattering  cross-section  of  3A2/27t. 
Therefore,  we  expect  that  the  TNF  should  have  a  dra¬ 
matically  increased  interaction  with  an  atomic  vapor. 

However,  in  an  open  environment  such  as  the  Rb-TNF 
system,  this  mode  area  is  so  miniscule  that  a  single  atom 
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FIG.  1:  Calculated  mode  area  and  external  energy  fraction  for 
a  TNF  (A=780  nm).  Inset  shows  the  electric  field  distribution 
for  the  fundamental  mode  of  a  400  nm  diameter  TNF. 


(with  a  velocity  typically  found  in  a  hot  vapor  cell)  passes 
through  the  optical  field  of  the  TNF  in  a  few  nanosec¬ 
onds.  This  time  is  much  shorter  than  the  intrinsic  decay 
time  of  the  atom  («  25  ns),  leading  to  an  effect  known 
as  transit-time  broadening.  Fundamentally,  transit-time 
broadening  can  be  understood  by  considering  the  opti¬ 
cal  interaction  in  the  atom’s  reference  frame.  Here,  the 
stationary  atom  experiences  a  temporally-varying  opti¬ 
cal  field,  with  a  temporal  shape  dependent  on  the  trans¬ 
verse  optical  mode  field.  In  this  case,  the  optical  absorp¬ 
tion  linewidth  is  given  by  a  combination  of  the  intrinsic 
atomic  decay  rate  and  the  transit  time  decay  rate  (de¬ 
termined  from  the  Fourier  transform  of  the  time- varying 
optical  field),  leading  to  a  broader  absorption  line  when 
transit-time  effects  dominate.  Additionally,  the  presence 
of  transit-time  broadening  can  significantly  increase  the 
threshold  to  observe  nonlinear  interactions.  Figure  2 
shows  a  calculation  of  the  transit  time  dephasing  rate 
expected  for  a  Rb  vapor  interacting  with  a  TNF  of  var¬ 
ious  diameters.  For  this  calculation  we  used  the  exact 
TNF  optical  field  profile  to  determine  the  dephasing  rate 
for  an  atom  passing  through  the  evanescent  field,  as¬ 
suming  an  isotropic  angular  distribution  and  that  the 
atoms  which  struck  the  TNF  bounced  off  with  no  de¬ 
lay  and  without  dephasing  (this  gives  an  estimate  of  the 
smallest  transit  time  linewidth  possible) .  This  procedure 
was  performed  for  all  possible  atomic  speeds,  with  the  fi¬ 
nal  linewidth  determined  by  integrating  over  a  Maxwell- 
Boltzmann  velocity  distribution  for  a  vapor  with  a  tem¬ 
perature  of  100°C.  The  calculation  shows  that  dephasing 
rates  are  on  the  order  of  100  MHz,  which  is  much  larger 
than  the  intrinsic  atom  linewidth  («6  MHz),  with  a  min¬ 
imum  value  of  89  MHz  for  a  TNF  diameter  of  400  nm. 

An  easy  way  to  describe  the  suitability  of  a  system  for 
nonlinear  optical  effects  is  to  consider  the  power  needed 
to  saturate  an  atomic  vapor.  This  is  given  by  the  relation 
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FIG.  2:  Calculated  transit  time  linewidth  and  saturation 
power  for  a  hot  rubidium  vapor  (vapor  temperature  of  100°  C) 
coupled  to  a  TNF.  The  saturation  calculation  assumes  an  in¬ 
trinsic  saturation  intensity  of  1.5  mW/cm2  and  a  Doppler 
width  of  575  MHz  for  the  bulk  vapor. 


Psat  =  4a*Am(7T/7o)2,  where  Isat  is  the  saturation  in¬ 
tensity  for  an  atomic  transition,  Am  is  the  optical  mode 
area,  and  7t/7o  represents  the  ratio  of  the  total  linewidth 
(including  transit  time  dephasing  and  Doppler  broaden¬ 
ing)  and  the  intrinsic  linewidth.  Figure  2  shows  the  cal¬ 
culated  nonlinear  saturation  power  as  a  function  of  di¬ 
ameter  corresponding  to  an  isotropic  Rb  vapor  of  100°C, 
with  Isat  =  1.5mW/cm2  for  the  D2  manifold.  In  this  case 
the  Doppler  width  of  575  MHz  is  much  larger  than  transit 
time  dephasing,  so  that  the  calculated  saturation  power 
is  effectively  independent  of  transit  time.  The  data  show 
that  due  to  the  small  mode  area  the  nonlinear  saturation 
power  is  on  the  order  of  tens  of  nanowatts  even  in  the 
presence  of  significant  Doppler  broadening. 

In  our  experimental  setup  (Fig.  3(a)),  a  TNF  was  in¬ 
serted  into  a  chamber  containing  hot  rubidium  vapor, 
which  contained  a  natural  mixture  of  Rb85  and  Rb8 '  iso¬ 
topes.  A  Newport  single  mode  fiber  at  780  nm  was  drawn 
adiabatically  into  a  TNF  with  a  waist  diameter  of  ap¬ 
proximately  400  nm  using  the  ’flame-brush’  technique[17] 
with  a  hydrogen-air  torch.  The  resulting  taper  profile  is 
nearly  exponential  with  a  1/e  length  of  ~  3  mm.  During 
the  adiabatic  tapering  process,  monitoring  of  fiber  trans¬ 
mission  showed  negligible  loss  (ranging  from  1  —  10%) 
for  final  taper  waist  diameters  of  approximately  400  nm. 
The  TNF  is  mounted  to  a  copper  chuck  using  UV-curable 
epoxy  and  inserted  into  the  vacuum  chamber  using  a 
set  of  Teflon  fiber  feedthroughs  [18],  which  maintain  fiber 
continuity  into  and  out  of  the  vacuum  chamber.  A  gate 
valve  isolates  the  optical  fiber  during  loading/unloading 
from  a  rubidium  metal  source.  During  experiments,  the 
vacuum  chamber  was  heated  to  ss  100°C  (to  minimize 
Rb  condensation) ,  and  Rb  vapor  was  created  by  heating 
the  source  region  to  ss  200°C,  with  the  chamber  pressure 
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FIG.  3:  (a)  Diagram  of  the  experimental  setup.  The  tapered 
nanofiber  passes  through  the  Rb  vapor-containing  chamber 
using  teflon  fiber  feedthroughs.  A  polarization  beamsplitter 
is  used  to  inject  a  strong  counterpropagating  cross-polarized 
pump  for  the  saturated  absorption  and  EIT  measurements, 
(b)  Transmitted  spectrum  for  a  TNF  and  reference  cell  for 
the  D2  manifold  of  Rb  vapor  at  780  nm. 

maintained  at  5  mTorr  by  a  roughing  pump.  For  these 
measurements,  we  made  no  attempt  to  prevent  dephasing 
for  atoms  which  struck  the  TNF.  Assuming  the  vapor  in¬ 
side  the  chamber  is  at  an  average  temperature  of  100°C, 
the  atomic  density  is  estimated  to  be  6  x  1012/cm3. 

For  the  first  set  of  experiments,  a  weak  probe  beam 
(approximately  10  nW)  obtained  from  a  Ti:Sapphire 
laser  was  transmitted  through  the  TNF,  and  the  optical 
transmission  was  monitored  with  an  avalanche  photodi¬ 
ode.  The  frequency  of  the  probe  was  scanned  over  the 
Doppler  broadened  spectrum  of  the  D2  manifold.  Fig¬ 
ure  3(b)  shows  the  transmission  spectrum  through  the 
TNF  (lower  spectrum),  and  a  reference  Rb  vapor  cell 
kept  at  100°C  (upper  spectrum).  The  dip  on  the  left 
corresponds  to  the  F=2  to  F’=2  transition  in  the  D2 
manifold  of  Rb87 ,  while  the  dip  on  the  right  corresponds 
to  the  F=3  to  F’=3  transition  in  the  same  manifold  of 
Rb85.  Far  away  from  resonance,  the  probe  was  found 
to  be  attenuated  by  about  20%  due  to  a  combination 
of  losses  at  the  input  couplers  and  by  TNF  absorption. 
The  data  clearly  show  transmission  dips  representative 
of  rubidium  vapor  absorption,  with  peak  locations  sim¬ 
ilar  to  that  of  the  reference  cell.  However,  the  shape  of 
the  transmission  dips  for  the  Rb-TNF  system  is  slightly 
different  than  that  for  the  reference  cell,  which  is  due  to 
both  differences  in  Doppler  broadening  and  transit  time 
dephasing  present  in  this  experimental  measurement. 

In  optical  saturation,  the  nonlinear  absorption  coeffi¬ 
cient  ajvL  =  ao/(l  +  P/Psat),  where  «o  is  the  linear  ab¬ 
sorption  coefficient  and  P  is  the  incident  optical  power. 
As  the  incident  optical  power  increases  in  a  nonlinear 
medium,  the  relative  absorption  decreases,  and  therefore 
the  relative  transmission  increases.  Figure  4  shows  the 


FIG.  4:  Transmission  spectrum  for  a  TNF  for  increasing 
probe  powers  of  {2,  4,  8, 10}nW  for  the  Rb  D2  transition.  The 
spectra  show  optical  saturation  with  a  saturation  power  of  «  8 
nW. 


transmitted  power  spectrum  measured  for  a  variety  of  in¬ 
cident  power  levels  for  a  TNF  with  a  waist  diameter  of  ap¬ 
proximately  400  nm.  On  the  left  edge,  the  output  power 
is  close  to  its  value  observed  far  away  from  resonance. 
The  bottom  (green)  trace  corresponds  to  an  off-resonant, 
output  probe  power  of  about  2  nW.  At  this  power  level, 
the  probe  is  found  to  be  absorbed  very  strongly.  As  the 
probe  power  is  increased,  the  relative  transmission  in¬ 
creases,  as  expected  for  a  nonlinear  medium.  The  ab¬ 
sorption  decrease  with  incident  power  can  be  used  to  ex¬ 
tract  the  saturation  power,  which  is  approximately  8  nW. 
This  value  is  lower  than  the  theoretical  value  which  we 
attribute  to  a  nonisotropic  atomic  vapor  in  our  experi¬ 
ment.  Assuming  an  effective  Doppler  width  of  210  MHz 
(inferred  from  measurements  described  below) ,  we  obtain 
a  predicted  saturation  power  of  4  nW.  Accounting  for  the 
TNF  profile,  the  saturation  power  increases  to  «  8  nW, 
in  good  agreement  with  the  measurement. 

Next,  a  pump-probe  measurement  was  performed, 
with  the  strong  counterpropagating  pump  beam  cross- 
polarized  with  respect  to  the  probe.  This  measurement 
was  performed  on  the  Di  manifold  in  order  to  delin¬ 
eate  more  clearly  the  spectral  features  in  the  presence 
of  large  transit  time  broadening.  Figure  5  illustrates  the 
probe  transmission  spectrum  (power  ss  1  nW)  for  coun¬ 
terpropagating  pump  powers  of  10  nW  (middle,  green) 
and  30  nW  (lower,  red),  for  a  TNF  with  a  waist  diam¬ 
eter  of  approximately  400  nm.  For  comparison,  a  ref¬ 
erence  trace  is  also  shown  for  a  conventional  vapor  cell 
(blue).  The  peaks  on  the  left  and  the  right  correspond 
to  self-induced  transparency  (SIT)  for  Rb85  atoms  with 
zero  axial  velocities,  corresponding  to  the  F=3-F’=2  and 
the  F=3-F’=3  hyperfine  transitions,  respectively.  This 
measurement  is  similar  to  that  of  conventional  saturated 
absorption  spectroscopy,  with  a  V-type  three  level  EIT 
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FIG.  5:  Saturated  absorption  spectrum  with  cross-polarized 
pump  and  probe  beams  for  a  Rb  vapor  cell  (blue)  and  a  Rb- 
TNF  system  (green  and  red)  for  the  Di  manifold  of  Rb85.  The 
two  side  peaks  correspond  to  the  Doppler- free  F=3-F‘=2  and 
F=3-F‘=3  hyperhne  lines,  which  are  composed  of  atoms  with 
near  zero  longitudinal  velocity.  These  widths  indicate  the 
power-broadened  hyperfine  linewidth,  including  transit  time 
dephasing.  The  center  peak  corresponds  to  V-system  EIT. 

peak  (with  orthogonally  polarized  pump  and  probe)  lo¬ 
cated  at  the  cross-over  resonance  location.  This  peak  is 
comprised  of  V-type  EIT  measurements  of  two  groups 
of  atoms  (illustrated  in  figure  5):  (a)  one  with  an  axial 
velocity  such  that  the  probe  excites  the  F=3-F’=3  transi¬ 
tion  and  the  pump  excites  the  F=3-F’=2  transition,  and 
(b)  the  other  with  an  equal  but  opposite  axial  velocity 
so  that  the  probe  excites  the  F=2-F’=3  transition  and 
the  pump  excites  the  F=3-F’=2  transition.  The  EIT  sig¬ 
nal  for  the  TNF  is  readily  apparent  for  the  lower  pump 
power  measurement,  where  the  peak  is  clearly  separated 
from  the  SIT  peaks.  The  amplitude  of  this  peak  is  mea¬ 
surably  greater  than  the  background  signal  for  both  the 
low  and  high  power  spectra,  even  though  the  high  power 
spectrum  is  less  resolvable  due  to  power  broadening. 

For  the  conventional  vapor  cell,  the  transparency  seen 
at  the  EIT  resonance  is  larger  than  that  at  the  SIT  peaks 
due  to  the  larger  number  of  atoms  in  the  combined  ±180 
MHz  Doppler-shifted  groups  (there  is  approximately  1.5 
times  the  number  of  atoms  in  an  isotropic  vapor  cell, 
which  arises  from  the  575  MHz  width  of  the  Doppler  dis¬ 
tribution)  than  the  individual  zero  axial  velocity  groups. 
Furthermore,  the  peak  widths  are  the  power  broadened 
individual  hyperfine  transitions  («  6  MHz  power  broad¬ 
ened  to  «  18  MHz).  However,  for  the  TNF  signal,  there 
is  a  noticeable  difference  between  the  linewidths  and  rel¬ 
ative  amplitudes  of  the  SIT  and  EIT  peaks.  Here,  the 
line  broadening  for  the  lower  (middle)  power  curve  is 
due  primarily  to  the  transit  time  effect.  The  transit 
time  broadening  is  about  110  MHz,  close  to  the  theo¬ 
retical  prediction.  The  dramatically  smaller  amplitude 
of  the  V-system  EIT  peak  in  the  TNF  measurement  is 


likely  due  to  a  nonisotropic  atomic  velocity  distribution 
in  the  Rb-TNF  cell,  which  leads  to  a  smaller  population 
of  atoms  with  the  correct  longitudinal  Doppler-shifted 
velocity  component  for  the  EIT  interaction.  The  rel¬ 
ative  magnitude  of  the  three  absorption  peaks  is  con¬ 
sistent  with  a  210  MHz  wide  effective  Doppler  distribu¬ 
tion,  which  also  gives  good  agreement  with  the  satura¬ 
tion  power  data  in  figure  4.  The  effective  reduction  of 
the  Doppler  width  may  be  attributable  to  the  geometry 
of  the  cell  and  the  relative  position  of  the  TNF  and  the 
vacuum  pump.  Further  investigation  will  be  carried  out 
in  the  future  to  study  this  effect. 

In  conclusion,  we  have  shown  that  a  TNF  coupled  to  a 
Rb  vapor  exhibits  optical  saturation  and  EIT  at  very  low 
pump  power  levels  below  10  nW,  comparable  to  the  best 
results  in  any  Rb-waveguide  system.  Numerical  model¬ 
ing  indicates  that  the  TNF’s  small  mode  area  results  in 
a  significant  linewidth  broadening  due  to  the  short  tran¬ 
sit  time  of  hot  Rb  atoms  through  the  evanescent  field, 
which  combined  with  Doppler  broadening  increases  the 
threshold  for  nonlinear  optical  effects.  However,  this  is  a 
limitation  of  the  large  atomic  velocities  in  a  high  temper¬ 
ature  Rb  vapor,  which  can  be  reduced  in  principle.  In 
particular,  the  use  of  cooled  atoms  from  a  magnetoop¬ 
tic  trap[15,  16]  will  allow  the  elimination  of  transit  time 
and  Doppler  broadening,  which  should  lead  to  nonlinear 
interactions  with  pW  level  optical  powers. 
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B.5.  Identification  and  Elimination  of  Driver  Phase  Correlated  Fluctuations  in  the 
Rotation  of  a  Strongly  Driven  Quantum  Bit 

Summary: 

The  need  to  maximize  the  number  of  operations  of  a  quantum  bit  within  it’s 
decoherence  time  may  require  the  ratio  of  Rabi  frequency  to  transition  frequency  to  be 
large  enough  to  invalidate  the  rotating-wave  approximation.  This  is  particularly  important 
for  the  quantum  computing  model  in  NV-diamond,  and  is  also  relevant  for  the  atomic 
ensemble  based  quantum  computing  model  summarized  in  section  B.3.  We  have  shown, 
under  work  perfonned  for  this  grant,  that  the  state  of  the  quantum  bit  under  any  initial 
condition  in  such  a  case  depends  explicitly  on  the  phase  of  the  driving  field,  resulting  in 
driver-phase-correlated  fluctuations  and  a  violation  of  the  rule  that  the  degree  of 
excitation  depends  only  on  the  pulse  area.  This  is  due  to  the  interference  of  the 
excitations  caused  by  the  co-rotating  and  counter-rotating  fields,  and  is  a  significant 
source  of  error,  correctable  only  by  controlling  the  driver  phase. 

This  effect  also  shows  how  one  can  measure  the  amplitude  and  absolute  (i.e., 
temporal  and  initial)  phase  of  a  monochromatic  microwave  field  at  a  specific  point  of 
space  and  time.  In  addition  to  suppressing  error  in  high-speed  quantum  computing,  this 
measurement  process  has  other  potential  applications,  including  wavelength  quantum 
teleportation.  We  have  shown  how  such  a  measurement  can  indeed  be  made  using 
resonant  atomic  probes  via  detection  of  incoherent  fluorescence  induced  by  a  laser  beam. 
This  measurement  is  possible  due  to  self-interference  effects  between  the  positive-  and 
negative-frequency  components  of  the  field.  In  effect,  the  small  cluster  of  atoms  here  act 
as  a  highly  localized  pickup  coil,  and  the  fluorescence  channel  acts  as  a  transmission  line. 


Reproduction  of  published  journal  papers:  (4+4  pages) 

This  work  has  been  published  in  two  articles  in  Physical  Review  A.  These  are 
reproduced  below. 
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49 


PHYSICAL  REVIEW  A  69,  032308  (2004) 


Driver-phase-correlated  fluctuations  in  the  rotation  of  a  strongly  driven  quantum  bit 


19  |  9  9 

M.  S.  Shahriar,  ’  Prabhakar  Pradhan,  ’  and  Jacob  Morzinski 

1  Department  of  Electrical  and  Computer  Engineering,  Northwestern  University,  Evanston,  Illinois  60208,  USA 
2 Research  Laboratory  of  Electronics,  Massachusetts  Institute  of  Technology,  Cambridge,  Massachusetts  02139,  USA 

(Received  20  July  2002;  published  16  March  2004) 


The  need  to  maximize  the  number  of  operations  of  a  quantum  bit  within  its  decoherence  time  may  require 
the  ratio  of  Rabi  frequency  to  transition  frequency  to  be  large  enough  to  invalidate  the  rotating-wave  approxi¬ 
mation.  The  state  of  the  quantum  bit  under  any  initial  condition  then  depends  explicitly  on  the  phase  of  the 
driving  held,  resulting  in  driver-phase-correlated  fluctuations  and  a  violation  of  the  rule  that  the  degree  of 
excitation  depends  only  on  the  pulse  area.  This  is  due  to  the  interference  of  the  excitations  caused  by  the 
corotating  and  counterrotating  fields,  and  is  a  significant  source  of  error,  correctable  only  by  controlling  the 
driver  phase.  We  present  a  scheme  for  observing  this  effect  under  currently  realizable  parameters. 
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In  order  to  minimize  the  decoherence  rate  of  a  two-state 
quantum  bit  (qubit)  embodied  in  a  massive  particle,  one  of¬ 
ten  chooses  to  use  low-energy  transitions.  In  general,  one  is 
interested  in  performing  these  transitions  as  fast  as  possible 
[1-5]  which  demands  a  strong  Rabi  frequency.  The  ratio  of 
Rabi  frequency  to  qubit  transition  frequency  is  therefore  not 
necessarily  very  small,  thus  invalidating  the  so-called 
rotating-wave  approximation  (RWA).  A  key  effect  due  to 
violation  of  the  RWA  (VRWA)  is  the  so-called  Bloch-Siegert 
shift  [6-9]  which  is  negligible  in  optical  transitions,  but  is 
manifested  in  nuclear  magnetic  resonance  [10],  Here,  we 
show  that  VRWA  leads  to  another  important  effect,  which 
can  lead  to  controllable  errors  that  are  significant  on  the  scale 
of  precisions  envisioned  for  a  functioning  quantum  computer 
[11].  Specifically,  we  show  that  under  VRWA  the  population 
difference  between  the  two  levels  of  the  quantum  bit,  with 
any  initial  condition,  depends  explicitly  on  the  phase  of  the 
driving  field  at  the  onset  of  an  excitation  pulse,  which  is  a 
violation  of  the  rule  [6]  that  for  a  two-level  system  starting  in 
the  ground  state,  the  population  difference  is  a  function  of 
the  integral  of  the  field  amplitude  over  the  pulse  duration  and 
does  not  depend  on  the  phase  of  the  field.  We  provide  a 
physical  interpretation  of  this  effect  in  terms  of  an  interfer¬ 
ence  of  the  excitations  caused  by  the  corotating  and  counter¬ 
rotating  fields,  and  present  a  scheme  for  observing  this  effect 
under  currently  realizable  parameters. 

To  see  the  implication  of  this  result,  consider  a  scenario 
where  one  has  a  qubit,  initialized  to  the  ground  state,  and 
would  like  to  prepare  it  to  be  in  an  equal  superposition  of  the 
ground  and  excited  states.  To  this  end,  one  would  apply  a 
resonant  pulse  with  an  area  of  77/ 2  starting  at  a  time  f=  fo¬ 
under  the  RWA,  one  does  not  have  to  know  what  the  abso¬ 
lute  phase  of  the  field,  <f>P ,  is  at  f0,  and  the  population  differ¬ 
ence  for  the  qubit  would  be  zero.  Under  VRWA,  however, 
the  desired  excitation  would  only  occur  if  <pP=0.  Otherwise, 
the  population  difference  would  have  a  component  varying 
as  77  sin(2  (f>P),  where  77  is  a  parameter  that  is  proportional  to 
the  ratio  of  Rabi  frequency  to  transition  frequency.  Suppose 
one  has  to  apply  this  pulse  to  many  such  qubits,  with  a  po¬ 
tentially  different  cf>P  for  each  (e.g.,  because  the  pulses  are 
applied  at  different  times  or  the  qubits  are  spatially  sepa¬ 


rated),  but  with  identical  pulse  areas.  The  population  differ¬ 
ence  for  the  qubits  will  then  exhibit  a  fluctuation,  correlated 
to  their  respective  values  of  <bP.  For  a  quantum  computer, 
this  variation  would  represent  a  source  of  error.  For  some 
experiments  (e.g..  Ref.  [5]),  the  value  of  77  is  already  close  to 
0.01,  so  that  the  magnitude  of  this  error  is  much  larger  than 
the  ultimate  accuracy  (10-6)  desirable  for  a  large-scale  quan¬ 
tum  computer  [11]  and  must  be  controlled. 

To  illustrate  this  effect,  we  consider  an  ideal  two-level 
system  where  a  ground  state  |0)  is  coupled  to  a  higher-energy 
state  |l).  We  also  assume  that  the  0<->l  transition  is  mag¬ 
netic  dipolar,  with  a  transition  frequency  a>,  and  the  magnetic 
field  is  of  the  form  B=B0  cos(cor+  cf>).  We  now  summarize 
briefly  the  two-level  dynamics  without  the  RWA.  In  the  di¬ 
pole  approximation,  the  Hamiltonian  can  be  written  as 


H  =  e(cr0  -  <u)/2  +  g(t)(Tx,  (1) 

where  g(r)  =  -g0[exp(/cor+;^>)  +  c.c.]/2,  07  are  the  Pauli  ma¬ 
trices,  and  £=10  corresponds  to  resonant  excitation.  The 
state  vector  is  written  as 


\m= 


c0(t ) 

C1(t)_  • 


(2) 


We  perform  a  rotating-wave  transformation  by  operating  on 
|£(f))  with  the  unitary  operator  Q,  where 


Q  =  ((To  +  crz)/2  +  exp(/ cor  +  itjf) {a0  -  cr,)/2 .  (3) 

The  Schrodinger  equation  then  takes  the  form  (setting  h=  1) 
\f)=-iH{t)\f(t))  where  the  effective  Hamiltonian  is  given  by 

H=  a(t)a+  +  a(t)cr_,  (4) 

with  a(f)  =  -(g0/2)[exp(-/2cor-i2(A)  + 1],  and  in  the  rotating 
frame  the  state  vector  is 
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c0(i) 


(5) 


Now,  one  may  choose  to  make  the  RWA,  corresponding  to 
dropping  the  fast-oscillating  term  in  a{t).  This  corresponds 
to  ignoring  effects  (such  as  the  Bloch-Siegert  shift)  of  the 
order  of  (gQ/ co),  which  can  easily  be  observable  in  an  experi¬ 
ment  if  g0  is  large  [6-10].  On  the  other  hand,  by  choosing  g0 
to  be  small  enough,  one  can  make  the  RWA  for  any  value  of 
co.  We  explore  both  regimes  in  this  paper.  As  such,  we  find 
the  general  results  without  the  RWA. 

From  Eqs.  (4)  and  (5),  one  gets  two  coupled  differential 
equations 


C0d)  =  -  (go/2)[l  +  exp(-  i2a>t—  i2cj)j\Cx(t),  (6a) 


C,d)  =  -  (go/2)[l  +  exp(+  i2cot+  i2cj)]C0(t) .  (6b) 

We  assume  |C0d)|2=  1  to  be  the  initial  condition  and  proceed 
further  to  find  an  approximate  analytical  solution  of  Eq.  (6). 
Given  the  periodic  nature  of  the  effective  Hamiltonian,  the 
general  solution  to  Eq.  (6)  can  be  written  in  the  form 


oo 

\m=  2 

n=- oo 


exp[/z(-  i2a>t  -  i2<t>)\. 


(7) 


Inserting  Eq.  (7)  into  Eq.  (6)  and  equating  coefficients  with 
same  the  frequencies,  one  gets  for  all  n, 

a„  =  i2ncoa„  +  igo(bn  +  bn_x)/2,  (8a) 


bn  =  i2ncobn  +  igo(an  +  a„+x)/2.  (8b) 

Here,  the  coupling  between  a0  and  bQ  is  the  conventional  one 
present  when  the  RWA  is  made.  The  couplings  to  the  nearest 
neighbors,  a±1  and  b±x,  are  detuned  by  an  amount  2 co  and  so 
on.  To  the  lowest  order  in  (g0/ co),  we  can  ignore  terms  with 
\n  \  >  1,  thus  yielding  a  truncated  set  of  six  equations 


ao=ig0(b0  +  b_x)/2. 

(9a) 

b0  =  ig0(a0  +  ax)/2. 

(9b) 

dx=i2coax  +  igo(bx  +  b0)/2. 

(9c) 

bx  =  i2cobx  +  ig0ax/2. 

(9d) 

d_x  =  -  i2coa_x  +  ig0b_x/2. 

(9e) 

1  =  -  i2cob_x  +  /go(a-i  +  flo)/2. 

(9f) 

We  consider  g0  to  have  a  time  dependence  of  the  form 
god)=<?OM[l  -exp(-f/rJM,)],  where  the  switching  time  con¬ 
stant  TiwS>  o>-1  ,gQXM.  We  can  solve  these  equations  by  em¬ 
ploying  the  method  of  adiabatic  elimination,  which  is  valid 
to  first  order  in  rj=  (g0/ 4oi).  Note  that  77  is  also  a  function  of 
time  and  can  be  expressed  as  77(f)  =  ?70[1 -exp(-f/Tiw)], 


where  ?70=  {gmIAco).  We  consider  first  Eqs.  (9e)  and  (9f).  In 
order  to  simplify  these  two  equations  further,  one  needs  to 
diagonalize  the  interaction  between  a_x  and  .  Define  fi_ 
=  (a_x-b_x)  and  /r+  =  (a_x  +  b_x),  which  now  can  be  used  to 
reexpress  these  two  equations  in  a  symmetric  form  as 

=  -  i(2co  +  g0/2)fjc_  -  ig0a0/2,  (10a) 

A+  =  -  i(2co  -  g0/2)fj,+  +  ig0a0/2.  (10b) 

Adiabatic  following  then  yields  (again,  to  lowest  order  in  rj) 
/j,_~-rja0  and  /jl+~  7}a0,  which  in  turn  yields  cz^  —  O  and 
b_x  ~  77 a0.  In  the  same  manner,  we  can  solve  Eqs.  (9c)  and 
(9d),  yielding  ax~-  77/7,,  and  bx~  0. 

Note  that  the  amplitudes  of  a_x  and  bx  are  vanishing  (each 
proportional  to  rj1)  to  lowest  order  in  77  and  thereby  justify¬ 
ing  our  truncation  of  the  infinite  set  of  relations  in  Eq.  (9).  It 
is  easy  to  show  now 

d0  =  ig0b0/2  +  z'A(f)fl0/2,  (11a) 

bo  =  igoaQ/2-  iA(t)bo/2,  (lib) 

where  A(f)=go(f)/4<w  is  essentially  the  Bloch-Siegert  shift. 
Equation  (11)  can  be  thought  of  as  a  two-level  system  ex¬ 
cited  by  a  field  detuned  by  A.  For  simplicity,  we  assume  that 
this  detuning  is  dynamically  compensated  for  by  adjusting 
the  driving  frequency  &>.  This  assumption  does  not  affect  the 
essence  of  the  results  to  follow,  since  the  resulting  correction 
to  77  is  negligible.  With  the  initial  condition  of  all  the  popu¬ 
lation  in  |0)  at  t= 0,  the  only  nonvanishing  (to  lowest  order  in 
rj)  terms  in  the  solution  of  Eq.  (9)  are 

a0(t)  »  cos[go(f)f/2],  b0(t)  =  i  sin[^g(f)f/2], 

ax{t)  ~  -  ir]  sin[go(f)f/2],  b_x{t)  =  77  cos[go(f)f/2], 
where 


g{x(t)  =\lt\  g0(t')dt'  =  gQ[l  -  (t/rsw)  1  exp(-f/rJW)]. 
do 

We  have  verified  this  solution  via  numerical  integration  of 
Eq.  (6)  as  shown  later.  Inserting  this  solution  into  Eq.  (6)  and 
reversing  the  rotating-wave  transformation,  we  get  the  fol¬ 
lowing  expressions  for  the  components  of  Eq.  (2): 

C0(t)  =  cos[go(f)f/2]  -  2?7 X  sin[go(f)f/2],  (12a) 

Cx{t)  =  ie-l(“,+*\s\n[g'0{t)tl2\  +  277S*cos[^(f)f/2]}, 

(12b) 

where  we  have  defined  2  =  (z72)exp[-z(2&>r+2<^>)].  To  low¬ 
est  order  in  77,  this  solution  is  normalized  at  all  times.  Note 
that  if  one  wants  to  carry  this  excitation  on  an  ensemble  of 
atoms  using  a  tt/2  pulse  and  measure  the  population  of  the 
state  |l)  after  the  excitation  terminates  [at  1=  t  when 
g'(r)r/2=  7t/2],  the  result  would  be  an  output  signal  given 
by 
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FIG.  1.  Schematic  illustration  of  an  experimental  arrangement 
for  measuring  the  phase  dependence  of  the  population  of  the  excited 
state  |l):  (a)  The  microwave  field  couples  the  ground  state  (|0))  to 
the  excited  state  (|l)).  A  third  level,  1 2),  which  can  be  coupled  to  |l) 
optically,  is  used  to  measure  the  population  of  1 1)  via  fluorescence 
detection,  (b)  The  microwave  field  is  turned  on  adiabatically  with  a 
switching  time  constant  tsw,  and  the  fluorescence  is  monitored  after 
a  total  interaction  time  of  r. 

|Ci(#o(7'),0)|2=  ^  +2r/  sin(20r)]>  (13) 

where  we  have  defined  the  phase  of  the  field  at  /  =  r  to  be 
4>t=  oj7+  <f>.  This  signal  contains  information  of  both  the  am¬ 
plitude  and  phase  of  the  driving  field. 

This  result  can  be  appreciated  best  by  considering  an  ex¬ 
perimental  arrangement  of  the  type  illustrated  in  Fig.  1.  Con¬ 
sider,  for  example,  a  collection  of  S7Rb  atoms,  caught  in  a 
dipole  force  trap,  where  the  states  |0)  =  52S1/2:|F=l,m=l) 
and  |l)  =  52S1/2:  \F=2,m=2)  form  the  two-level  system. 
These  states  differ  in  frequencies  by  6.683  47  GHz.  When 
illuminated  by  resonant  right-circularly  polarized  light  at  a 
frequency  of  3.844  X  1014  Hz,  state  |l)  couples  only  to  the 
state  |2)  =  52F>3/2: |F=3,m  =  3),  which  in  turn  can  decay  only 
to  state  |l).  This  cycling  transition  can  thus  be  used  to  pump 
the  system  into  state  |l).  When  a  right-circularly  polarized 
microwave  field  at  6.683  47  GHz  is  applied,  state  |l)  couples 
only  to  state  |0),  even  when  the  RWA  approximation  breaks 
down.  The  strong-coupling  regime  (e.g.,  r/(l  of  the  order  of 
0.1)  can  be  reached,  for  example,  by  using  a  superconduct¬ 
ing,  high-Q  (1010)  microwave  cavity  [12].  The  theoretical 
model  developed  above  is  then  a  valid  description  of  the 
coupling  between  |0)  and  |l). 

The  strong  microwave  field  is  turned  on  adiabatically  with 
a  switching  time  constant  tsw,  starting  at  t= 0.  After  an  inter¬ 
action  time  of  t,  chosen  so  that  go(r)T=7r/2,  the  population 
of  state  |l)  can  be  determined  by  coupling  this  state  to  the 
state  \2)  with  a  short  (faster  than  1  / a>  and  1/goAf)  laser  pulse 
and  monitoring  the  resulting  fluorescence  [13]. 

We  have  simulated  this  process  explicitly  for  the  follow¬ 
ing  parameters:  w=2irX  6.683  47  X  109  sec  *,  g0(V7=0.1,  and 
tjh,=0.1.  These  numbers  are  easily  achievable  experimen¬ 
tally.  The  laser  pulse  width  tlp  is  chosen  to  be  10-12  sec  in 
order  to  satisfy  the  constraint  that  rL<i\ /&>  and  TL<il/g0M. 


PHYSICAL  REVIEW  A  69,  032308  (2004) 


FIG.  2.  Illustration  of  the  Bloch-Siegert  oscillation  (BSO):  (a) 
the  population  of  state  |l),  as  a  function  of  the  interaction  time  r, 
showing  the  BSO  superimposed  on  the  conventional  Rabi  oscilla¬ 
tion;  (b)  the  BSO  oscillation  (amplified  scale)  by  itself,  produced  by 
subtracting  the  Rabi  oscillation  from  the  plot  in  (a);  and  (c)  the  time 
dependence  of  the  Rabi  frequency.  Inset:  BSO  as  a  function  of  the 
absolute  phase  of  the  field. 

In  order  to  optimize  the  signal,  the  laser  Rabi  frequency  ilL 
is  chosen  to  be  such  that  £Iltl=  n r,  so  that  all  the  populations 
of  state  |l)  are  excited  to  state  |2)  at  the  end  of  the  pulse.  For 
the  cycling  transition  (1-2)  and  a  pulse  focused  to  an  area  of 
25  /rm2,  the  power  needed  for  achieving  this  Rabi  frequency 
is  1 .2  W,  which  is  achievable  experimentally.  After  the  laser 
pulse  is  turned  off,  the  fluorescence  is  collected  for  a  dura¬ 
tion  longer  than  the  spontaneous-decay  lifetime  (32  nsec)  of 
state  |2).  Under  this  condition,  our  simulation  verifies  that 
the  detector  signal  is  essentially  proportional  to  the  popula¬ 
tion  of  state  1 1 ),  as  given  by  Eq.  (13),  with  the  proportional¬ 
ity  constant  determined  by  the  efficiency  of  the  detection 
system.  If  106  atoms  are  used  (a  number  easily  achievable  in 
a  dipole  trap),  the  signal -to-noise  ratio  can  be  more  than  100 
for  the  parameters  considered  here,  assuming  a  detector  solid 
angle  of  0.1 7r  and  a  quantum  efficiency  of  0.8.  In  Fig.  2(a), 
we  have  shown  the  evolution  of  the  excited-state  population 
as  a  function  of  the  interaction  time  r  using  the  analytical 
expression  of  Eq.  (12).  Under  the  RWA,  this  curve  would 
represent  the  conventional  Rabi  oscillation.  However,  we  no¬ 
tice  here  some  additional  oscillations,  which  are  magnified 
and  shown  separately  in  Fig.  2(b),  produced  by  subtracting 
the  conventional  Rabi  oscillation  (sin2[gQ(7)T/2])  from 
Fig.  2(a).  That  is.  Fig.  2(b)  corresponds  to  what  we 
call  the  Bloch-Siegert  oscillation  (BSO),  given  by 
V  sin[^o(T)7‘]sin(2</'T).  The  dashed  curve  (c)  shows  the  time 
dependence  of  the  Rabi  frequency.  These  analytical  results 
agree  closely  with  the  results  obtained  via  direct  numerical 
integration  of  Eq.  (7).  Consider  next  a  situation  where  the 
interaction  time  t  is  fixed  so  that  we  are  at  the  peak  of  the 
BSO  envelope.  The  experiment  is  now  repeated  many  times, 
with  a  different  value  of  cf>  each  time.  The  corresponding 
population  of  |l)  is  given  by  r;  sin(2 <^>T)  and  is  plotted  as  a 
function  of  4>  in  the  inset  of  Fig.  2.  This  dependence  of  the 
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population  of  |l)  on  the  initial  phase  <f>  (and,  therefore,  on  the 
final  phase  cf>T )  makes  it  possible  to  measure  these  quantities. 

Note,  of  course,  that  the  speed  of  the  detection  system  is 
limited  fundamentally  by  the  spontaneous  decay  rate 
y-1  (—32  nsec  in  this  case)  of  state  \2).  As  such,  it  is  impos¬ 
sible  in  this  explicit  scheme  to  monitor  the  phase  of  the 
microwave  field  on  a  time  scale  shorter  than  its  period.  If  one 
were  interested  in  monitoring  the  phase  of  a  microwave  field 
of  a  lower  frequency  (so  that  aT1  S>  y-1),  it  would  be  possible 
to  track  the  phase  on  a  time  scale  much  shorter  than  its 
period.  One  possible  set  of  atomic  levels  that  can  be  used  for 
this  purpose  is  the  Zeeman  sublevels  (e.g.,  those  of  the 
5 2 3 4 5.S' | /2 :  /'  =  1  hyperfine  level  of  87Rb  atoms),  where  the  energy 
spacing  between  the  sublevels  can  be  tuned  by  a  dc  magnetic 
field  to  match  the  microwave  field  to  be  measured.  However, 
the  number  of  sublevels  that  get  coupled  is  typically  more 
than  2.  A  simple  extension  of  our  theoretical  analysis  shows 
that  the  signature  of  the  phase  of  the  microwave  field  still 
appears  in  the  population  of  any  of  these  levels  and  can  be 
used  to  measure  the  phase.  More  generally,  the  phase  signa¬ 
ture  is  likely  to  appear  in  the  population  of  the  atomic  levels, 
no  matter  how  many  levels  are  involved,  as  long  as  the  Rabi 
frequency  is  strong  enough  for  the  RWA  to  break  down. 

A  recent  experiment  by  Martinis  et  al.  [5]  is  an  example 
where  a  qubit  is  driven  very  fast.  In  this  experiment,  a  qubit 
is  made  using  the  two  states  of  a  current-biased  Josephson 
junction,  the  resonance  frequency  is  a>=6.9  GHz,  and  the 
Rabi  frequency  is  g=80  MHz.  If  this  experiment  is  carried 
out  without  keeping  track  of  the  phase  of  the  driving  field, 
the  degree  of  qubit  excitation  will  fluctuate  due  to  the  BSO, 
leading  to  an  error  which  is  of  the  order  of  g/w=0.01 — i.e., 


nearly  1%.  This  error  is  much  larger  than  the  permissible 
error  rate  of  10-6  for  an  error-correcting  quantum  computer 
that  would  consist  of  106 7 8 9 10 11  qubits  [11].  In  order  to  eliminate 
the  BSO-induced  error,  one  can  design  the  driving  system 
such  that  the  phase  is  measured,  e.g.,  by  using  an  auxiliary 
cluster  of  bits  located  close  to  the  qubit  of  interest,  at  the 
onset  of  the  qubit  excitation,  and  the  measured  value  of  the 
phase  is  used  to  determine  the  duration  of  the  excitation 
pulse,  in  order  to  ensure  the  desired  degree  of  excitation  of 
the  qubit  [14,15].  Finally,  we  point  out  that  by  making  use  of 
distant  entanglement,  the  BSO  process  may  enable  teleporta¬ 
tion  of  the  phase  of  a  field  that  is  encoded  in  the  atomic  state 
amplitude,  for  potential  applications  to  remote  frequency 
locking  [16-19]. 

In  conclusion,  we  have  shown  that  when  a  two-level 
atomic  system  is  driven  by  a  strong  periodic  field,  the  Rabi 
oscillation  is  accompanied  by  another  oscillation  at  twice  the 
transition  frequency,  and  this  oscillation  carries  information 
about  the  absolute  phase  of  the  driving  field.  One  can  detect 
this  phase  by  simply  measuring  only  the  population  of  the 
excited  state.  This  leads  to  a  phase-correlated  fluctuation  in 
the  excitation  of  a  qubit  and  violation  of  the  rule  that  the 
degree  of  excitation  depends  only  on  the  pulse  area.  We  have 
shown  how  the  resulting  error  may  be  significant  and  must 
be  controlled  for  low-energy  fast  qubit  operations. 
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Measuring  the  amplitude  and  absolute  (i.e.,  temporal  and  initial)  phase  of  a  monochromatic  microwave  field 
at  a  specific  point  of  space  and  time  has  many  potential  applications,  including  precise  qubit  rotations  and 
wavelength  quantum  teleportation.  Here  we  show  how  such  a  measurement  can  indeed  be  made  using  resonant 
atomic  probes  via  detection  of  incoherent  fluorescence  induced  by  a  laser  beam.  This  measurement  is  possible 
due  to  self-interference  effects  between  the  positive-  and  negative-frequency  components  of  the  field.  In  effect, 
the  small  cluster  of  atoms  here  act  as  a  highly  localized  pickup  coil,  and  the  fluorescence  channel  acts  as  a 
transmission  line. 


DOI:  10. 1103/PhysRevA.7 1.063408 

Measurement  of  the  amplitude  and  the  absolute  (i.e.,  tem¬ 
poral  and  initial)  phase  of  a  monochromatic  wave  is  chal¬ 
lenging  because  in  the  most  general  condition  the  spatial 
distribution  of  the  field  around  a  point  is  arbitrary.  Therefore, 
one  must  know  the  impedance  of  the  system  between  the 
point  of  interest  and  the  detector,  and  ensure  that  there  is  no 
interference  with  the  ambient  field.  It  is  recently  shown  in 
the  literature  that  the  absolute  phase  measurement  can  be 
used  for  accurate  qubit  rotations  [1-3]  and  quantum  wave¬ 
length  teleportation  [4-6]. 

Before  we  describe  the  physics  behind  this  process,  it  is 
instructive  to  define  precisely  what  we  mean  by  the  term 
“absolute  phase.”  Consider,  for  example,  a  microwave  field 
such  that  the  magnetic  field  at  a  position  R  is  given  by 
B(t)=B0  cos(u>t+  <p)x,  where  a>  is  the  frequency  of  the  field 
and  <[>  is  determined  simply  by  our  choice  of  the  origin  of 
time.  The  absolute  phase  is  the  sum  of  the  temporal  and 
initial  phases — i.e.,  a> t+<p.  In  order  to  illustrate  how  this 
phase  can  be  observed  directly,  consider  a  situation  where  a 
cluster  of  noninteracting  atoms  is  at  rest  at  the  same  location. 
For  simplicity,  we  assume  each  atom  to  be  an  ideal  two-level 
system  where  a  ground  state  |0)  is  coupled  to  an  excited  state 
|l)  by  this  field  B(f),  with  the  atom  initially  in  state  |0).  The 
Hamiltonian  for  this  interaction  is 

H  =  s(cr0  -  crz)/2  +  g(t)crx,  (1) 

where  g(t)  =  -go  cos (wt+cft),  g0  is  the  Rabi  frequency,  cr,  are 
the  Pauli  matrices,  and  the  driving  frequency  to=s  corre¬ 
sponds  to  resonant  excitation.  We  consider  g0  to  be  of  the 
form  £o(f)=£oM[l-exp(-t/TJvt,)]  with  a  switching  time  tsw 
relatively  slow  compared  to  other  time  scales  in  the 
system— i.e.,  tsw>  or'  and  g^'M. 

As  we  have  shown  before  [2,3],  without  the  rotating-wave 
approximation  (RWA)  and  to  the  lowest  order  in  rj 
=  (g0/4<u),  the  amplitudes  of  |0)  and  |l)  at  any  time  1  are  as 
follows: 

C0(f)  =  cos[go(f)t/2]  -  2rj%  sin[go(f)f/2],  (2) 


PACS  number(s):  32.80.Qk,  03.67.Hk,  03.67.Lx 

Cid)  =  i'e-,(“'+<W{sm[god)f/2]  +  2r) BX  cos[go(?)r/2]}, 

(3) 

where  S  =  (i72)exp[-/(2wf+2</>)]  and  g'Q(t) 
=  (l/f)/o£od'VC  =£0{l-d/ TsJ-1[l-exp(-f/ tJ]}.  If  we 
produce  this  excitation  using  a  7r/2  pulse  [i.e.,  £0(7)7=77/2] 
and  measure  the  population  of  state  |l)  after  the  excitation 
terminates  (at  t=f),  we  get  a  signal 

\C\{g'a{T),4>)\2  =  1/2+  77  sin[2(wr  +  <; b )].  (4) 

This  signal  contains  information  of  both  the  amplitude 
and  phase  of  the  field  B(t).  The  second  term  of  Eq.  (3)  is 
related  to  the  Bloch-Siegert  shift  [7,8],  and  we  have  called  it 
the  Bloch-Siegert  oscillation  (BSO)  [2,3].  It  is  attributable  to 
an  interference  between  the  so-called  corotating  and  counter¬ 
rotating  parts  of  the  oscillating  field,  with  the  atom  acting  as 
the  nonlinear  mixer.  For  r/=  0,  we  have  the  conventional 
Rabi  flopping  that  is  obtained  with  the  RWA.  For  a  stronger 
coupling  field,  where  the  RWA  is  not  valid,  the  second  term 
of  Eq.  (3)  becomes  important  [2,3],  and  the  population  will 
depend  now  both  on  the  Rabi  frequency  and  the  phase  of  the 
driving  field.  In  recent  years,  this  effect  has  also  been  ob¬ 
served  indirectly  using  ultrashort  optical  pulses  [9-11]  under 
the  name  of  carrier-wave  Rabi  flopping.  However,  to  the  best 
of  our  knowledge,  the  experiment  we  report  here  represents 
the  first  direct,  real-time  observation  of  this  effect. 

From  the  oscillation  observed,  one  can  infer  the  value  of 
2 (cuf+</>),  which  represents  the  absolute  phase  of  the  second 
harmonic.  This  is  equivalent  to  determine  the  absolute  phase 
of  the  fundamental  field,  ( u>t+cf> ),  modulo  77.  In  principle,  a 
simple  modification  of  the  experiment  can  be  used  to  elimi¬ 
nate  the  modulus  77  uncertainty.  Specifically,  if  one  applies  a 
dc  magnetic  field  parallel  to  the  rf  field,  it  leads  to  a  new 
oscillation  (in  the  population  of  either  level)  at  the  funda¬ 
mental  frequency,  with  exactly  the  same  phase  as  that  of  the 
driving  field.  In  the  experiment  described  here,  we  have  re¬ 
stricted  ourselves  to  the  case  of  determining  the  absolute 
phase  of  the  second  harmonic  only. 
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FIG.  1.  BSO  amplitude  versus  time  t  (in  units  of  gg1)  plots  for 
all  the  levels  of  a  three-level  system.  The  initial  densities  of  the 
levels  are  p00(f=0)  =  0.5,  pn(r=0)=0.3,  and  p22(?=0)  =  0.2,  the  Rabi 
frequency  g0=l,  and  the  resonant  frequencies  (o02=co2l  =  10. 

While  the  above  analytical  model  presented  here  is  based 
on  a  two  level  system,  practical  examples  of  which  are  pre¬ 
sented  in  Ref.  [2],  the  effect  is  more  generic,  and  is  present 
even  in  three-level  or  multilevel  systems.  In  particular,  we 
employed  a  three-level  system  to  observe  this  effect,  due 
primarily  to  practical  considerations.  The  specific  system 
used  consists  of  three  equally  spaced  Zeeman  sublevels  of 
87Rb  (5  2Sm\  F=  1:  mF=— 1,  0,  and  1,  denoted  as  states  |0), 
|l),  and  |2),  respectively),  where  the  degeneracy  can  be  lifted 
by  applying  an  external  bias  field.  We  have  performed  nu¬ 
merical  simulations  to  confirm  the  presence  of  the  BSO  sig¬ 
nature  in  the  population  dynamics  of  such  a  system  as  de¬ 
scribed  below. 

Consider  an  equally  spaced,  ladder-type  three-level  sys¬ 
tem  (|0),  1 1),  and  |2)).  The  transition  frequencies  for  |0)-|l) 
and  |l)-|2)  are  of  the  same  magnitude  e.  We  also  consider 
that  a  direct  transition  between  |0)  and  |2)  is  not  allowed. 
Now,  let  the  system  be  pumped  by  the  same  field  at  a  fre¬ 
quency  &>.  Consider  also  that  the  Rabi  frequency  for  the  |0)- 
|l)  transition  is  g0  and  that  for  1 1  )-|2)  is  also  g0.  Then,  the 
Hamiltonian  of  the  three-level  system  in  a  rotating  frame  can 
be  written  as 

H'  =-g0[l  +exp(-?2&>f-/20)](|O)<l|  +  |l)<2|)  +  c.c., 

(5) 

where  io=s.  The  amplitudes  of  the  three  levels  are  calculated 
numerically  by  solving  the  Schrodinger  equation  for  the 
above  Hamiltonian.  The  BSO  amplitudes  are  then  calculated 
by  subtracting  the  population  amplitude  of  each  level  with 
the  RWA  from  the  population  amplitude  without  the  RWA. 
The  BSO  oscillations  for  all  the  levels  of  such  a  system  are 
shown  in  Fig.  1. 

The  experimental  configuration,  illustrated  schematically 
in  Fig.  2,  uses  a  thermal,  effusive  atomic  beam.  The  rf  field  is 
applied  to  the  atoms  by  a  coil,  and  the  interaction  time  t  is 
set  by  the  time  of  flight  of  the  individual  atoms  in  the  rf  field 


FIG.  2.  Experimental  setup.  The  1-mm  cross  section  rubidium 
atomic  beam  passes  through  the  symmetry  axis  of  the  rf  coil  whose 
magnetic  field  is  along  the  beam.  The  rf  field  of  frequency  co  is  fed 
by  a  power  amplifier  connected  to  the  resonant  coil.  A  circularly 
polarized  probe  laser  beam  is  focused  down  to  30  /rm  in  diameter 
through  a  gap  in  the  middle  of  the  if  coil  and  perpendicularly  to  the 
atomic  beam.  The  atomic  fluorescence  is  collected  by  the  lens  and 
detected  by  an  avalanche  photodiode  (APD).  The  phase  signature 
appears  in  the  fluorescence  signal  encoded  in  an  oscillation  at  a 
frequency  2a>  due  to  the  Bloch-Siegert  oscillations.  In  the  picture,  S 
is  an  additional  phase  delay  due  to  the  APD  circuits  and  cabling. 
Inset:  Diagram  of  the  relevant  sublevels  of  the  D2  line  of  87Rb.  The 
numbers  on  the  left  represent  the  total  angular  momentum  of  the 
respective  levels.  The  strong  driving  rf  field  couples  to  the  ground- 
state  Zeeman  sublevels.  The  probe  beam  must  be  resonant  with  an 
appropriate  optical  transition  for  the  observation  of  the  phase- 
locked  signal,  as  discussed  in  the  main  text. 

before  they  are  probed  by  a  strongly  focused  and  circularly 
polarized  laser  beam.  The  rf  field  couples  the  sublevels  with 
|  Am  |  =  1,  as  detailed  in  the  inset  of  Fig.  2.  Optical  pumping  is 
employed  to  reduce  the  populations  of  states  |l)  and  |2)  com¬ 
pared  to  that  of  state  |0)  prior  to  the  interaction  with  the 
microwave  field. 

A  given  atom  interacts  with  the  rf  field  for  a  duration  t 
prior  to  excitation  by  the  probe  beam  that  couples  state  |0)  to 
an  excited  sublevel  in  5  2Py2-  The  rf  field  was  tuned  to 
0.5  MHz,  with  a  power  of  about  10  W,  corresponding  to  a 
Rabi  frequency  of  about  4  MHz  for  the  |0)— >|l)  as  well  as 
the  1 1)  — >■  |2)  transition.  The  probe  power  was  0.5  mW  fo¬ 
cused  to  a  spot  of  about  30  pirn  diameter,  giving  a  Rabi 
frequency  of  about  60  1.  where  F(6.06  MHz)  is  the  lifetime 
of  the  optical  transition.  The  average  atomic  speed  is 
500  m / s,  so  that  the  effective  pulse  width  of  the  probe,  tlp, 
is  about  60  ns,  which  satisfies  the  constraint  that  tlp  1  /  a>. 
Note  that  the  resolution  of  the  phase  measurement  is  essen¬ 
tially  given  by  the  ratio  of  min[riP,r_1]  and  1  /co,  and  can  be 
increased  further  by  making  the  probe  zone  shorter.  The  fluo¬ 
rescence  observed  under  this  condition  is  essentially  propor¬ 
tional  to  the  population  of  level  |0),  integrated  over  a  dura¬ 
tion  of  tlp,  which  corresponds  to  less  than  0.3  Rabi  period  of 
the  rf  driving  field  [for  g(tM/ (2  77)  =  4  MHz].  Within  a  Rabi 
oscillation  cycle,  the  BSO  signal  is  maximum  for  g0(T)r/2 
=  (2n+ l)7r/2,  where  n  =  0,l,2,...,  so  that  there  is  at  least 
one  maximum  of  the  BSO  signal  within  the  region  of  the 
probe. 
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FIG.  3.  Bloch-Siegert  oscillation  spectra,  rf  at  0.5  MHz  and 
Rabi  frequency  around  4  MHz.  (a)  Probe  beam  resonant  with  the 
5Si/2,  F=  1  <->553/2,  F'  =0  transition.  The  signal  appears  at  1  MHz 
with  a  linewidth  less  than  1  kHz  (resolution  limited  by  the  spectrum 
analyzer),  (b)  Probe  beam  blocked.  The  dip  structure  around 
100  kHz  is  an  artifact  due  to  the  amplifier  gain  curve.  Inset:  Spec¬ 
trum  for  same  configuration  and  rf  Rabi  frequency  around  10  MHz. 
Notice  the  2-MHz  harmonic  which  corresponds  to  the  higher-order 
BSO  at  4co. 

Note  that  atoms  with  different  velocities  have  different 
interaction  times  with  the  rf  field  and  produce  a  spread  in  the 
BSO  signal  amplitude  within  the  probe  region.  However,  the 
phase  of  the  BSO  signal  is  the  same  for  all  the  atoms,  since 
it  corresponds  to  the  value  of  (wt+c/))  at  the  time  and  loca¬ 
tion  of  interaction.  Thus,  there  is  no  washout  of  the  BSO 
signal  due  to  the  velocity  distribution  in  the  atomic  beam. 

Figure  3  shows  the  spectrum  of  the  observed  BSO  signal. 
In  Fig.  3(a),  we  show  that  the  BSO  stays  mainly  at  2a>.  When 
the  probe  beam  is  blocked,  there  is  no  signal  [Fig.  3(b)]. 
When  the  rf  intensity  is  increased  a  component  of  the  BSO  at 
4w  begins  to  develop,  as  predicted.  For  the  data  in  Fig.  4,  the 
second  harmonic  of  the  driving  field  is  used  to  trigger  a 
100-MHz  digital  oscilloscope  and  the  fluorescence  signal  is 
averaged  256  times.  When  the  probe  beam  is  tuned  to  the 
F=1<->F'=0  transition,  the  population  at  in  =  - 1  state  is 
probed.  When  the  probe  is  tuned  to  F=  1<->F'  =  1,  the  com¬ 
bined  populations  of  m  =  —  1  and  m  =  0  states  are  probed.  That 
results  in  an  effective  detection  of  the  complement  of  the 
population  of  m  =  1 .  On  the  other  hand,  when  the  probe  beam 
is  locked  to  the  F=1<->F'=2  transition,  all  three  Zeeman 
sublevels  of  F=  1  are  simultaneously  probed  and  the  phase 
information  is  not  clearly  present,  since  the  total  population 
of  level  F=  1  is  a  constant.  The  observed  residual  phase  in¬ 
formation  is  a  result  of  different  coupling  efficiencies  for 
each  of  the  three  ground  Zeeman  sublevels.  We  observed  that 
the  BSO  signal  amplitude  varies  as  a  function  of  an  external 
magnetic  field  applied  in  the  £  direction,  with  a  peak  corre¬ 
sponding  to  a  Zeeman  splitting  matching  the  applied  fre¬ 
quency  of  0.5  MHz. 

In  Fig.  5,  we  show  that  the  fluorescence  signal  is  phase 
locked  to  the  second  harmonic  of  the  driving  field.  First,  we 
placed  a  delay  line  of  0.4  ps  on  the  cable  of  the  reference 
field  used  to  trigger  the  oscilloscope  and  recorded  the  fluo¬ 


FIG.  4.  Time  dependence  of  the  fluorescence  signal  at  2 co  when 
the  probe  beam  is  resonant  to  different  excited  states.  The  lines  (a), 
(b),  and  the  noisy  line  (c)  correspond  to  the  probe  locked  to  the 
transitions  F=l— >F'=0,  F=\—*F'  =  \,  and  F=l  — >F'=2,  respec¬ 
tively,  of  the  5.S'i/2— > 5R3/2  transition  in  S7Rb. 

rescence  [Fig.  5(a)].  Then,  we  put  the  0.4-ps  delay  line  on 
the  BSO  signal  cable  and  recorded  the  fluorescence  [Fig. 
5(b)].  The  phase  difference  between  the  signals  recorded  in 
Figs.  5(a)  and  5(b)  is  approximately  0.8  ps,  as  expected  for  a 
phase  locked  fluorescence  signal.  The  data  presented  were 
for  the  probe  resonant  with  the  transition  F=  1  <->  F'  =  1 ,  but 
the  same  results  were  observed  for  F=1«->F'=0. 

To  summarize,  we  report  the  first  direct  observation  of  the 
absolute  phase  of  the  second  harmonic  of  an  oscillating  elec- 


Time(|4S) 

FIG.  5.  Demonstration  of  phase-locked  fluorescence.  T  is  the 
period  of  the  Bloch-Siegert  oscillation,  (a)  Population  vs  time  when 
a  0.4T  delay  line  was  inserted  in  the  reference  field  cable,  (b)  Popu¬ 
lation  vs  time  when  the  same  0.4T  delay  line  was  placed  in  the 
fluorescence  signal  cable.  The  figure  shows  that  signal  (b)  is  about 
0.8T  ahead  of  signal  (a),  confirming  that  the  atomic  fluorescence 
carries  phase  information  which  is  locked  to  the  absolute  rf  field 
phase.  The  solid  and  dashed  sinusoidal  smooth  curves  are  fittings  to 
the  experimental  data  and  were  used  for  period  and  delay 
determination. 
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tromagnetic  field  using  self-interference  in  an  atomic  reso¬ 
nance.  This  process  is  important  in  the  precision  of  quantum 
bit  rotations  at  a  high  speed.  The  knowledge  of  the  absolute 
phase  of  a  rf  field  at  a  particular  point  of  space  may  also  be 
useful  for  single-atom  quantum  optics  experiments.  For  ex¬ 
ample,  an  extension  of  this  concept  may  possibly  be  used  to 
teleport  the  wavelength  of  an  oscillator,  given  the  presence 
of  degenerate  distant  entanglement,  even  in  the  presence  of 
unknown  fluctuations  in  the  intervening  medium  [4-6,12]. 
Finally,  this  localized  absolute  phase  detector  may  prove  use¬ 


ful  in  mapping  of  radio-frequency  fields  in  microcircuits.  Al¬ 
though  a  particular  alkali-metal  atom  was  used  in  the  present 
experiment,  the  mechanism  is  robust  and  could  be  observed 
in  virtually  any  atomic  or  molecular  species. 
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